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□ □ 


Experimentally  determined  wall  pressure  distributions,  local  surface 
shear  stresses  and  their  directions,  and  detailed  turbulent  boundary-layer  tra- 
verses in  near  zero  heat  transfer  conditions,  are  presented  through  skewed 
shock /boundary-layer  interaction  regions  generated  by  a wedge  standing 
normal  to  a test  wall.  The  mainstream  Mach  numbers  were  2 and  4,  while  the 
Reynolds  number  based  on  the  undisturbed  test  boundary-layer  thickness  of 
0.2-in.,  growing  along  the  nozzle  sidewall  of  the  NAE  5 X 5-in.  blowdown 
wind  tunnel,  was  ~2  X 105 . 

Tangential  air  injection  at  a jet  exit  Mach  number  of  3 was  then 
introduced  into  the  3D  shock  separated  Mach  2 boundary  layer,  to  control 
the  separation.  The  optimum  direction  of  blowing  was  found  to  be  along  a 
line  somewhere  between  the  deflected  surface  of  the  wedge  and  the  line  of 
the  oblique  shock  wave. 


SOMMAIRE 


Les  pressions  statiques  sur  la  paroi  d’essai,  les  frottements  de  surface 
avec  les  directions,  et  quelques  profils  de  vitesse,  sont  presentes  aux  inter- 
actions en  trois  dimensions,  dans  lesquelles  une  onde  de  choc  oblique 
(produite  par  un  diedre  debout  sur  la  paroi)  a rencontre  une  couche  limite 
turbulente.  Les  nombres  de  Mach  a l’ecoulement  non-perturbe  etaient  2 et  4, 
en  meme  temps  le  nombre  de  Reynolds  etait  2 X 105 , calcule  avec  l’epaisseur 
non-perturbe  de  la  couche  limite  (—  0.2-inch).  La  couche  limite  existait  sur 
la  paroi  du  col  convergent-divergent  dans  la  soufflerie  a jet  d’air  compnme 
“NAE  5 X 5-inch”. 

L’injection  d’air  tangent  a la  paroi,  au  decollement  en  trois  dimen- 
sions de  la  couche  limite  (produit  par  une  onde  de  choc  oblique  dans  un 
ecoulement  non-perturbe  au  nombre  de  Mach  2)  afin  de  controler  le  decolle- 
ment. La  meilleure  direction  de  la  paroi-jet  etait  trouvee  le  long  d’une 
ligne  entre  la  surface  du  diedre  et  la  projection  sur  la  paroi  de  l’onde  de  choc 
oblique. 
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THREE  DIMENSIONAL  SWEPT  SHOCK/TURBULENT 
BOUNDARY-LAYER  SEPARATIONS  WITH  CONTROL  BY  AIR  INJECTION 

1.0  INTRODUCTION 

1.1  Objective 

The  goal  of  the  present  study  is  the  control  by  tangential  air  injection  of  a three-dimensional 
turbulent  boundary-layer  separation  in  a Mach  2 mainstream  flow,  that  is  caused  by  a swept  shock  wave 
standing  normal  to  a flat  test  wall.  The  injectant  air  issues  at  a nominal  Mach  number  of  3,  with  an  ex- 
cess of  momentum  in  the  jet  somewhat  more  than  the  momentum  deficit  of  the  undisturbed  boundary 
layer,  0 [ , . The  Reynolds  number  based  upon  0 , j at  the  slot  location  and  the  mainstream  flow  condi- 
tions is  nominally  14,000.  The  stagnation  temperatures  of  the  jet  and  mainstream  are  the  same  while 
the  direction  of  the  jet  may  be  altered  relative  to  the  three-dimensional  separation  line. 

A secondary  objective  is  the  detailed  measurement  of  the  three-dimensional  interaction  between 
the  swept  shock  wave  and  the  turbulent  boundary  layer  at  Mach  numbers  of  2 and  4,  at  shock  pressure 
ratios  near  and  greater  than  the  pressure  ratio  for  incipient  three-dimensional  separation. 

The  course  of  presentation  of  the  results  will  proceed  from  a discussion  of  the  swept  inter- 
action phenomena  to  the  mechanism  of  the  boundary-layer  control  by  jet  blowing. 

1.2  Overview 

The  potential  usefulness  of  artificially  controlling  the  behaviour  of  the  boundary  layer  has 
been  recognized  from  the  earliest  days  of  research  in  viscous  flow.  It  was  in  1904  that  Ludwig  Prandtl, 
to  support  his  discovery  of  the  boundary  layer,  described  several  experiments  to  show  the  great  influ- 
ence that  suction  boundary-layer  control  could  exert  upon  the  flow  about  a circular  cylinder. 

We  may  conveniently  classify  the  subject  of  boundary-layer  control  under  three  sub-headings 
concerned  respectively  with:  the  prevention  of  separation  and  the  increase  of  lift  in  subsonic  flow; 
control  for  low  drag;  and  the  problem  of  shock-induced  separation  (and  heat  transfer)  at  supersonic 
speeds  and  its  prevention  by  boundary-layer  control  with  cool  gas  injection.  Most  practical  schemes  so 
far  have  been  concerned  with  suppressing  separation  over  leading-edge  und  trailing-edge  flaps  or  wing 
top  surfaces  in  subsonic  flow  — the  ‘Buccaneer’  aircraft  is  one  prominent  example  — although  a well- 
known  application  in  supersonic  flow  has  been  the  suction  control  of  boundary  layers  in  supersonic 
air  inlets,  upstream  and  at  the  location  of  the  intake  throat  normal  shock  wave. 

Boundary-layer  control  (blc)  may  be  accomplished  with  either  ‘active’  or  ‘passive’  systems. 
‘Active’  or  powered  systems  may  be  moving  surfaces,  discrete  or  distributed  suction,  or  blowing  at  some 
angle  to  the  surface.  ‘Passive’  systems  include  slats,  ots,  fences  and  small  plate-type  vortex  generators, 
where  some  part  of  the  higher  energy  local  mainstream  is  redirected  by  the  passive  device  to  mix  with 
the  sluggish  viscous  flow,  and  fill  out,  on  average,  the  mean  velocity  profile. 

On  aerofoils,  boundary- layer  control  by  blowing  may  be  employed  to  give  high  lift,  but  usu- 
ally only  for  short  periods  of  time  to  improve  take-off  and  landing  performance.  The  injected  air,  in 
re-energizing  the  turbulent  boundary  layers,  delays  the  stall,  and  enables  high  angles  of  attack  and  large 
flap  deflexions  to  be  used.  In  contradistinction,  to  provide  low  drag  at  cruise,  distributed  suction  over 
the  entire  aerofoil  surface  may  be  installed  as  in  the  speculative  laminar  flow  wing.  Such  low  drag  aero- 
foils are  usually  provided  with  suction  slits  for  artificially  maintaining  a laminar  flow  across  the  aerofoil 
chord.  The  suction  ensures  a thin  stable  boundary  layer  while  laminarization  removes  the  relatively 
high  skin  friction  drag  associated  with  transition  to  a turbulent  boundary  layer.  These  two  categories 
of  boundary-layer  control,  specifically  for  either  high  lift  or  for  >rw  drag,  are  quite  distinct;  and  a sys- 
tem for  one  will  not  in  general  be  applicable  to  the  other.  The  research  and  full  scale  flight  work  on 
suction  and  blowing  blc  on  subsonic  aircraft  up  to  1961  is  fully  and  lucidly  described  in  Lachmann’s 
two  volumes1  on  the  subject. 
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It  is  worthy  of  note  that  the  present  and  future  high  cost  of  aviation  fuel,  combined  with  the 
interest  in  flying  large  commercial  and  cargo  transport  aircraft  at  speeds  approaching  Mach  1,  are  causing 
a reappraisal  of  employing  ‘active’  boundary-layer  control  during  cruise.  One  recent  study2  suggests 
that  the  concept  of  the  laminar  flow  wing  should  be  resurrected  to  be  coupled  with  the  new  supercritical 
wing  technology;  and  that  tangential  blowing  should  be  used  for  boundary-layer  control  at  several  sta- 
tions along  the  fuselage  particularly  near  the  wing  fuselage  junction,  where  separations  usually  exist. 

The  penalty  for  the  blc  installation  might  well  be  offset,  if  considerable  drag  savings  were  to  result.  An 
earlier  investigation  by  Whitcomb3  at  NASA  Langley,  yielded  results  that  a ‘p ass've’  device  incorporated 
into  a supercritical  wing  could  be  beneficial  for  controlling  shock/boundary-layer  interactions.  A wing 
section  with  a pronounced  flattened  upper  surface  was  equipped  with  a fixed  slot  near  the  trailing-edge, 
at  about  the  80  percent  chord  station.  Supersonic  flow  was  maintained  across  the  top  surface  of  the 
aerofoil  as  far  as  the  slot,  where  the  flow  then  diffused  through  a normal  shock  wave,  to  recover  to  a 
level  of  positive  (subsonic)  pressure  coefficient  at  the  trailing-edge.  The  slot  directed  high  pressure  air 
from  the  lower  surface  to  mix  with  the  low-energy  air  of  the  top  surface  boundary  laye-  The  advan- 
tages claimed  for  the  supercritical  wing  with  boundary-layer  control  included  a significant  delay  in  the 
drag  rise  corresponding  with  a notable  increase  in  cruise  Mach  number.  A noteworthy  example  given 
in  Reference  4 was  that  the  incorporation  of  such  a wing  into  the  Boeing  707  aircraft  would  potentially 
improve  the  cruise  efficiency  by  15  percent  and  permit  cruise  Mach  numbers  in  excess  of  0.9  to  be  used 
without  requiring  any  increase  in  engine  thrust. 

The  advent  of  the  supersonic  commercial  jet  transport  aircraft,  with  a continuing  emphasis 
on  reducing  direct  operating  costs,  is  encouraging  the  further  development  of  efficient  boundary-layer 
control  schemes  in  the  diffusion,  from  supersonic  to  subsonic  speeds,  of  the  air  to  the  propulsion  en- 
gines. Control  by  bleed  or  suction  is  currently  favoured,  but  there  is  a good  potential  for  control  by 
blowing5 . The  advantages  of  re-energizing  the  momentum  deficient  fluid  on  the  intake  wetted  areas 
compared  with  removing  it  by  bleed  are  several.  Only  a small  volume  is  required  for  the  ducting  of 
high  pressure  air,  relatively  large  pressure  ratios  are  available  by  tapping  a compressor  stage,  duct  losses 
are  low  and  there  is  no  loss  of  the  intake  capture  flow  upstream  of  the  compressor  face  as  there  is  in 
the  case  of  bleed. 

Boundary-layer  control  has  been  recognized  as  an  acceptable  and  practical  means  of  achieving 
performance  gains  on  aerodynamic  configurations  especially  in  subsonic  flow.  However,  there  is  a dearth 
of  information  on  the  control  of  boundary-layer  separation  in  supersonic  flow,  especially  in  a typical 
three-dimensional  flow  configuration,  and  it  is  hoped  that  the  present  study  will  at  least  rectify  a small 
part  of  this  deficiency. 

1.3  Previous  Work  on  Supersonic  B'owing 

To  provide  a background  to  the  present  study,  we  shall  present  a brief  appraisal  of  work-to- 
date  on  wall  jets  mixing  with  boundary  layers  in  adverse  pressure  gradients.  Such  work  has  almost 
exclusively  dealt  with  quasi/two-dimensional  flows  and  two-dimensional  separations,  whereas  on  air- 
craft we  are  usually  dealing  with  complex  three-dimensional  flow-fields.  The  work  herein  attempts  to 
rectify  previous  omissions  by  dealing  with  blowing  boundary-layer  control  of  shock-induced  three- 
dimensional  separation,  and  is  one  topic  in  the  continuing  investigation  of  three-dimensional  boundary 
layers  by  the  senior  authors  of  Reference  6. 

The  performance  of  a supersonic  speed  blowing  jet  is  based  on  the  premise  that  for  a given 
configuration,  a prescribed  value  of  jet  excess  momentum  flux  compared  with  the  momentum  defi- 
ciency of  the  boundary  layer  to  be  re-energized,  will  result  in  virtually  the  same  result,  regardless  of 
the  mass  flow  and  jet  velocity  chosen  for  the  excess  momentum.  The  excess  momentum  flux  at  the 
jet  exit,  is  defined  as  [riij  (Uj  - uE  )]cxit,  in  order  to  correspond  with  the  momentum  deficiency  of  a 
boundary  layer.  The  definitions  are  virtually  the  same  except  for  sign  for  pE  uE  2 0 1 p = / pu(uE  - u)dy 
= /(uF.  - u)drh.  Clearly,  unless  the  jet  velocity  exceeds  the  local  velocity  outside  the  boundary  layer, 
there  will  be  a decrease  in  the  boundary-layer  momentum,  the  reverse  of  the  desired  goal  of  blc.  The 
consequence  will  be  an  increase  in  form  drag,  although  the  skin  friction  may  be  reduced.  We  should 
note  that  the  majority  of  recent  publications  on  tangential  injection7-14  with  supersonic  mainstreams, 
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both  analytical  and  experimental,  have  been  to  determine  the  magnitudes  of  reduction  in  skin  friction 
that  might  be  possible  with  low  energy  air  injection  from  discrete  slots.  Such  reduction  is  achieved  by 
altering,  distorting  and  thickening  the  boundary-layer  profile  such  that  the  gradient  of  mean  velocity 
at  the  surface,  which  is  a direct  measure  of  local  skin  friction,  is  reduced.  A thorough  review  of  con- 
cepts for  reducing  surface  friction  and  reducing  wall  temperature  at  high  supersonic  speeds  is  presented 
in  Reference  15.  It  is  claimed  that  for  effectiveness  of  cooling  the  surface,  a slow  mixing  rate  between 
the  original  boundary  layer  and  the  injected  fluid  is  desirable;  and  the  effects  of  compressibility  are 
such  that  when  compared  with  equivalent  incompressible  flow  injectant  velocity  to  freestream  ratios, 
the  mixing  in  compressible  flow  takes  a much  longer  downstream  distance. 

We  will  emphasize  that  we  shall  be  dealing  with  an  injectant  fluid  velocity  greater  than  tne 
external  stream,  where  the  injected  gas  behaves  as  a wall  jet  and  an  increase  in  local  skin  friction  occurs. 
Such  an  increase  in  skin  friction  is  obviously  desirable  if  control  of  boundary -layer  separation  is  to  be 
achieved  in  strong  adverse  pressure  gradient  conditions. 

Let  us  now  define  the  terms  ‘jet  peak’  and  ‘wake  trough’,  for  a high  energy  wall  jet  mixing 
with  a boundary  layer,  where  the  velocity  of  the  injected  air  is  larger  than  the  external  stream.  Let  us 
refer  to  Figure  1.  At  the  exit  plane  of  the  blowing  nozzle,  the  jet  velocity  is  zero  at  the  nozzle  lip  and 
at  the  lower  wall  because  of  the  internal  boundary  layers  on  the  nozzle  wetted  areas.  Immediately 
downstream  of  the  lip,  a mixing  region  forms  between  the  jet  and  the  original  boundary  layer.  The 
mean  velocity  profile  downstream  of  the  jet  exit  has  consequently  a maximum  and  a minimum:  the 
maximum  is  called  the  ‘jet  peak’  and  the  minimum,  the  ‘wake  trough’.  These  maxim';  and  minima 
vary  with  distance  downstream  of  the  jet  exit  to  eventually  merge 

Stratford16  in  1953  suggested  that  the  concept  of  high  energy  air  injection  for  boundary-layer 
control  in  adverse  pressure  gradients,  showed  considerable  promise  although  more  basic  experimental 
evidence  was  desirable  It  was  proposed  that,  ideally,  instantaneous  mixing  between  the  jet  and  the 
boundary  layer  would  permit  an  adverse  pressure  gradient  to  be  applied  immediately  without  fear  of 
separating  the  boundary  layer.  For  this  ideal  conception  it  was  reckoned  the  excess  momentum  of  the 
injection  fluid  should  equal  the  momentum  deficiency  of  the  original  boundary-layer  fluid. 

It  was  pointed  out,  however,  that  for  a practical  flow  the  mixing  of  the  two  streams  required  a 
finite  distance;  also,  that  a new  wall  boundary  layer  grew  under  the  jet.  Neither  the  extent  of  the  mixing 
region  nor  the  history  of  the  jet  peak  downstream  of  the  injection  point  appeared  calculable  at  that 
time,  and  for  a quantitative  estimate,  experimental  results  were  required.  It  was  thought  that  about 
2.5  times  the  ideal  quantity  of  jet  excess  momentum  would  be  necessary  for  total  removal  of  the 
boundary-layer  momentum  defect  and  to  give  a reserve  of  total  pressure  in  the  jet  for  any  subsequent 
control  needs. 

It  was  mentioned  that  the  idealized  treatment  might  no  longer  be  appropriate  for  flow  with 
lower  energy  injection,  that  is,  when  the  jet  velocity  was  only  slightly  more  than  the  local  mainstream 
velocity.  If  an  adverse  pressure  gradient  were  applied  too  quickly  before  adequate  mixing  had  occurred, 
then  the  low  energy  air  in  the  mixing  region  would  be  slowed  considerably  or  reversed  in  flow  direction. 

Some  experiments  were  performed  by  staff  of  the  National  Physical  Laboratory  at  about  the 
same  time17'18;  these  were  to  estimate  the  use  of  blowing  to  reduce  the  fall  in  control  effectiveness 
associated  with  separation  of  the  boundary  layer  that  was  caused  by  an  upper  surface  normal  shock 
wave.  These  tests  were  made  at  transonic  speeds.  It  was  stated  that  air  injection  reduced  the  effects 
of  shock-induced  separation  of  the  boundary  layer. 

A blowing  slot  was  located  it  about  2/3  chord  in  the  convex  top  surface  of  a half  section 
aerofoil  (chord  of  6-in.).  The  boundary  layer  on  the  aerofoil  was  turbulent  and  local  jet  conditions 
were  sonic  at  the  blowing  slot.  Two  sizes  of  blowing  slot  were  used:  these  were  0.009-in.  and  0.005-in. 
in  height.  The  former  height  was  about  0.1  of  the  thickness  of  the  original  boundary  layer  at  the 
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blowing  slot.  Blowing  stagnation  pressures  up  to  10  times  the  mainstream  stagnation  pressure  were 
used,  with  blowing  momentum  coefficients,  CM,  up  to  0.05*. 

The  authors1 7 suggested  that  the  injected  air  altered  the  flow  development  by  preventing 
the  original  boundary  layer  from  growing  normally.  They  proposed  that  the  effects  of  blowing  with 
the  narrower  slot  were  the  same  as  those  obtained  with  the  wider  one,  but  smaller  quantities  of  air 
were  required  for  the  narrower  slot  to  produce  similar  effects  on  the  flow  development. 

It  was  also  reported  that  with  the  slot  immersed  in  the  dead  air  or  separation  bubble  region 
of  the  normal  shock /boundary-layer  interaction,  the  injection  of  small  quantities  of  air  produced  a 
large  effect  cn  the  separation.  In  fact,  for  a blowing  total  pressure  of  about  2.3  times  the  mainstream 
total  pressure,  the  boundary  layer  was  able  to  withstand  the  full  pressure  rise  through  a normal  shock 
without  separating.  From  the  figures  in  Reference  17,  we  infer  that  the  static  pressure  ratio  across  the 
shock  was  of  the  order  of  two.  With  the  slot  upstream  of  the  separation  and  the  mixing  region  length 
about  15  original  boundary-layer  thicknesses  (estimates  from  the  data  given)  larger  quantities  of  air 
were  required  to  produce  an  appreciable  effect  on  flow  development. 

Although  in  Reference  17  it  was  proposed  that  the  overall  effect  of  injecting  into  the  sepa- 
rated region  of  the  boundary  layer  was  beneficial,  no  pitot  traverses  were  taken  to  indicate  any  flow 
irregularities  (with  consequent  losses)  in  the  mixing  region  — the  flow  characteristics  described  in 
Reference  17  were  obtained  by  considering  the  static  pressure  distribution  on  the  aerofoil  surface  and 
Schlieren  flow  visualization. 

The  potential  advantages  of  slot  blowing  on  aerofoils  claimed  by  the  NPL  some  twenty  years 
ago,  have  been  proved  correct  by  Peake  and  Yoshihara19  in  some  recent  tests  at  NAE,  to  determine 
the  transonic  lift  augmentation  of  two-dimensional  supercritical  aerofoils  by  means  of  aft  camber,  slot 
blowing  and  jet  flaps,  in  high  Reynolds  number  flow.  Figure  2(a)  illustrates  the  type  of  supercritical 
aerofoil  examined,  having  a ‘peaky’  nose,  an  aft  camber  of  about  20-degrees  at  the  trailing-edge  and  a 
chord  of  15-in.  The  aerofoil  was  provided  with  a sonic  blowing  slot  on  the  upper  surface  at  the  80  per- 
cent chord  station,  as  we  see  on  Figure  2(b).  The  0.045-in.  slot  gap  was  approximately  25  percent  of 
the  local  boundary-layer  thickness. 

The  usefulness  of  slot  blowing  is  seen  in  the  chord-wise  pressure  distributions  on  Figure  3,  in 
some  tests  at  a mainstream  Mach  number  of  0.85,  at  a Reynolds  number  of  30  million  based  on  the 
15-in.  chord.  With  no  boundary-layer  control,  the  aerofoil  is  close  to  a condition  of  maximum  lift, 
and  at  or  on  the  verge  of  buffet  onset.  The  upper  surface  shock  wave  is  positioned  near  mid-chord, 
with  shock-induced  separated  flow  extending  from  the  shock  to  the  trailing-edge.  The  flow  has  attained 
a local  Mach  number  of  about  1.5  upstream  of  the  shock,  and  the  static  pressure  ratio  across  the  shock 
wave  is  grea'  er  than  two. 

On  applying  slot  blowing,  with  a momentum  coefficient,  CM , equal  to  0.01  (that  is  a jet 
excess  momentum  about  equal  to  the  momentum  deficiency  of  the  boundary  layer)  the  upper  surface 
pressure  distribution  changes  dramatically.  The  overall  effect  of  the  slot  blowing  is  particularly  bene- 
ficial for  it  increases  the  effective  aft  camber,  as  evidenced  by  the  increased  plateau  loading  and  the 
large  rearward  displacement  of  the  upper  surface  shock  wave,  by  almost  30  percent  of  the  wing  chord. 
In  spite  of  a local  separation  bubble  downstream  of  the  jet  slot  position,  the  flow  almost  recovers  to  a 
positive  pressure  coefficient  at  the  trailing-edge.  Even  allowing  a penalty  on  drag  equal  to  CM  for  sup- 
plying the  blowing  air,  we  note  a better  than  30  count  reduction  in  profile  drag,  and  an  increase  in  lift 
coefficient  if  0.27. 

Bradshaw  and  Gee20  explain  some  interesting  flow  features  in  their  work  on  turbulent  wall 
jets  with  and  without  an  external  stream.  Although  their  mainstream  and  jet  velocities  were  subsonic, 


Cfl  is  defined  as  rrijUj/q^  Sw , where  ir»j  and  Uj  are  the  mass  flow  rate  and  jet  velocity  (assuming 
isentropic  expansion  to  mainstream  static  pressure),  and  q„  and  Sw  are  the  mainstream  dynamic 
pressure  and  wing  area  respectively. 
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their  results  bear  quoting  because  of  their  observations  of  jet  excess  momentum  required.  The  main- 
stream speed  was  varied  up  to  120-ft/s  and  the  jet  velocities  were  650-ft/s  and  350-ft/s  for  slot  sizes 
of  height  0.018-in.  and  0.040-in.  respectively.  A calculation  suggests  that  the  slot  heights  were  approx- 
imately 0.01  and  0.10  of  the  original  boundary-layer  thicknesses  at  the  blowing  slot.  Pressure  traverses 
were  made  from  which  the  mean  velocity  distributions  were  obtained.  Distances  downstream  from  the 
injection  point  of  up  to  15  original  boundary-layer  thicknesses  were  investigated  for  flow  with  and  with- 
out adverse  pressure  gradient.  No  comment  was  offered  on  any  effect  of  slot  size. 

It  was  suggested  from  the  experimental  measurements,  that  for  a tangential  wall  jet  to  ‘absoio' 
the  initial  boundary  layer  in  zero  pressure  gradient,  the  jet  excess  momentum  must  be  about  five  times 
the  original  boundary-layer  momentum  deficiency.  The  main  part  of  the  original  boundary-layer  ve- 
locity defect  was  absorbed  quickly  and  the  jet  peak  was  still  pronounced  at  about  15  original  boundary  - 
layer  thicknesses  from  the  blowing  slot.  A non-zero  value  of  the  turbulent  component  of  the  shear 
stress  at  the  jet  peak  was  also  found. 

With  positive  pressure  gradient,  a value  of  jet  excess  momentum  of  about  twice  that  of  the 
original  boundary-layer  momentum  deficiency  appeared  adequate  to  prevent  separation,  but  not 
enough  to  entrain  completely  the  original  boundary  layer.  The  value  of  jet  excess  momentum  was 
suggested  as  being  lower  for  the  flow  with  adverse  pressure  gradient  since  the  high  speed  jet  suffered 
relatively  less  retardation  than  the  mainstream  — also  with  adverse  pressure  gradient,  the  skin  friction 
coefficient  was  less.  Stratford’s16  proposal  that  the  jet  excess  momentum  should  be  about  2.5  times 
that  of  the  original  boundary- layer  momentum  deficiency  was  thus  reasonable  when  compared  with 
the  results  of  Bradshaw  and  Gee. 

Neale  and  Stratford21  related  an  attempt  to  re-energize  an  intake  sidewall  boundary  layer 
whose  flow  direction  had  been  subject  to  three-dimensional  shock  wave/boundary-layer  interactions*. 

A convergent-divergent  blowing  nozzle  of  slot  height  0.5-in.  (about  0.4  of  a boundary-layer 
thickness)  on  each  sidewall,  injected  air  tangential  to  the  wall  into  the  three-dimensional  separated 
flow  region.  The  jet  flow  expanded  in  the  nozzle  to  match  the  local  condition  of  static  pressure  in  the 
rig.  Pressure  traverses  were  made  further  downstream  at  the  throat  where  the  mainstream  Mach  number 
was  1.75.  Two  locations  of  jet  nozzle  exit  were  tried:  the  first  was  at  about  3.5  and  the  second  at 
roughly  two  original  boundary-layer  thicknesses  upstream  of  the  throat  traverse  position. 

For  the  first  jet  nozzle  exit  position,  the  wake  trough  total  pressure  at  the  throat  traverse 
location  was  slightly  more  than  the  local  value  of  static  pressure  and  equalled  0.15  of  rig  inlet  total 
pressure.  The  Mach  number  in  the  wake  trough  was  0.27,  so  that  the  flow  in  the  mixing  region  was 
still  considerably  retarded.  The  jet  peak  Mach  number  was  not  known,  but  was  greater  than  2.5. 

Blowing  nearer  to  the  throat  traverse  position  produced  a wake  trough  total  pressure  of  0.06 
tunnel  entry  pressure,  which  was  equal  to  the  local  s*°tk  pressure.  The  minimum  Mach  number  in  the 
wake  was  measured  as  zero  and  it  was  therefore  likely  that  the  wake  flow  was  reversed.  As  before,  the 
jet  peak  Mach  number  was  greater  than  2.5. 

For  both  blowing  positions,  therefore,  the  throat  traverses  revealed  little  or  no  mixing 
between  the  jet  and  the  original  boundary  layer.  With  the  application  of  further  static  pressure  rise 
downstream  of  the  injection  point,  the  intake  flow  system  broke  down  entirely. 

The  extent  of  mixing  of  about  two  to  3.5  boundary -layer  thicknesses  was  evidently  insuffi- 
cient to  destroy  the  three-dimensional  separated  flow  of  the  original  bound ^-y  layer.  The  inference 


The  intake  was  an  all-internal  compression  type  designed  to  operate  at  Mach  3.  The  upper  and 
lower  surfaces  consisted  of  a series  of  ramps,  the  oblique  shock  waves  from  which  focussed  at 
the  intake  throat.  These  shock  waves  interacted  with  the  boundary  layers  on  the  flat  sidewalls 
to  cause  significant  three-dimensionality  and  eventual  separation  of  the  sidewall  flows  near  the 
intake  throat. 
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to  be  obtained  from  the  work  of  Neale  and  Stratford2 1 was  that  blowing  into  a three-dimensional 
separated  boundary  layer  did  not  produce  the  favourable  effects  as  reported  in  References  17  and  18 
for  two-dimensional  interactions.  Reference  21  concluded  that  a fundamental  comprehension  of  the 
injection  phenomena  in  two  and  three  dimensions  was  necessary  before  re-energization  techniques 
could  be  justified  in  intake  applications. 

The  foundation  stones  were  hence  laid  for  a new  experiment  to  elucidate  the  qualitative  and 
quantitative  aspects  of  a supersonic  wall  jet  mixing  with  a retarded  and  supersonic,  turbulent  boundary 
layer.  Peake,  in  Reference  5,  attempted  such  an  investigation  in  a quasi/two-dimensional  flow,  with  a 
mainstream  Mach  number  of  1.8  at  the  jet  exit.  The  Reynolds  number  at  the  slot  based  upon  the 
momentum  deficit  thickness  of  the  undisturbed  boundary  layer  was  about  7000.  High  energy  air,  of 
variable  total  pressure  and  momentum  was  injected  into  the  boundary  layer,  tangential  to  the  wall,  to 
minimize  mixing  losses.  The  injection  nozzle  was  a convergent/divergent  design,  with  straight  walls 
downstream  of  the  injection  throat,  yielding  a Mach  number  of  roughly  2.4  at  design  pressure  ratio. 

The  slot  height  of  0.1-in.  was  less  than  one  half  the  original  boundary-layer  thickness.  The  combined 
boundary  layer  and  the  jet  with  and  without  the  imposition  of  adverse  pressure  gradient  was  studied 
optically  via  Schlieren  records  and  by  pressure  measurements.  An  incident  oblique  shock  wave  and 
its  reflection  created  the  region  of  adverse  pressure  gradient,  whose  strength  could  be  varied  by  altering 
the  flow  deflection  angle  of  a wedge  in  the  tunnel  ceiling.  Provision  was  made  for  changing  the  injection 
position  relative  to  the  position  of  the  applied  adverse  pressure  gradient.  Pitot  traverses,  static  pressure 
traverses,  wall  static  pressure  measurements  and  Schlieren  photographs  enabled  the  properties  of  the 
new  wall  boundary  layer  under  the  jet  and  the  mixing  between  the  jet  and  the  original  boundary  layer 
to  be  measured.  A check  was  made  on  the  quasi/two-dimensionality  of  the  flow  by  establishing  that 
the  tunnel  centre-line  flow  quantities  with  and  without  injection,  more  or  less  satisfied  the  momentum 
integral  equation. 

For  a physical  evaluation  of  the  flow,  it  was  found  possible  to  subdivide  the  flow  of  the 
jet/boundary-layer  combination  into  (a)  the  new  wall  boundary  layer  under  the  jet  and  a consideration 
of  jet  peak  characteristics;  and  (b),  the  wake  or  mixing  region  between  the  jet  and  the  original  boundary 
layer.  To  prevent  flow  reversal  in  an  adverse  pressure  gradient,  in  either  of  (a)  or  (b)  above,  it  was 
necessary  to  choose  a length  downstream  for  mixing  so  that  (i)  the  wake  trough  total  pressure  was 
greater  than  the  local  static  pressure  from  the  applied  pressure  gradient  — otherwise  reversal  of  the 
wake  flow  resulted;  and  (ii),  the  falloff  in  jet  peak  total  pressure  was  not  severe.  The  indication  from 
Reference  5 was  that  at  best,  the  wake  trough  total  pressure  approached  a constant  value  of  about  0.6 
of  the  mainstream  total  pressure.  The  initial  rate  of  mixing  between  the  jet  and  the  original  boundary 
layer  within  the  range  tested  was  much  the  same  for  all  injection  stagnation  pressures  and  jet  excess 
momentum  rates. 

Tests  showed  that  with  a jet  excess  momentum  slightly  more  than  the  momentum  deficiency 
>_f  the  original  boundary  layer  and  a mixing  region  extent  of  six  to  ten  reference  boundary-layer  thick- 
nesses, or  approximately  100  undisturbed  momentum  deficit  thicknesses,  control  (of  shock-induced 
separation)  could  be  maintained  at  least  up  to  static  pressure  ratios  of  three.  (The  Mach  number  of  the 
mainstream  approaching  the  jet  slot  was  1.8.)  The  injection  of  air  actually  at  the  location  of  the  shock- 
induced  separation  produced  an  overall  beneficial  effect  in  that  the  flow  satisfactorily  negotiated  the 
shock  system:  although  it  was  probable  that  the  wake  flow  or  the  wall  boundary  layer  was  reversed. 

The  wall  boundary  layer  separation,  if  present,  was  always  short  in  extent  and  was  followed  by  a rapid 
reattachment.  The  use  of  injection  decreased  the  upstream  distance  over  which  the  pressure  rise  oc- 
curred and  increased  the  pressure  rise  obtained  compared  with  the  no-injection  case.  This  latter  feature 
would  appear  plausible  because  the  jet  Mach  number  is  greater  than  the  freestream  value;  while  the 
former  observation  results  from  the  increased  supersonic  velocity  region  close  to  the  wall. 

An  empirically  derived  equation  was  presented  that  correlated  the  minimum  velocity  of  the 
wake  with  distance  downstream  of  the  injection  point.  The  velocity  defect  in  the  wake  was  shown  to 
be  roughly  proportional  to  the  inverse  square  root  of  the  distance  downstream  of  the  injection  point. 
With  reasonable  satisfaction  of  the  momentum  integral  equation  as  was  demonstrated  in  Reference  5, 
the  equation  might  thus  be  utilized  as  a first  step  in  predicting  the  characteristics  of  a jet  mixing  with 
either  a boundary  layer  or  a wake. 
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A subsequent  quasi/two-dimensional  investigation  on  blowing  in  supersonic  flow  to  alleviate 
shock-induced  separation  was  reported  by  Howell  and  Tatro22 . They  developed  an  analysis  based  on 
previously  observed  experimental  separation  pressure  ratios  and  the  boundary-layer  momentum  defi- 
ciency, to  predict  criteria  for  control.  In  this  author’s  view,  the  analysis  contained  several  anomalies, 
a discussion  of  which  is  presented  in  Reference  129.  A written  critique  from  this  author  to  Howell 
and  Tatro  elucidated  no  response.  They  obtained  results  at  mainstream  Mach  numbers  of  3.2  and  4.2 
and  at  a flow  deflection  angle  of  10-degrees.  No  momentum  quantities  were  quoted  in  terms  of  the 
undisturbed  boundary-layer  deficiency,  and  this  author  would  disagree  with  Howell  and  Tatro’s  obser- 
vations from  wall  static  pressure  distributions  of  ‘separation’  and  ‘no  separation’. 

Alzner  and  Zakkay10  investigated  the  shock/boundary-layer  interaction  in  axisymmetric 
flow  at  Mach  6,  where  a 10-degree  axisymmetric  wedge-shaped  ramp  was  fitted  around  an  axisymmetric 
centre-body,  which  was  used  for  the  test  boundary  layer.  Pitot  and  static  pressure  profiles,  stagnation 
temperature  profiles,  wall  static  pressure  and  heat  transfer  measurements,  and  optical  data  were  obtained 
in  this  comprehensive  investigation.  Both  air  and  hydrogen  were  used  as  injeetant  fluids  with  choked 
jet  exit  Mach  numbers.  The  slot  height  varied  between  one  and  three  percent  of  the  undisturbed 
boundary-layer  thickness,  5 = 0.85-in. 

An  inspection  of  the  results  reveals  that  separation  was  not  destroyed,  although  the  prime 
objective  of  these  tests  was  to  realise  the  potential  benefit  of  wall  cooling  from  the  jet. 

In  summary,  the  work  of  Reference  5 appeared  to  offer  a sound  base  for  determining  whether 
the  control  characteristics  would  be  carried  over  to  a three-dimensional  swept  shock/turbulent  boundary- 
layer  interaction  at  Mach  2.  Such  a flow  situation  is  pertinent  to  swept  wing  flows,  to  flows  in  wing/body 
junctions,  and  to  sidewall  boundary-layer  interactions  in  quasi/two-dimensional  supersonic  inlets  currently 
installed  on  present-day  aircraft.  In  addition,  no  work  had  been  reported  on  the  effectiveness  or  other- 
wise of  changing  the  wall  jet  direction  relative  to  the  three-dimensional  separation  and/or  the  mainstream 
direction.  An  investigation  of  the  three-dimensional  turbulent  flow  downstream  of  swept  slot  injection 
in  hypersonic  Mach  6 flow,  in  zero  pressure  gradient,  was  reported  recently,  however,  by  Hefner  et  al23, 
in  which  the  film-cooling  effectiveness  of  the  injeetant  fluid  was  the  essential  target  of  the  study,  as 
well  as  the  reduction  in  skin  friction  because  of  the  low  (sonic)  Mach  number  of  injection.  Due  to  the 
latter,  the  jet  flow  injected  tangential  to  the  wall,  at  22.5  or  45  degrees  to  the  mainstream  direction, 
also  appeared  to  not  penetrate  more  than  typically  20  slot  heights  downstream  before  returning  to  the 
mainstream  direction. 

The  present  investigation  was  begun  in  1966  and  the  5-in.  X 5-in.  blowdown  wind  tunnel 
modifications  completed  in  1968.  The  tests  were  completed  during  1974  and  the  review  of  previous 
work  compiled  at  that  time. 

1.4  Two-Dimensional  Interactions  Between  Shock  Waves  and 
Supersonic  Turbulent  Boundary  Layers 

Most  of  the  investigations  of  shock  wave/turbulent  boundary-layer  interactions  during  the 
last  twenty  years  have  focussed  on  the  largely  qualitative,  and  sometimes  quantitative  mean  flow  aspects 
of  nominal  two-dimensional  flows.  There  is  no  firm  indication  that  there  is  quantitative  similarity  be- 
tween two  and  three-dimensional  interactions30  and  we  should  preferably  refer  to  the  nominal-2D 
interaction  as  the  limiting  case  of  the  more  general,  and  certainly  more  often  found,  3D  swept  inter- 
action. 

To  illustrate  the  conceptual  difference  between  the  nominal-2D  shock  wave/boundary-layer 
interaction  and  the  3D  swept  shock  interaction,  let  us  turn  to  Figure  4.  We  see  that  both  types  of 
interactions  may  be  produced  by  a wedge  in  a rectangular  duct.  Such  a geometry  is  representative  of 
the  initial  compression  surface  of  a supersonic  inlet  currently  employed  on  military  and  commercial 
supersonic  aircraft,  where  the  3D  interaction  is  formed  on  the  flat  sidewall  of  the  inlet.  The  boundary 
layer  on  the  sidewall  skews  in  the  interaction  region  and  tends  to  flow  in  a downstream  direction  close 
to  the  line  of  intersection  of  the  shock  with  the  sidewall. 


. 
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In  deference  to  the  historical  development  of  the  subject  of  interactions  between  shock  waves 
and  turbulent  boundary  layers,  however,  we  comment  briefly  on  some  of  the  important  findings  by  in- 
vestigators of  the  nominal-2D  flows. 

It  is  well-known  that  the  presence  of  a boundary  layer  developing  along  a solid  surface  causes 
an  oblique  shock  wave  to  reflect  as  a complex  pattern  of  waves,  rather  than  as  a single  reflected  wave  in 
inviscid  flow.  The  boundary  layer  is  essentially  composed  of  two  strata,  an  outer  supersonic  region  that 
merges  with  the  mainstream,  and  a subsonic  inner  layer  adjacent  to  the  wall.  The  subsonic  nature  of  the 
latter  implies  that  the  pressure  is  continuous  there,  and  that  disturbances  may  be  transmitted  in  both 
upstream  and  downstream  directions.  In  the  shock  interaction  region,  large  shear  stress  gradients  per- 
pendicular to  the  wall  are  generated  in  the  inner  layer,  which  assist  in  dragging  low  energy  air  downstream 
into  areas  of  high  pressure.  The  outer  supersonic  layer,  on  the  other  hand,  may  be  thought  of  as  more 
nearly  inviscid,  and  may  be  qualitatively  described  by  a pattern  of  waves.  The  incoming  shock  wave  is 
refracted  as  it  passes  through  the  rotational  supersonic  part  of  the  boundary  layer,  becoming  negligible 
in  strength  as  the  sonic  line  is  approached.  The  refraction  mechanism,  however,  involves  both  a change 
in  strength  of  the  incoming  wave  and  the  generation  of  outgoing  secondary  waves,  the  latter  merging  to 
produce  the  ‘reflected’  shock.  In  addition,  other  compression  waves  are  generated  by  the  outwardly 
concave  shaping  of  the  sonic  line,  due  to  thickening  of  the  subsonic  layer  in  the  adverse  pressure  gradient. 

If  the  inner  layer  does  not  separate  from  the  wall,  the  streamwise  extent  of  the  interaction  is 
typically  two  or  three  undisturbed  thicknesses  for  turbulent  boundary  layers.  Once  separation  com- 
mences, the  growth  of  the  separated  flow  region  causes  the  outer  flow  pattern  co  depart  increasingly 
from  that  which  would  be  seen  in  inviscid  flow.  With  large  regions  of  separation,  there  is  strong  coupling 
between  the  viscous  and  inviscid  flow  regimes  and  it  is  debatable  whether  the  boundary  layer  and  main- 
tream  may  still  be  considered  as  separate  entities,  with  a suitably  defined  ‘edge’  to  the  viscous  flow. 

The  qualitative  aspects  and  quantitative  mean-flow  measurements  in  nominal  two-dimensional 
shock  wave/turbulent  boundary-layer  interactions  have  been  assessed  by  many  investigators  in  the  last 
two  decades;  notable  authors  in  the  subject  being  Gadd,  Holder,  Pearcey  and  Green  in  the  United  Kingdom, 
and  Bogdonoff  and  his  co-workers  at  Princeton  University  in  the  United  States.  References  24  to  29, 
for  instance,  provide  the  essential  constituents  to  our  understanding  of  the  complexities  of  shock/boundary- 
layer  interactions.  Green30  has  compiled  a well-organized  review  of  the  work  in  two-dimensional  shock/ 
boundary-layer  interactions  investigated  prior  to  1969,  which  almost  exclusively  have  involved  mean-flow 
measurements.  To  advance  the  state  of  knowledge,  measurements  of  the  turbulence  structure  in  the 
interaction  region  are  required,  and  some  recent  experimental  investigations31’32  have  attempted  this. 

The  experiment  of  Grande  and  Oates3 1 was  performed  using  constant  temperature  hot  wire 
anemometry  and  miniature  pressure  transducers  in  nominal  two-dimensional  shock/ boundary-layer 
interactions  in  a Mach  2.4  mainstream,  with  combined  incident  and  reflected  shock  strengths  high 
enough  to  cause  separation.  The  Reynolds  number  based  on  the  undisturbed  boundary-layer  thickness 
of  0.2-in.  was  equal  to  R6  of  3.4  X 104 . Grande  and  Oates  found  that  the  turbulence  levels  in  the 
boundary  layer  increased  as  a result  of  the  interaction  and  increased  with  increasing  shock  wave  strength. 

Rose32  made  a similar  experiment  to  Grande  and  Oates,  measuring  the  rms  intensities  of  fluctuating 
mass  flux  and  total  temperature,  and  their  correlation  coefficients,  through  a Mach  4,  axisymmetric 
shock  wave/turbulent  boundary-layer  interaction.  The  Reynolds  number  based  on  the  test  boundary- 
layer  thickness  of  0.2-in.  was  R6  = 8.7  X 104,  again,  a relatively  low  dynamic  pressure  flow  to  prevent 
destruction  of  the  hot  wires.  In  agreement  with  Reference  31,  Rose  also  found  that  the  intensity  of 
the  boundary-layer  turbulence  increased  as  a result  of  the  interaction,  while  the  fluctuations  of  mass 
flux  and  total  temperature  remained  highly  correlated  over  most  of  the  boundary  layer  throughout 
the  interaction.  Rose  also  noticed  that  the  intensity  of  the  turbulence  increased  at  certain  frequencies 
more  than  others  at  the  same  h/6  above  the  wall.  On  the  basis  of  the  above  findings,  we  might  expect 
that  a shear  layer  passing  through  a Prandtl-Meyer  expansion  fan  would  provide  large  decreases  in 
turbulence  intensity;  and  the  experiments  of  Lewis  and  Behrens33  appear  to  support  this  contention. 

t j 


In  closing  this  brief  discussion  on  nominal  two-dimensional  shock/boundary -layer  interactions, 
we  should  note  that  the  mean  value  of  surface  shear  stress  generated  by  a strong  adverse  pressure  gradient 


in  supersonic  flow  has  not  been  readily  amenable  to  calculation.  For  example,  Peake  et  al34  demon- 
strated that  the  surface  shear  stress  beneath  an  axisymmetric  decelerating  flow  from  Mach  4 to  Mach  2 
could  not  be  duplicated  by  present  day  calculation  methods.  Bradshaw35  used  the  results  of  Refer- 
ences 34  and  36  to  demonstrate  that  the  consistent  discrepancies  between  calculation  and  experiment 
were  due  to  unexpectedly  large  effects  of  mean  dilatation  (or  compression)  on  the  Reynolds  stresses. 
Bradshaw  found  that  agreement  between  calculation  and  experiment  was  substantially  improved  by 
representing  the  effects  of  mean  dilatation  in  the  same  way  as  that  used  satisfactorily  to  represent  the 
effects  of  extra  rates  of  strain  in  low-speed  flows.  (We  recollect  that  according  to  Morkovin’s  hypo- 
thesis37 the  effects  of  density  fluctuations  produced  by  the  turbulence  itself  are  relatively  small  at 
speeds  less  than  hypersonic.) 

Finally,  the  separation  point  and  the  length  of  separation  in  supersonic  flow  have  been  shown 
dependent3  8 on  the  manner  of  diagnosis,  whether  by  a mean  flow  or  fluctuating  flow  technique.  Pitot 
pressure  surveys  and  surface  skin  friction  data  indicated  the  furthest  downstream  separation  and  the 
shortest  separation  bubble  in  an  axisymmetric  flow  with  adverse  pressure  gradient  at  Mach  7.2,  6 = 1-in. 
and  R6  = 0.2  X 106.  Both  surface  pressure  data  (between  inflection  points)  and  oil  flow  results  indi- 
cated a slightly  earlier  separation,  while  power  spectra  and  correlation  measurements  (using  thin  platinum 
films  flush  mounted  with  the  wall  on  pyrex  glass)  indicated  a still  earlier  separation.  These  results  imply, 
therefore,  that  the  onset  and  reattachment  locations  of  separation  in  quasi- 2D  flow  are  intermittent. 

1.5  Three-Dimensional  Separations  in  Swept  Interactions 

Although  the  fundamentals  of  skewed  boundary  layers  and  three-dimensional  separations 
were  reported  almost  twenty  years  ago,  it  is  perhaps  opportune  to  reflect  upon  the  physical  aspects  of 
three-dimensional  separations.  We  must  desist  from  thinking  only  in  terms  of  two-dimensional  criteria 
for  separation  and  endeavour  to  promote  a suitable  and  general  concept  of  flow  separation.  Ideas  based 
on  two-dimensional  flow  where  separation  lines  must  be  perpendicular  to  the  external  flow  direction 
are  of  limited  value  in  synthesizing  flow  separations  in  three  dimensions  For  the  sake  of  simplicity, 
the  calculation  of  boundary  layers  still  tends  to  exclude  the  third  space  dimension,  leading  to  overt 
emphasis  by  experimenters  on  flows  with  similar  (it  is  hoped)  two-dimensional  restrictions  in  order  to 
allow  comparison  with  theory. 

Considering  the  subject  of  three-dimensional  viscous  flows*  in  its  historical  context,  Sears  in 
1948,  utilized  the  very  important  concept  of  ‘limiting  streamlines’  in  a paper39  discussing  the  laminar 
boundary  layer  on  an  infinite  yawed  cylinder.  Each  streamline  near  the  surface,  in  the  inner  collinear 
part  of  the  boundary  layer,  is  one  of  a continuous  pattern  of  streamlines  that  approaches  the  surface 
at  an  attachment  point  or  line  and  leaves  the  surface  at  some  point  or  line  of  separation.  (We  remember 
that  an  aerofoil  with  attached  flow  over  its  surface  still  suffers  separation  at  a cusped  trailing-edge. ) 
Clearly,  if  the  distance  of  a streamline  above  the  surface  is  infinitesimal,  its  direction  is  everywhere 
tangential  to  the  vanishing  fluid  velocity  as  the  surface  is  approached.  It  is  then  identifiable  as  a limiting 
streamline,  or,  in  other  words,  it  possesses  the  same  direction  as  a skin  friction  line  or  surface  shear 
stress  trajectory.  These  limiting  streamlines  at  the  base  of  the  boundary  layer  must  not  be  confused 
with  the  streamlines  of  the  external  inviscid  flow  over  the  surface,  and  they  will  in  general  follow  path, 
that  are  different  in  direction  to  those  of  the  external  (inviscid)  streamlines.  The  terms:  limiting  stream- 
line and  surface  shear  stress  trajectory  are  hence  synonymous  for  all  practical  purposes. 

Sears,  moreover,  included  a sketch39  which  showed  the  convergence  of  limiting  streamlines 
towards  a three-dimensional  separation  line  (and  from  either  side  of  it)  although  he  did  not  call  it  by 
name;  and  later,  Wild40  substantiated  Sears’  physical  interpretation  of  separation,  in  his  investigation 
of  the  laminar  boundary  layer  on  yawed  ‘infinite’  wings.  In  subscribing  to  the  Cornell  University  work, 
Eichelbrenner  and  Oudart41-42  proposed  that  a line  of  three-dimensional  separation  was  the  envelope 
of  the  limiting  streamlines  (conversely,  a line  of  attachment  is  the  origin  of  divergent  limiting  stream- 
lines) and  this  was  also  supported  by  Susan  Brown43  in  a mathematical  treatise  of  the  Royal  Society. 

term  viscous  region  or  viscous  flow  is  used  instead  of  boundary  layer  because  both  separated 
unseparated  flows  may  be  discussed  together,  without  ambiguity 
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Lighthill44,  on  the  other  hand,  disputed  the  term  ‘envelope’,  in  spite  of  flow  visualization  experiments 
apparently  supporting  the  previous  authors.  Coalescence  of  the  limiting  streamlines  is  perhaps  an  accept- 
able description  to  avoid  argument. 

It  remained  for  Maskell45 , however,  to  lay  the  ground  rules  for  the  precise  interpretation  of 
the  physical  composition  of  three-dimensional  viscous  flows  (including  the  effects  after  the  flow  leaves 
the  surface)  by  demonstrating  that  the  limiting  streamlines  provide  a ‘three-dimensional  skeleton  struc- 
ture of  the  viscous  flow’.  Provided  a flow  visualization  indicator  is  very  thin  on  a wind  tunnel  model, 
it  has  been  shown46  that  the  direction  indicated  by  the  streaks  on  the  surface  is  very  close  to  the  direc- 
tion of  surface  shear  and  the  limiting  direction  of  the  velocity  when  the  indicator  is  absent.  Thus  the 
many  surface  flow  visualization  techniques  yield  a powerful  means  of  diagnosing  the  qualitative  features 
of  three-dimensional  viscous  flows. 

Two  types  of  separation  were  conceived  by  Maskell45 : a bubble  type  and  a free  shear  layer 
type,  both  of  which  may  exist  together.  The  latter  would  appear  to  be  the  more  common  flow  found 
in  practice.  Fluid  accumulating  at  a three-dimensional  separation  line,  which  in  general  is  set  obliquely 
to  the  direction  of  the  external  and  essentially  inviscid  flow  field,  leaves  the  separation  line  as  a free 
shear  layer  and  rolls  up  in  the  process  of  passing  downstream.  The  shear  stress  is  usually  finite  along 
such  a three-dimensional  separation  line,  although  it  may  be  zero  at  a singular  point.  High  local  suction 
pressures  are  frequently  induced  on  the  surface  beneath  the  resulting  viscous  vortices,  close  to  which 
are  noted  inflection  points  ir.  the  limiting  streamline  patterns.  Some  examples  are  the  delta  wing  flow 
(the  ONERA  and  RAE  work  being  the  most  extensive  — see  the  subject  list  in  Reference  47)  and 
cones47'48'49,  the  flows  on  long  pointed  slender  bodies  at  incidence,  those  about  upswept  fuselages50, 
and  those  about  bluff  protuberances47 . In  the  bluff  or  slender  protuberance  case,  for  example,  which 
applies  to  many  junction  problems,  the  rolled-up  shear  layer  departing  from  the  three-dimensional 
separation  line  occurs  within  the  depth  of  the  oncoming  boundary  layer,  but  the  changes  to  the  pressure 
distribution  are  still  significant. 

In  the  swept  shock  wave/turbulent  boundary-layer  interaction,  the  flow  in  the  plane  normal 
to  both  the  shock  wave  and  the  wall  is  frequently  quoted  as  resembling  a two-dimensional  flow  through 
a normal  shock.  In  Figure  5(a),  we  show  schematically,  the  result  of  a weak  shock/boundary-layer 
interaction,  where  although  the  limiting  streamlines  are  deflected  substantially  more  than  the  inviscid 
streamlines,  there  is  no  coalescence  of  them,  and  no  three-dimensional  separation.  The  small  sketch  on 
the  right-hand  side  is  a view  along  the  shock  wave,  showing  a pattern  and  interaction  analogous  to  that 
for  a weak  normal  shock  wave  in  nominal  two-dimensional  flow. 

With  increasing  shock  strength,  from  small  to  larger  static  pressure  ratios,  there  reaches  a stage 
in  the  flow  development  where  three-dimensional  separation  just  begins,  and  important  qualitative  flow 
changes  start  to  occur.  When  the  pressure  rise  is  sufficiently  strong  to  cause  coalescence  of  the  limiting 
streamlines,  the  picture  might  be  that  shown  on  Figure  5(b).  The  view  along  the  shock  wave  might 
provide  us  with  an  ‘open’  or  ‘closed’  separation  bubble,  with  a consequent  lambda  foot  to  the  shock 
and  a vortical  slip  line  passing  downstream  from  the  triple  point  of  shock  intersection.  In  the  left-hand 
sketch,  the  ‘open’  or  ‘closed’  bubble  commences  at  the  3D  separation  line  upstream  of  the  inviscid  shock 
wave  position.  Limiting  streamlines  emanating  from  a reattachment  line  downstream  of  the  shock  pass 
through  the  line  of  the  inviscid  shock  to  run  with  the  separation  line.  The  upstream  limiting  streamlines 
run  with  the  separation  line  from  ahead  of  the  shock*.  A more  detailed  cross-section  of  the  flow  when 
viewed  along  the  shock  might  follow  the  picture  produced  by  Seddon52  in  2D  flow,  as  we  see  on 
Figure  6. 

The  shear  stress  is  finite  along  the  separation  line,  and  the  flow  is  generally  in  the  direction  of 
the  shock  wave.  In  the  immediate  neighbourhood  of  the  junction  of  a shock  generator  with  the  wall, 
where  the  interaction  is  spreading  quickly  the  scale  of  the  separation  increases  with  distance  outwards 


* 


If  the  bubble  is  ‘open’  the  flow  leaves  the  surface  along  the  separation  line  as  a free  shear  layer  and 
rolls  up  to  form  a vinous  vortex. 
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along  the  shock,  tending  it  is  felt,  from  a conical  type  flow  to  one  of  cylindrical  similitude  further  away 
from  the  shock  generator.  The  reattachment  stream  surface  passes  over  the  stream  surface  springing 
from  the  separation  line. 


When  an  incident  oblique  shock  wave  collided  with  a turbulent  boundary  layer  in  nominal 
two-dimensional  flow,  Kuehn5 1 observed  that  incipient  separation  coincided  with  the  appearance  of 
three  inflection  points  in  the  surface  static  pressure  distribution  (Fig.  8)  although  Seddon52  did  not 
notice  the  same  feature  in  a 20  normal  shock/turbulent  b jundary -layer  interaction.  On  the  other 
hand,  Stanbrook53  did  observe  the  phenomenon  in  his  3D  swept  shock  investigation.  Clearly,  if  the 
boundary  layer  is  relatively  thin,  and  the  pressure  rise  is  very  steep,  then  the  physical  extent  of  the 
interaction  region  will  be  short,  and  inflection  points  will  not  be  discerned5 . 

Green30  in  other  incident  oblique  shock  wave  experiments  in  2D  flow,  detected  abrupt 
changes  in  the  reflected  wave  pattern,  that  he  linked  with  the  appearance  of  separation.  He  noted  that 
there  was  a sudden  doubling  of  strength  of  the  outgoing  shock  and  as  a consequence  an  almost  eight 
fold  increase  in  stagnation  pressure  loss  through  it. 

We  have  seen  that  the  characteristic  feature  of  a 3D  separation  line  is  that  the  limiting  stream- 
lines or  surface  shear  stress  trajectories  approach  the  line  asymptotically  from  both  sides,  or  as  appear 
in  experiment,  they  coalesce  to  form  an  envelope  (Fig.  7(a)).  McCabe54,  utilizes  the  envelope  concept 
in  his  study  of  a three-dimensional  sidewall  interaction  by  means  of  a shock  wave  generator  standing 
normal  to  the  test  wall.  He  obtained  the  shear  stress  trajectories  in  and  about  the  3D  interaction  region 
by  overall  coating  of  the  surface  with  an  oil  base  emulsion.  Coagulation  of  large  amounts  of  the  emul- 
sion occurred  in  the  vicinity  of  separation  due  to  local  boundary-layer  thickening  and  the  consequent 
reduced  shear  stresses  at  the  wall.  Jn  principle,  it  would  appear  difficult  to  define  incipient  separation 
from  such  a technique,  unless  the  mixture  can  be  constituted  to  provide  feather-fine  streaks55 . 

Lowrie56  used  a criterion  devised  by  Rogers  and  Hall57  for  separations  on  swept  wings 
(Fig.  7(b))  that  defined  incipient  separation  as  occurring  when  the  upstream  surface  shear  stress  trajec- 
tories converge  into  one  line  that  is  parallel  to  the  surface  shear  stress  trajectories  immediately  just 
downstream  of  it,  and  that  still  further  downstream,  they  turn  progressively  away  from  it.  Again,  this 
criterion  is  probably  adequate,  provided  the  flow  visualization  pattern  is  sufficiently  detailed.  It  is  for 
this  reason  that  the  oil  dot  technique58  is  to  be  generally  preferred  for  limiting  streamline  definition, 
for  no  messy  coagulation  of  surface  oil  occurs. 

Stanbrook5  3 used  the  observation  to  define  incipient  separation,  that  at  their  maximum  deflec- 
tion, the  limiting  streamlines  run  parallel  to  the  shock  wave. 


So  far,  it  has  not  been  proven  which  criterion  is  best,  but  it  is  probable  that  the  pressure  ratios 
obtained  at  incipient  separation  from  any  of  the  techniques,  would  differ  little.  We  should  note,  however, 
that  for  a general  three-dimensional  flow  case,  say  on  a body,  the  commencement  of  a three-dimensional 
separation  is  difficult  to  pin-point.  Figure  9,  for  example,  shows  a typical  result  of  the  surface  shear 
stress  trajectories  obtained  with  the  oil  dot  technique,  where  there  are  three-dimensional  separations 


and  attachments  on  the  upswept  afterbody  of  a cargo  transport  aircraft  model.  The  test  was  conducted 
at  a Mach  number  of  0.72,  and  at  a Reynolds  number  of  15  million  based  on  the  afterbody  length  of 
1.5-ft.  with  turbulent  boundary  layers  present.  The  convergence  of  the  limiting  streamlines  implying 
local  thickening  of  the  boundary  layers  followed  by  the  formation  of  an  ‘envelope’  is  well  illustrated  — 
but  it  should  be  emphasized  that  we  can  only  guess  the  absolute  point  of  envelope  commencement. 
These  remarks,  perhaps,  provoke  some  thoughts  about  how  three-dimensional  separation  should  be 
defined  in  a general  case.  The  features  of  boundary-layer  convergence  and  hence  thickening  into  lobes 
of  flow,  resulting  finally  in  shedding  of  fluid  from  the  region  of  the  separation  line  might  better  be 
perceived  together  as  a region  of  ‘three-dimensional  separation’.  We  will  reiterate  that  along  the  en- 
velope, the  resultant  surface  shear  stress  is  non-zero  in  general,  whereas  for  a nominal  two-dimensional 
separation,  the  converse  is  true.  We  should  also  note47  that  the  effects  of  three-dimensional  separation 
may  be  beneficial  and  not  necessarily  catastrophic  as  is  usual  in  the  two-dimensional  case.  It  appears 
that  the  shedding  of  the  vortices  in  the  three-dimensional  situation  is  quasi-steady,  and  unlike  the  ran- 
dom, pulsating  separated  flow  in  nominal  two-dimensional  cases  such  as  stall  of  high  aspect  ratio  wings. 
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Indeed,  the  very  interaction  problem  under  investigation  herein  has  been  used  to  direct  the  boundary 
layers  on  intake  sidewalls  into  the  ramp  throat  bleed  on  the  compression  surface  perpendicular  to  the 
sidewall. 


The  length  of  the  swept  interaction  in  the  streamwise  sense  depends  upon  the  pressure  ratio 
across  the  shock  wave,  the  nature  of  the  disturbance  causing  it,  and  the  upstream  Mach  number,  Reynolds 
number  and  oncoming  boundary-layer  thickness,  60.  The  latter  has  been  commonly  used  as  a scaling 
parameter  for  streamwise  distance  through  the  interaction,  utilizing  the  inviscid  shock  wave  position  as 
origin.  Equally  well,  the  streamwise  displacement  thickness  Sj  *,  would  also  be  suitable.  The  stream- 
wise  momentum  thickness,  j , is  also  a characteristic  parameter  for  scaling  purposes,  for  it  changes 
rapidly  with  distance  through  the  interaction,  but  varies  only  slowly  in  regions  of  constant  pressure 
upstream  and  downstream  of  the  interaction. 

Korkegi64  has  found  that  a simple  correlation  exists  for  incipient  turbulent  boundary-layer 
separation  due  to  a skewed  shock  wave.  He  took  the  incipient  separation  data  of  McCabe54  and  Lowrie56 , 
along  with  the  corresponding  shock  wave  angles,  Mach  numbers  and  pressure  rises,  and  concluded  that 
the  flow  deflection  angle  at  incipient  separation  6wi,  was  inversely  proportional  to  the  upstream  main- 
stream Mach  number,  where 


0.30 


Korkegi  also  associated  incipient  separation  occurring  at  an  approximate  constant  pressure  rise,  Pj/p,  = 1.5, 
independent  of  Mach  number,  at  least  over  the  range  of  3.44  > M,  > 1.96.  tWe  note  that  such  behaviour 
is  contrary  to  that  for  cwo-dimensional  turbulent  boundary-layer  separation,  where  and  Pj/pj  both  in- 
crease with  increase  of  M, .) 

By  comparison,  Goldberg6  5 has  reported  that  the  angle  for  3D  incipient  separation  at  M j =5.9 
is  virtually  insensitive  to  Reynolds  number  based  on  boundary -layer  thickness  in  the  range  10s  < R6  < 106, 
and  is  about  9°.  This  is  not  consistent  with  an  extrapolation  of  Korkegi’s  correlation. 

1.6  Some  Previous  Experiments  in  Swept  Shock/Boundary-Layer  Interactions 

As  we  noted  previously,  there  is  little  quantitative  understanding  of  interactions  between  swept 
shock  waves  and  turbulent  boundary  layers.  Stalker59  and  Stan  brook 5 3 introduced  and  discussed  some 
qualitative  features  of  the  swept  interaction  region,  while  McCabe54  at  Manchester  University  and  Lowrie56 
at  Cambridge  provided  some  quantitative  insight  into  the  cross  flows  developing  in  the  interaction.  Much 
of  the  existing  analysis  has  been  developed  for  so-called  ‘infinite’  swept  flows,  and  so  is  of  a quasi  two- 
dimensional  nature;  while  most  experiments  have  been  conducted  using  the  nozzle  boundary  layers  in 
relatively  small  supersonic  wind  tunnels. 

In  1960,  Stalker59  published  the  results  of  some  experiments  using  the  sidewall  turbulent 
boundary  layer  of  a 6 X C-in.  intermittent  type  blowdown  wind  tunnel.  He  investigated  the  separation 
regions  produced  by  swept  forward  and  swept  rearward-facing  steps,  and  in  a flow  generated  by  an 
oblique  shock  wave  standing  normal  to  the  test  surface  Stalker  demonstrated  his  results  to  be  broadly 
consistent  with  simple  extensions  of  two-dimensional  flow  analyses.  In  his  forward-facing  step  experi- 
ment, the  plateau  pressure  rise  ahead  of  the  step  decreased  as  the  angle  of  sweep  of  the  step  increased. 

The  2D  free-interaction*  analysis  of  Crocco  and  Probstein60  was  uced  for  the  flow  normal  to  the  step 
face  and  appeared  to  correlate  Stalker’s  results.  The  Mach  number,  normal  to  the  step  face,  was  the 
flow  parameter  that  most  affected  the  plateau  pressure  rise,  as  it  did  the  extent  of  upstream  influence. 
Miller  and  Redeker66  also  reported  that  a simple  modification  of  Chapman’s67  ord  .i  cl  magnitude 


* 


‘free-interaction’  implies  that  the  region  of  disturbed  flow  up  to  the  separation  point  (in  two- 
dimensional  flow)  is  free  from  the  influence  of  flow  downstream. 
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analysis  provided  an  effective  means  to  account  for  sweep  effects  on  plateau  and  separation  pressures. 
One  evident  difficulty  in  likening  the  swept  interaction  to  the  two-dimensional  interaction  is  that  the 
boundary-layer  integrals  of  flow  quantities  resolved  in  a direction  perpendicular  to  the  swept  shock  are 
not  typical  of  their  edge  Mach  number.  This  is  because  density  (like  temperature)  is  a scalar,  and  its 
distribution  across  the  boundary  layer  characterizes  the  flow  in  the  streamwise  direction.  We  see, 
therefore,  that  even  though  the  behaviour  of  the  flow  perpendicular  to  the  swept  shock  may  appear 
quasi/two-dimensional,  the  various  integral  thicknesses  calculated  in  a direction  normal  to  the  shock, 
will  be  distorted  by  an  uncharacteristic  density  distribution. 

Stanbrook53  made  some  exploratory  investigations  of  a swept  shock/turbulent  boundary-layer 
interaction  on  the  sidewall  of  a 4 X 4- in.  wind  tunnel  at  Mach  numbers  1.6, 1.8  and  2.0,  at  Reynolds 
numbers  based  on  the  0.14-in.  boundary-layer  thickness  of  0.25  X 10s . The  shock  was  generated  by  a 
wedge  standing  perpendicular  to  the  test  wall.  It  was  found  that  incipient  separation  occurred  at  a 
pressure  ratio  across  the  shock  wave  roughly  equal  to  1.5,  which  corresponded  to  a flow  deflection  of 
between  7.5  and  8-degrees.  This  ratio  was  substantially  lower  than  that  found  in  2D  flow,  where  the 
value  was  between  1.6  and  1.8  in  the  same  Mach  number  range.  At  distances  of  more  than  ten  boundary- 
layer  thicknesses  from  the  shock  generator,  Stanbrook  expected  the  flow  to  approach  asymptotically 
a cylindrical  (that  is,  infinite  swept  back)  form. 

McCabe54  used  a wind  tunnel  configuration  not  unlike  Stanbrook’s,  to  conduct  experiments 
at  Mach  numbers  of  1.96  and  2.94  at  Reynolds  numbers  close  to  0.4  X 10s  based  on  the  boundary-layer 
thickness  of  0.14-in.  To  determine  the  onset  of  3D  separation,  Maskell’s  criterion  of  the  formation  of 
an  envelope  of  the  limiting  streamlines  was  used,  and  above  a pressure  ratio  of  1.5,  McCabe  also  found 
3D  separation  to  exist  at  both  test  Mach  numbers.  Like  Stanbrook,  McCabe  indicated  that  at  a distance 
of  about  20  boundary-layer  thicknesses  from  the  wedge  apex,  measured  along  the  shock,  the  flow  pat- 
tern tended  towards  a cylindrically  symmetric  form.  He  also  presented  an  approximate  theory  to  relate 
the  deflection  of  the  flow  close  to  the  surface  to  the  changes  in  the  mainstream.  Upstream  of  the  inter- 
action, an  element  of  vortex  tube  was  defined  whose  height  was  equal  to  the  boundary-layer  thickness. 
The  circulation  around  the  element  was  determined  and  ssumed  to  be  conserved  as  it  was  transported 
downstream  through  the  interaction  at  the  mainstream  velocity.  The  flow  was  assumed  inviscid  and 
slip  allowed  at  the  wJl.  McCabe  sliced  the  element  downstream  of  the  shock  in  a place  perpendicular 
to  the  local  mainstream  direction,  to  determine  a ‘cross  flow’  slip  velocity.  He  added  this  vectorially  to 
the  streamwise  velocity  in  the  ‘boundary  layer’,  that  was  taken  equal  to  the  local  mainstream  velocity, 
to  obtain  the  flow  direction  at  the  surface.  This  latter  assumption  of  uniform  streamwise  velocity  in 
the  downstream  boundary  layer  is  a gross  simplification  of  the  real  flow,  although  McCabe’s  measured 
directions  of  the  wall  limiting  streamlines  agreed  satisfactorily  with  his  theory  for  flow  deflections  up 
to  those  when  separation  was  approached  (see  Fig.  10). 

Lowrie56  performed  some  experiments  in  a 3 X 8-in.  working  section  of  a blowdown  wind 
tunnel,  the  8-in.  width  being  used  for  his  test  surface.  A shock  generator  perpendicular  to  the  test  wall 
provided  the  swept  interaction.  Tests  were  conducted  at  Mach  numbers  of  2.50  and  3.44  at  nominal 
Reynolds  numbers  of  0.3  X 106  based  on  an  undisturbed  boundary -layer  thickness  of  0.3-in.  Lowrie’s 
measurements  included  wall  static  pressures,  surface  oil  flow  patterns  and  pitot  and  (chamfered  two 
tube)  yaw  probe  traverses  through  the  test  boundary  layer.  Incipient  separation  was  quoted  as  occur- 
ring at  flow  deflection  angles  of  about  IV2  and  5-degrees  at  Mach  numbers  of  2.50  and  3.44.  These 
conditions  corresponded  with  overall  static  pressure  ratios  across  the  shock  wave  of  about  1.60  and 
1.50  respectively.  Lowrie  provided  some  analysis  of  the  development  of  cross  flows  within  the  inter- 
action, by  extending  the  compressible  flow,  Johnston’s61"63  suggestion  that  over  most  of  the  boundary 
layer,  the  cross  flow  velocity  is  proportional  to  the  streamwise  velocity  defect. 

Johnston’s  polar  plot  of  velocity  is  shown  on  Figure  11(a).  The  velocity  u is  parallel  to  the 
mainstream  at  the  boundary-layer  edge,  and  the  cross  flow  velocity,  v,  is  in  a transverse  direction  within 
the  boundary  layer  and  normal  to  the  local  mainstream.  The  line,  PE,  of  the  outer  region  represents 
the  virtually  ‘inviscid’  response  to  transverse  pressure  gradients  of  the  flow  between  the  viscous  wall 
region  and  the  mainstream.  The  inner  line,  OP,  is  more-or-less  straight  signifying  co-planar  velocity 
profiles  in  the  wall  region  that  is  dominated  by  viscous  effects.  The  junction  point  where  the  inner  and 
outer  regions  join  is  not  found  to  be  discontinuous  in  experiments  (see  Ref.  56  for  example). 
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Provided  the  pressure  gradient  was  the  dominant  feature  of  the  flow,  and  hence,  the  cross 
flow  angle  was  increasing,  Lowrie  derived  a criterion  for  the  junction  by  considering  the  equilibrium 
of  the  inner  layer  on  the  basis  of  a balance  between  shear  and  pressure  forces  (assuming  inertial  effects 
were  negligible).  If  in  the  outer  region  the  shear  stress  is  in  the  direction  of  the  shear  strain  (along  PE), 
an  elemental  balance  between  pressure  and  shear  forces  perpendicular  to  this  line  permits  the  height  of 
the  inner  region  to  be  found  as  that  over  which  the  pressure  force  is  just  equal  to  the  normal  component 
of  only  the  wall  shear  stress.  If  the  law  of  the  wall  is  applicable  to  the  inner  region,  the  wall  shear  stress 
direction  may  be  determined  in  relation  to  the  slope  of  PE,  the  local  pressure  gradient  perpendicular  to 
PE,  and  the  wall  friction  velocity,  urw . Figure  11(b)  demonstrates  satisfactory  agreement  between 
Lowrie’s  predicted  magnitude  of  the  velocity,  qp,  at  the  junction  point,  when  tested  against  some  low 
speed  data. 

Green30  quotes  that  the  direction  of  the  surface  friction  can  be  adequately  estimated  by 
letting  the  vector,  qp , be  proportional  to  ur  w (based  on  local  wall  conditions)  and  the  ratio  of  the 
local  to  upstream,  undisturbed  mainstream  velocity,  uE1 , in  the  form: 

Qp  w 

— = 19.5  

UE  UE 1 

Using  the  sine  rule  on  Figure  11(a),  the  angle,  /3,  between  the  surface  shear  stress  direction 
and  the  local  flow  direction  u.  the  boundary-layer  edge  is 


j8  = arcsin 


1 

(Qp/ue) 


sin  <p 


- <t> 


where  Green30  suggests  that  tan  <j>  is  the  ratio  of  the  transverse  volume  flow  integral  or  cross  flow  dis- 
placement thickness  to  the  stream  wise  displacement  thickness  integral  (=  62  */5 , *). 

The  polar  plot  profiles  of  cross  flow  velocity  versus  streamwise  velocity  that  Lowrie  illustrates 
in  his  thesis56 , provide  some  evidence  of  triangularity.  Nevertheless,  Lowrie  omitted  altogether  any 
boundary -layer  pitot  profiles  where  there  were  indications  of  wave  patterns  in  the  outer  region,  and 
so  perhaps  lost  some  of  the  interaction  details  near  the  position  of  the  calculated  shock  wave  location*. 
In  any  event  it  would  seem  doubtful  that  any  profile  form  based  on  equilibrium  boundary-layer  results 
would  apply  in  the  non-equilibrium  interaction  region. 

Finally,  in  this  review  of  three-dimensional  interaction  phenomena  we  will  mention  the  super- 
sonic interference  flow  along  the  corner  of  intersecting  wedges  discussed  at  length  in  References  68-7  4, 
and  the  side-mounted  supersonic  intake  problem,  discussed,  for  example,  in  References  75  and  76. 

The  work  on  the  symmetrical  streamwise  corner  problem68'74,  where  two  wedges  of  equal 
deflection  to  the  oncoming  stream  are  jointed,  has  identified  the  complex  shock  structure,  and  the 
inviscid  flow  features.  There  are  strong  cross-stream  pressure  and  entropy  gradients  both  at  the  comer 
and  away  from  it,  which  presumably  should  not  be  disregarded  in  analyses  of  the  boundary  layers. 

The  importance  of  comer  interactions  in  high  supersonic  and  hypersonic  flows  is  that  they  can  give 
rise  to  high  local  heating  rates  at  reattachment  zones.  The  wedge  bow  shock  waves  do  not  intersect, 
but  are  joined  by  a third  comer  shock.  The  internal  flow  includes  two  embedded  shocks  which  ter- 
minate at  the  wedge  surfaces,  and  a triangular  region  bounded  by  two  slip  surfaces  and  the  comer 
shock.  The  embedded  shocks  interact  with  the  boundary  layers  on  the  wedge  surfaces  and,  depending 


* The  calculated  location  of  the  oblique  shock  wave  in  an  inviscid  supersonic  mainstream  (upstream 
and  downstream  of  the  shock  wave)  is  sometimes  referred  to  as  ‘inviscid  shock  wave  location’. 

The  nomenclature  is  admittedly  odd,  but  it  does  provide  a shorthand  notation. 
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on  the  wedge  angles,  may  be  strong  enough  to  cause  separation.  Kutler7  3 has  modelled  the  conical  flow 
generated  by  two  intersecting  wedges  immersed  in  a supersonic  stream,  and  obtained  a numerical  solu- 
tion for  the  inviscid  comer  flow  surface  pressure  distributions,  which  compare  quite  favourably  with 
the  experiments  of  West  and  Korkegi7 1 . 


In  flow  field  measurements  in  an  asymmetric  comer  at  Mach  12.5,  where  laminar  boundary 
layers  were  present,  and  where  the  wedge  angle  of  incidence  was  15-degrees,  Cooper  and  Hankey  re- 
vealed two  large  vortices  within  the  floor  boundary  layer  that  were  responsible  for  high  local  heating. 
These  appeared  to  emanate  from  three-dimensional  separation  lines.  A complex  shock  pattern  domi- 
nated by  a triple  point  structure  was  determined,  essentially  conical  in  nature.  Cooper  and  Hankey’s 
asymmetric  corner  corresponds  with  the  wedge  mounted  normal  to  the  flat  test  wall  that  we  have  pre- 
viously discussed,  except  now  the  test  boundary  layer  commences  at  the  leading-edge  of  the  wall,  which 
is  the  same  station  as  the  leading-edge  of  the  shock  generator. 


A side-mounted  supersonic  half-cone  intake  centre  body,  causes  junction  problems76  not 
unlike  those  already  discussed.  It  was  demonstrated  in  Reference  75  that  the  quasi-conical  half-cone 
pressure  field  impressed  on  the  adjacent  surrounding  wall  of  a fuselage,  could  be  modelled  experimen- 
tally by  choosing  a suitable  deflection  angle  of  a wedge  mounted  normal  to  a flat  plate.  Pressure 
distributions  between  the  half -cone  and  its  shock  wave  (calculation)  and  the  wedge  and  its  shock  wave 
(experiment)  were  compared  as  shown  in  Figure  9 of  Reference  75,  with  relatively  good  agreement. 
Downstream  of  the  wedge  apex,  the  interacting  viscous  and  non-viscous  flows  produced  a quasi-conical 
flow. 

1.7  Prospectus 

We  have  noted  then,  some  agreement  and  some  contradictory  findings  in  the  three-dimensional 
swept  shock/turbulent  boundary-layer  work  already  reported.  We  shall  endeavour  in  our  experiments 
to  present  a fuller  picture  of  the  interaction  region,  by  measuring  at  close-spaced  intervals,  boundary- 
layer  pitot,  total  temperature  and  yaw  angle  profiles,  wail  temperatures  and  pressures,  and  surface  shear 
stress  directions.  Such  measurements  will  be  conducted  with  a swept  shock  wave  standing  normal  to 
the  test  surface  of  strength  near  to  that  producing  incipient  separation  at  Mach  2,  to  provide  a back- 
ground for  interpreting  the  later  results  of  blowing  boundary-layer  control;  and  at  shock  pressure  ratios 
slightly  more  and  well  above  those  required  for  incipient  separation  at  Mach  4.  Where  appropriate,  the 
present  results  will  be  tested  against  those  obtained  by  previous  experimenters. 


2.0  APPARATUS 

2.1  The  5 X 5-IN.  BLOWDOWN  WIND  TUNNEL 

The  5 X 5-in.  blowdown  wind  tunnel  was  constructed  as  a 1-in.  to  1-ft.  scale  aerodynamic 
model  to  evolve  the  aerodynamic  design  of  the  NAE  5 X 5-ft.  trisonic  blowdown  wind  tunnel.  The 
5-in.  tunnel  has  interchangeable  transonic  and  supersonic  test  sections,  but  only  the  supersonic  nozzle 
block  assembly  was  used  for  the  tests  described  herein. 

The  5-in.  tunnel  may  be  supplied  with  compressed  air  at  300  psia  from  the  50,000  ft.3  storage 
tanks  connected  to  the  5-ft.  facility;  or  alternatively,  it  can  use  air  from  smaller  reservoirs  fed  by  15-hp 
and  24-hp  compressor  sets  in  series.  Using  the  small  reservoirs,  the  run  time  is  restricted  to  about  20 
seconds,  but  the  main  tunnel  storage  tanks  provide  the  means  for  continuous  running.  The  absolute 
humidity  of  the  air  in  the  large  reservoirs  is  kept  below  0.0002  lb.  water  per  lb.  mass  of  air.  The  settling 
chamber  air  is  stabilized  via  thermal  matrices  and  is  typically  about  70°  F,  with  no  provision  for  heat 
addition,  so  tha*  the  limiting  Mach  number  to  which  the  flow  may  be  expanded  is  about  4.2  -►  4.5 
without  condensation  and  liquefaction  problems  occurring.  Upon  initiation  of  a tunnel  run,  the  air 
passes  from  either  of  the  reservoir  facilities  along  pipes  that  connect  to  the  automatic  pressure  control 
valve.  This  sleeve  valve  maintains  the  settling  chamber  pressure  at  a constant,  pre-selected  value.  From 
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the  settling  chamber,  the  air  exhausts  to  atmosphere  via  the  nozzle,  working  section  and  variable  geom- 
etry supersonic/subsonic  diffuser.  Figure  12  indicates  tht  scheme  of  the  5-in.  tunnel  installation. 

The  supersonic  nozzle  section  of  the  tunnel  is  compoased  of  flat  top  and  bottom  plates  (the 
‘sidewalls’)  that  are  separated  by  two  contoured  nozzle  blocks,  see  Figure  13.  The  supersonic  expansion 
of  the  flow  ends  about  6.49-in.  upstream  of  the  exit  of  the  nozzle  blocks;  the  length  of  working  section 
is  this  latter  dimension,  along  which  the  nozzle  blocks  possess  a straight,  divergent  taper  to  correct  for 
the  boundary-layer  displacement  thickness  and  maintain  a constant  Mach  number  in  the  test  section. 

One  of  the  nozzle  blocks  has  two  slots  to  accept  the  support  legs  of  the  wedge  shock  generator  (Fig.  13). 

In  the  opposite  nozzle  block,  there  is  a circular  window  to  provide  artificial  illumination  of  the  tunnel 
interior.  There  are  nozzle  blocks  available  for  mainstream  Mach  numbers  of  1.5,  2.0,  3.0  and  4.0. 

The  test  boundary-layer  flow  is  along  the  sidewall  that  is  also  the  ‘floor’  of  the  tunnel.  In 
the  region  of  the  working-section,  there  is  a rectangular  cut-out  in  the  test  floor  that  is  filled  with  an 
instrumented,  removable  insert  that  provides  the  necessary  instrumentation  and  equipment  (the  blowing 
nozzle  for  boundary-layer  control,  for  example)  for  any  particular  test  requirement.  Three  such  inserts 
are  illustrated  on  Figures  15,  22  and  24,  and  are  the  subjects  of  Paragraphs  2.3,  2.5  and  2.6. 

Tunnel  total  pressures  may  be  taken  up  to  at  least  12  atmospheres,  giving  a Reynolds  number 
of  approximately  7 X 106  based  on  the  tunnel  cross-section  dimension  of  5-in. 

An  additional  supply  pipe  for  compressed  air  to  the  blowing  nozzle  for  boundary-layer  control 
is  affixed  to  the  floor  beneath  the  wind  tunnel  working-section.  Stagnation  pressures  up  to  300  psia  are 
also  available  for  this  blc  supply. 

There  is  a large  glass  window  in  the  tunnel  top  plate.  Figure  13,  that  may  be  quickly  unclamped 
to  apply  flow  visualization  material  to  the  test  floor,  or  to  view  and  photograph  the  flow  visualization 
patterns.  A Schlieren  system  may  be  set  up  if  transparent  floor  panels  are  used.  The  settling  chamber 
stagnation  pressure  and  temperature  were  measured  respectively  on  50  psia  or  250  psia  pressure  trans- 
ducers, and  with  a glass  bulb  thermometer  and  calibrated  diode. 

2.2  The  Movable  Shock  Generator 

The  shock  (or  expansion)  generator  is  a very  thin  triangular  profile  supersonic  aerofoil  mounted 
vertically  between  the  tunnel  sidewalls,  whose  leading-edge  and  trailing-edge  included  angle  are  4.8-degrees. 
The  foot  of  the  generator  makes  firm  contact  with  the  test  floor  of  the  working  section  and  a brass  edge 
ensures  non-abrasion  of  the  steel  floor.  The  installation  of  the  generator,  hereinafter  called  ‘wedge’, 
and  its  dimensions  are  shown  on  Figures  13  and  16  respectively.  Slots  in  the  tunnel  wall  allow  passage 
for  the  two  support  struts  of  the  wedge,  which  are  attached  to  mechanisms  outside  of  the  wind  tunnel 
that  provide  wedge  movement  in  translation  and  iotation. 

The  wedge  may  be  rotated  about  its  leading-edge  between  the  angular  limits  of  20-degrees 
compression  and  15-degrees  expansion,  when  referenced  to  the  tunnel  axis.  It  may  also  translate 
parallel  to  the  tunnel  centre-line  and  0.985-in.  from  it  for  a total  excursion  of  6-in.  where  the  most 
downstream  position  of  the  leading-edge  that  may  be  selected  is  5.25-in.  ahead  of  the  exit  of  the 
nozzle  blocks.  The  useful  movement  of  the  wedge  is  indicated  on  Figures  15,  22,  and  24. 

The  cross-sectional  thickness  of  the  wedge  support  legs  and  their  spacing  were  chosen  to 
avoid  choking  in  the  small  channel  between  the  tunnel  wall  and  the  wedge  for  the  range  of  expansion 
angle  mentioned  above.  The  maximum  stagnation  blowing  pressures  that  were  used  in  the  tests  were 
chosen  so  that  with  a normal  shock  spanning  the  test  flow  between  the  sidewall  and  the  flat  face  of  the 
wedge,  the  design  stress  levels  in  the  wedge  supports  and  attachments  were  not  exceeded  (that  is,  the 
wedge  was  stressed  to  handle  the  so-called  starting  loads).  These  maximum  pressures,  for  the  four 
available  tunnel  mainstream  Mach  numbers  are  listed  on  Table  1.  The  region  within  the  hatching  on 
the  tunnel  performance  chart  (Fig.  17)  shows  the  range  of  flow  conditions  available  for  safe  operation 
of  the  wedge. 
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Hydraulic  jacks  provide  the  means  to  move  the  wedge  in  its  axial  and  rotation  modes,  as  we 
see  on  Figure  14.  Usually,  the  axial  location  of  the  wedge  was  set  prior  to  the  tunnel  run,  and  for  wedge 
deflection  angles  up  to  8-degrees,  the  wedge  angle  was  also  preset.  For  angles  much  larger  than  8-degrees, 
however,  the  internal  contraction  ratio  was  too  large  to  establish  supersonic  flow  throughout  the  working 
section  and  additional  movement  in  rotoation  was  only  attempted  once  the  tunnel  run  was  in  progress. 
Both  of  the  piston  strokes  were  calibrated  to  provide  direct  readings  of  the  wedge  angular  and  transla- 
tional positions,  and  were  checked  frequently  during  the  course  of  the  tests.  The  respective  movements 
could  be  guaranteed  to  better  than  ±0.1-degree  in  rotation  and  ±0.010-in.  in  distance  setting.  The  de- 
flection angle  of  the  wedge  was  also  checked  by  measuring  static  pressures  at  several  mid-chord  stations 
on  the  wedge  face  adjacent  to  the  tunnel  axis,  and  comparing  these  in  terms  of  the  local  pressure  ratio 
and  mainstream  Mach  number,  with  values  of  flow  deflections  listed  in  oblique  shock  tables.  (This 
latter  check  was  unreliable  when  a significant  region  of  3D  separation  existed.)  The  agreement  between 
the  wedge  angle  and  that  assessed  from  the  wedge  pressures  was  usually  within  the  movement  tolerance, 
up  to  a wedge  deflection  angle  of  8-degrees. 


2.3  Plain  Floor 


Figure  15  shows  a plan  view  of  the  dummy  floor  provided  to  measure  the  wall  static  pressures 
beneath  the  boundary  layer  in  the  working-section  without  the  shock  generator;  and  through  the  swept 
shock/ boundary-layer  interaction  with  the  wedge  in  situ.  There  are  more  than  200  static  pressure  orifices 
contained  in  the  insert,  and  in  the  fixed  part  of  the  tunnel  wall  further  downstream,  along  lines  parallel 
and  perpendicular  to  the  tunnel  centre-line.  The  spanwise  and  longitudinal  wall  static  distributions  may 
hence  be  obtained  virtually  anywhere  in  the  working  section,  while  provision  was  made  at  two  transverse 
stations  for  three  pitot  traverses  at  each  to  investigate  the  profiles  in  the  nominal  two-dimensional 
boundary  layers.  All  of  the  static  pressure  orifices  are  0.015-in.  in  diameter  excepting  those  through 
which  the  pitot  tube  stems  passed,  which  were  enlarged  to  0.032-in.  diameter  (see  Rows  4 and  18  on 
Fig.  15).  The  locations  of  the  static  pressure  orifices  were  referenced  to  a station  8.625-in.  upstream 
of  the  geometric  exit  of  the  tunnel  nozzle.  This  was  labelled  as  x = 0.0  and  v'as  the  centre-line  of  the 
turntable  blowing  assembly  to  be  described  in  Section  2.6. 


The  wall  static  pressure  points  are  connected  to  automatic  scan  wafer  valve  assemblies  and 
sampled  at  one  point  per  second.  The  pressures  are  measured  with  differential  or  absolute  ‘Statham’ 
unbonded  strain  gauge  type  transducers  and  thence  to  a digital  voltmeter  or  an  analogue  chart  recorder. 
All  static  pressures  were  referenced  to  that  measured  at  Station  1,  that  was  situated  upstream  in  the 
test  section  just  within  the  test  diamond  at  Mach  2,  as  we  see  on  Figure  15. 

2.4  Traverse  Gear  to  Investigate  the  Nominal  Two-Dimensional 
Boundary  Layer  in  the  Empty  Wind  Tunnel 


The  equipment  shown  on  Figure  18  permits  the  remote  positioning  by  an  electric  servo  of 
one  to  three  circular  pitot  probes  to  traverse  the  nominal  two-dimensional  boundary  layer  along  the 
tunnel  floor.  A traverse  height  of  0.5-in.  is  available  in  any  given  run.  The  probe  is  inserted  into  one 
of  the  enlarged  static  pressure  orifices  from  inside  the  working  section  and  clamped  to  the  probe  holder 
outside  of  the  wind  tunnel.  The  external  diameter  of  the  pitot  tube  is  0.0164-in.  The  pitot  pressure 
is  fed  to  either  a differential  or  an  absolute  pressure  transducer,  depending  on  the  blowdown  pressure 
and  Mach  number  of  the  run.  At  the  wall,  the  hook  shape  of  the  probe  enables  it  to  be  used  as  a 
Preston  tube  to  measure  local  skin  friction.  Figure  18(c)  shows  the  dimensions  of  a Preston  tube  device 
whose  stem  may  be  inserted  into  one  of  the  enlarged  static  pressure  orifices  along  Lines  4 and  18.  The 
hook  geometry  of  each  traversing  circular  pitot  probe  was  identical  with  that  depicted  on  Figure  18(c). 


The  servo  has  two  sections,  the  probe  unit  mounted  to  the  wall  that  we  see  on  Figure  18;  and 
the  control  chassis  to  which  the  probe  unit  is  connected  with  shielded  cables. 

The  probe  unit  consists  essentially  of  a two  rpm  reversible  dc  motor  with  an  integral  reduc- 
tion gearbox,  an  error  correcting  element  in  the  form  of  a ten-turn  potentiometer;  and  a worm  gear 
driving  a slider,  to  which  the  probes  are  attached.  The  probe  unit  may  be  located  at  various  lengthwise 
stations  along  the  tunnel  working-section  by  slackening  the  wing  nuts  shown  on  Figure  18(b)  and 
sliding  the  probe  unit  along  its  mounting  rails. 


The  control  servo  determines  the  direction  of  motor  shaft  rotation  and  the  position  of  the 
probe  relative  to  the  wall.  The  control  chassis  has  three  parts:  a position  control  and  two  power  chan- 
nels, each  channel  consisting  of  a low  frequency  amplifier  and  an  electronic  switch.  A voltmeter  indicates 
a decreasing  signal  amplitude  as  the  probe  approaches  its  dialled  position  set  on  the  control  panel.  Limit 
switches  are  incorporated  to  protect  the  error-correcting  potentiometer  from  damage  that  would  result 
from  motor  overshoot.  Point-by -point  measurements  are  taken  with  this  system,  rather  than  a continuous 
traverse. 

2.5  ‘Cobra’  Probe  Traverse  Gear  to  Investigate  the  Three-Dimensional  Shock 
Wave/Boundary-Layer  Interaction  Region 

The  pneumatic  sensor  part  of  this  tiny  boundary-layer  probe  comprises  a flattened  three-tube 
measuring  head,  where  the  two  outer  tubes  are  each  chamfered  back  at  45-degrees  to  form  a yawmeter, 
while  the  centre  tube  (0.028-in.  wide  X 0.012-in.  high)  is  a rectangular  pitot.  Figure  19  illustrates  the 
neat  and  tidy  design  appearance  of  the  probe  head,  the  shape  of  which  bears  some  resemblance  to  the 
reared  head  of  a cobra  snake  — hence  the  name,  ‘cobra’  probe.  Figure  20,  on  the  other  hand,  shows  the 
less  tidy  frontal  view  of  the  manufactured  probe  head  when  viwed  through  a microscope  at  approximately 
50  times  true  size. 

Each  time  a cleaning  tool  is  passed  across  the  brass  faces  of  the  pneumatic  tubes,  it  is  advisable 
to  re-check  the  probe  dimensions,  for  small  changes  occur. 

As  the  probe  is  traversed  through  the  viscous  flow,  from  an  initial  position  0.0005-in.  from 
the  wall  to  a maximum  extension  of  about  0.8-in.,  a yaw  servo-system  turns  the  probe  head  to  align  it 
with  the  oncoming  flow  direction.  Local  flow  angles  in  a plane  parallel  to  the  test  floor,  and  pitot 
pressures,  can  thus  be  measured  in  a skewed  boundary  layer.  Probe  height  above  the  surface  and  yaw 
angle  are  measured  by  geared  multi-turn  potentiometers.  The  pitot  pressure  was  recorded  by  a close- 
coupled  50  psia  pressure  transducer.  The  probe  extension  rate  can  be  preselected:  the  fastest  rate  was 
chosen  that  was  consistent  with  obtaining  the  same  pitot  profile  at  the  slowest  available  traverse  rate, 
and  a traverse  height  of  0.5-in.  was  usually  sampled  in  about  10-15  seconds. 

A 0.5-in.  diameter  flattened  steel  cylinder  into  which  the  triangular-section  probe  stem  is 
effectively  keyed  provides  the  means  of  yawing  the  probe.  Figure  21  shows  the  assembly  of  the  prc  be 
and  its  appearance  above  a test  wall. 

To  measure  the  variation  of  stagnation  temperature  through  the  viscous  flow,  an  unshielded 
iron-constantan  thermocouple,  of  bead  diameter  0.002-in.,  was  mounted  above  the  pitot  orifice  (Figs  19 
and  20).  The  distinct  advantage  of  this  installation  was  its  small  size  and  its  relative  robustness  in  a 
sometimes  contaminated  wind  tunnel  stream. 

As  the  probe  was  required  for  use  in  relatively  thin  boundary  layers  neither  the  shielded 
probes  (Winkler77 , Meier78 , East79 ) nor  the  hot-wire  type  of  Vas80  were  considered  appropriate. 

Instead,  the  design  of  stagnation  line  probe  introduced  by  Bradfield  et  al81  and  Harris82 , was  adapted 
to  the  cobra  probe  head  thereby  enabling  pneumatic  and  temperature  measurements  to  be  taken  simul- 
taneously. The  thermocouple  wires  were  ‘sandwiched’  between  two  pieces  of  cigarette  paper  and  the 
assembly  soaked,  secured  and  insulated  from  the  probe  head  by  epoxy  resin.  The  thickness  of  the 
sandwich  added  0.004-in.  to  the  probe  height.  The  hair- thin  wires  were  then  affixed  in  slots  at  the 
back  of  the  probe  support  stem  and  connected  in  series  with  a reference  temperature  source  at 
103-degrees  F. 

The  cobra  probe  assembly  was  mounted  to  a special  floor  insert,  at  one  of  four  available 
stations,  see  Figure  22.  The  Stations  T3  and  T4  were  located  along  a line  at  8-degrees  to  the  tunnel 
centre-line,  which  would  be  in  the  direction  of  the  calculated  inviscid  flow  streamline  for  the  same 
wedge  compression  angle  with  the  wedge  leading-edge  situated  at  XWLE  = 0.0-in. 

The  same  floor  insert  was  also  used  to  accommodate  embedded  thermocuples  to  measure  the 
local  near  adiabatic  wall  temperature.  Clearly,  the  principal  requirement  is  that  the  thermocouple  device 
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should  not  change  the  conditions  of  heat  transfer  that  would  exist  if  it  were  absent.  Hence,  the  wall 
temperature  measurement  was  made  by  setting  a thin  (0.010-in.)  steel  plate  into  a 0.5-in.  outside  diam- 
eter insulated  cylinder,  which  was  in  turn  inserted  into  Stations  T1  through  T4,  so  that  the  plate  surface 
was  flush  with  the  surrounding  wall.  To  the  backside  of  the  thin  plate  element  the  same  size  of  0.002-in. 
diameter  thermocouple  bead  was  secured  as  used  on  the  cobra  head.  The  wall  temperature  measuring 
device  is  illustrated  in  section  on  Figure  23.  Such  thin  wall  plates  attain  equilibrium  conditions  quickly83, 
and  data  were  taken  throughout  a tunnel  run  to  verify  that  a plateau  in  temperature  was  indeed  achieved. 
The  wall  temperature  inserts  were  coupled  to  the  same  electronic  circuitry  as  the  stagnation  probe,  a 
switching  device  permitting  the  selection  of  any  one  of  six  temperature  channels. 

All  of  the  cobra  probe  outputs  were  amplified  and  registered  on  X-Y  charts,  while  the  cobra 
and  wall  temperature  voltages  (again  amplified)  were  also  displayed  on  a Hewlett-Packard  digital  volt- 
meter. 

2.6  Floor  Insert  Containing  the  Supersonic  Blowing  Nozzle  for  Boundary-Layer  Control 

The  arrangement  of  the  blowing  slot  for  tangential  air  injection  at  a chosen  supersonic  Mach 
number  is  shown  on  Figures  13,  24  and  25.  The  two-dimensional  blowing  nozzle  profile  is  illustrated 
on  Figure  26.  The  straight  nozzle  surfaces  in  the  supersonic  expansion  region  are  at  4-degrees  to  the 
line  of  the  test  floor,  these  being  prescribed  to  overcome  the  problem  of  machining  an  accurately  pro- 
filed upper  wall  with  an  associated  very  thin  lip  (0.005-in.  radius),  yi  Iding  the  possibility  of  lip  deflec- 
tions. The  straight  chamfer  to  the  thin  lip  provides  rigidity  to  the  iip,  as  well  as  permitting  the  use  of  a 
small  tip  radius  to  minimize  wake  disturbances  downstream. 

Provision  is  included  in  the  design  of  the  nozzle  to  rotate  the  injection  slot  exit  at  different 
angular  settings  relative  to  the  shock  wave  angle.  The  blowing  nozzle  is  part  of  a turntable,  the  centre 
of  which  is  located  on  the  tunnel  centre-line,  8.625-in.  ahead  of  the  end  of  the  tunnel  nozzle  blocks, 
see  Figure  24.  All  of  the  locations  of  the  static  pressure  orifices  in  the  plain  floor  and  in  the  floor  con- 
taining the  blowing  nozzle  are  referenced  to  the  turntable  centre  of  rotation.  The  range  of  angular 
movement  of  the  turntable  is  ±60-degrees  about  the  transverse  normal  to  the  mainstream  flow  direction, 
a protractor  scale  on  the  underside  of  the  tunnel  indicating  the  angle. 

The  exit  height  of  the  blowing  nozzle  is  0.1-in.  but  may  be  doubled  using  a spacer  beneath 
the  upstream  nozzle  block.  Because  of  this  finite  slot  height  above  the  nominal  level  of  the  test  wall, 
a 3-degree  ramp  leads  the  undisturbed  boundary  layer  over  the  top  of  the  injection  nozzle,  as  we  see 
on  Figure  25.  The  span  of  the  slot  is  3.6-in. 

To  change  the  Mach  number  of  the  injectant  air,  a series  of  blocks  has  been  made  to  alter  the 
throat  dimension,  see  Figure  26.  The  nominal  injection  Mach  numbers  range  from  1.5  to  4.5,  the  ‘A’ 
and  ‘B’  dimensions  indicated  or  Figure  23  and  the  design  pressures  of  the  injected  flow  being  annotated 
on  Table  2.  The  contraction  ratio  from  the  entrance  of  the  blowing  nozzle  to  the  throat  varies  from 
10:1  at  Mj  = 1.5  to  145:1  for  M,  = 4.5.  The  throat  radius  of  curvature  changes  from  1.76  to  24.8  times 
the  throat  height  for  the  same  two  Mach  numbers,  both  values  exceeding  those  necessary  for  satisfactory 
operation84  of  an  effuser.  A transition  duct,  of  virtually  constant  area,  takes  the  blowing  air  from  the 
circular  supply  pipe  to  the  rectangular  entrance  of  the  blowing  nozzle  settling  chamber. 

The  layout  of  the  0.015-in.  diameter  static  pressure  orifices  about  and  downstream  of  the 
blowing  slot  are  shown  on  Figure  24.  There  are  155  static  pressure  points  in  the  floor  and  at  the  throat 
of  the  blowing  nozzle.  In  addition,  the  state  of  the  blowing  air  supply  is  measured  via  a pitot  and  static 
pressure  orifices,  and  a glass  bulb  thermometer  in  the  transition  duct.  As  before,  the  static  pressures 
were  collected  via  waft  r valves  to  pressure  transducers,  and  read  subsequently  from  analogue  charts. 


3.0  TEST  CONDITIONS 

3.1  Stagnation  Pressures  in  the  5 X 5-in.  Tunnel 


1 
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The  maximum  value  of  tunnel  stagnation  pressure  is  governed  by  the  strength  of  the  wedge 
shock  generator  supports.  Maximum  safe  values  of  blowing  pressure  versus  mainstream  Mach  number 
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are  shown  on  Table  1,  with  the  corresponding  nominal  starting  pressures  and  Reynolds  numbers  based 
on  the  tunnel  span  of  5-in.  for  the  four  available  supersonic  nozzle  blocks.  Only  the  Mach  2 and  Mach  4 
nozzles  were  used  in  the  current  investigation.  The  available  run  conditions  are  illustrated  graphically 
on  Figure  17  within  the  hatched  border,  the  usable  stagnation  pressures  increasing  from  two  to  twelve 
atmospheres  for  the  two  chosen  test  Mach  numbers  (q  is  dynamic  pressure). 

3.2  Blowing  Quantities  for  Boundary-Layer  Control 

The  jet  stagnation  pressures  listed  on  Table  3 indicate  those  required  to  give  the  static  pressure 
at  the  jet  nozzle  exit  equal  to  the  local  ambient  pressure  of  the  supersonic  mainstream  for  1.5  < Mj  < 4.5. 
The  jet  excess  momentum  values  for  each  blowing  condition  taking  the  span  of  the  blowing  slot  equal 
to  3.6-in.  (Table  3)  are  plotted  on  Figure  27,  with  reference  to  the  calculated  momentum  deficit  of  the 
undisturbed  tunnel  wall  boundary  layer  at  the  jet  slot  position  (Table  4)  for  the  listed  tunnel  stagnation 
pressures.  We  note  that  at  a mainstream  Mach  number  of  2,  a jet  nozzle  Mach  number  of  3 provides  a 
jet  excess  momentum  somewhat  more  than  the  ‘old’  boundary-layer  momentum  deficit.  This  injection 
Mach  number  of  3 was  used  throughout  the  M,  = 2 boundary-layer  control  experiment  pursued  herein. 

4.0  RESULTS 

The  essential  results  of  the  two  experiments  are  shown  respectively  on  Figures  43  to  87  for 
the  swept  interaction  phenomena  at  Mach  numbers  of  2 and  4;  and  on  Figures  88  to  103,  demonstrating 
the  results  of  air  injection  at  Mach  2.  To  provide  a datum  against  which  to  view  these  results,  the  Mach  2 
and  4 empty  tunnel  flows  are  presented  on  Figures  30  to  42.  Amongst  the  results,  we  find  static  pres- 
sure anr'  Mach  number  distributions  in  the  test  section;  streamwise  flow  and  crossflow  velocity  profiles 
in  three-dimensional  boundary  layers  along  with  their  integral  parameters,  surface  flow  visualization, 
and  the  distributions  of  local  skin  friction  coefficient  calculated  according  to  various  available  correla- 
tions for  Preston  tube  devices. 

Relevant  information  regarding  the  tunnel  settings  for  each  test  will  be  found  on  Table  5, 
while  integral  parameters  and  local  skin  friction  coefficients  from  the  circular  pitot  and  cobra  probe 
traverses  are  listed  on  Tables  6 and  7.  No  tabulations  of  point-by-point  data  through  the  viscous  flows 
will  be  found  herein,  but  are  available  from  the  NAE  upon  request. 

4 1 Accuracy  and  Calibrations 

4.1.1  Surface  Finish 

The  rms  surface  finish  on  the  test  wall  was  close  to  ten  micro  inches  with  surface  waviness 
better  than  0.0002-in.  per  in.  of  surface. 

4.1.2  Cobra  Probe  Pressures  and  Displacement  Effects 

The  pitot  pressures  collected  with  the  cobra  probe  were  accurate  to  within  ±0.04  psia  for  a 
typical  pitot  pressure  range  of  3 to  25  psia.  The  positioning  and  backlash  errors  in  the  probe  unit 
potentiometers  and  gearbox  were  assessed  by  means  of  a laboratory  workshop  comparator,  resulting 
in  translation  and  yaw  accuracies  within  ±0.0005-in.  and  ±0. 2-degree  respectively. 

On  the  two-dimensional  pitot  traverse  gear,  the  height  was  measurable  to  better  than  ±0. 001-in. 

The  flattened  three-tube  geometry  of  the  cobra  head  probe  results  in  a displacement  correc- 
tion that  is  somewhat  larger  than  that  normally  taken  for  circular  pitot  tubes85  in  incompressible  flow. 
Young  and  Maas86,  in  some  experiments  with  flattened  tubes,  made  a tentative  suggestion  of  an  outward 
displacement  correction  of  0.24  dc  towards  the  region  of  higher  velocity,  where  dc  is  the  width  of  the 
probe  face  in  the  direction  of  traverse.  This  value  was  used  in  experiments  utilizing  the  cobra  probe, 
that  were  reported  in  Reference  34  (for  h > 2 dc)  as  Marson  and  Lilley87  provide  some  evidence  that 
displacement  errors  are  not  affected  significantly  by  compressibility.  Since  the  publication  of  Reference  34, 
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however,  Ozarapoglu88  has  shown  that  the  displacement  effect  in  incompressible  flow  for  a flattened 
tube  of  similar  height  to  width  ratio  as  the  cobra  probe  is  closer  to  0.32  dc  for  h > 4 dc . For  distances 
closer  to  the  wall,  the  displacement  error  decreases  continuously  as  the  wall  is  approached,  and  Macmillan 
expressed  this  reduction  for  circular  pitots  in  terms  of  dc/h  and  uTw  dc/uw  . This  latter  quantity  in  the 
present  experiments  is  in  the  range  100  to  500.  Ozarapoglu’s  result  may  be  written  as 


A_ 

dc 


0.25 


where  the  second  term  is  included  only  when  h < 4 dc.  (For  the  circular  pitot  tube, 


A_ 

dc 


0.15  for 


h > 2 dc.)  For  distances  closer  to  the  wall  than  0.7  dc,  A/dc  was  negative  and  taken  to  be  -0.15,  a 
nominal  mean  value  of  the  scattered  data  of  Reference  88. 


We  should  comment  that  Allen  in  Reference  89,  reported  a circular  pitot  tube  displacement 
correction  found  in  a Mach  2 air  stream  that  was  very  much  larger  than  the  values  determined  in  incom- 
pressible flow.  His  assessment  yielded  A/dc  > 0.4  when  the  probe  was  away  from  the  wall.  On  this 
basis,  we  might  expect  the  rectangular  or  flattened  probe  to  yield  an  even  larger  correction.  However, 
since  the  substantial  body  of  data  so  far  is  towards  somewhat  lower  corrections  than  Allen’s*,  and 
since  Hall  and  Dickens90,  when  using  a similar  probe  form  to  that  of  the  cobra  found  a displacement 
correction  to  be  no  larger  than  0.15  dc  in  a Mach  1.6  to  2.0  airflow,  we  have  chosen  to  use  the  above 
equation  as  a reasonable  compromise.  In  fact,  the  displacement  correction  of  ^0. 005-in.  away  from 
the  wall  is  virtually  indiscernible  on  the  velocity  plots. 

No  corrections  for  the  effect  of  turbulence  on  the  pitot  readings  have  been  employed. 

4.1.3  Static  Hole  Errors 

Errors  due  to  large  static  pressure  holes  (i.e.  relative  to  the  boundary-layer  displacement 
thickness  6 , *)  can  be  important  in  the  reduction  of  boundary-layer  profile  measurements  and  when 
using  surface  Preston  tubes  to  determine  surface  shear  stress.  Rainbird91  found  that  in  subsonic  and 
supersonic  zero  pressure  gradient  flow  at  high  Reynolds  numbers,  large  static  orifices  (where  ds/5 1 * = 10) 
gave  pressure  errors,  AP,  typically  twenty  times  the  value  of  the  undisturbed  surface  shear  stress.  In  the 
present  experiments,  the  parameter  ds/6,  * was  < 1,  so  that  the  influence  of  ds/5,  * was  not  significant. 

In  fact,  we  may  compute  the  correction  to  static  pressure,  using  the  convenient  correlation  of  Franklin 
and  Wallace92  for  deep  static  holes,  but  extended  '■nmnressible  flow  by  inserting  the  local  wall  con- 
ditions of  density  and  temperature.  The  following  curve  fit  was  cho&cr.  to  represent  the  correlation92, 
namely 
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Wilson99,  in  fact,  pointed  out  that  the  displacement  effects  quoted  by  Allen  were  unnecessarily 
large  and  that  they  were  applicable  only  to  his  particular  pitot  tube  and  support  configuration. 
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The  functional  relationship  in  the  above  expression  follows  from  the  law  of  the  wall  and  the  assumption 
that  the  static  pressure  hole  error  will  be  a function  of  the  mean  velocity  u at  a distance  h from  the  wall, 
which  is  dependent  on  ds9 1 . Provided  a law  of  the  wall  region  exists  in  adverse  pressure  gradient  flow, 
we  may  argue  that  the  hole  error  correction  should  be  applicable  there,  as  in  the  zero  pressure  gradient 
flow  case,  noting  that  local  conditions  at  the  wall  should  be  used  in  the  hole  error  formula. 

In  the  present  experiment,  the  0.015-in.  diameter  static  orifice  yields  a static  hole  error  of 
amount  typically  twice  the  local  shear  stress.  This  value  is  the  approximate  uncertainty  in  measuring 
static  pressure,  and  hence  has  not  been  subtracted  from  the  experimental  measurement. 

4.1.4  Tunnel  Mainstream  Flow 

The  inaccuracies  in  measurement  of  tunnel  flow  quantities  and  repeatability  indicate  that 
the  maximum  likely  errors  in  local  Mach  number  and  skin  friction  are  AM  = ±0.02  and  ACp  = ±0.0001. 

The  tunnel  stagnation  pressure  is  controlled  typically  to  within  ±0.1  percent  of  the  mean 
stagnation  pressure  level  during  any  given  run,  which  is  to  within  ±0.03  psi  at  Mach  2 and  ±0.2  psi  at 
Mach  4. 

4.1.5  Tunnel  Noise 

Intense  pressure  fluctuations  have  been  found93  in  the  settling  chamber  of  the  5 X 5-in. 
blowdown  wind  tunnel.  These  fluctuations  are  a function  of  the  flow  conditions  at  the  control  valve, 
implying  that  they  are  aerodynamic  noise  originating  from  the  jet  of  the  control  valve.  Below  a working 
section  Mach  number  of  3,  the  fluctuations  in  the  working  section  were  found93  to  be  predominantly 
those  which  had  been  transmitted  from  the  settling  chamber.  The  attenuation  of  the  transmission  in- 
creased as  the  nozzle  throat  area  closed,  and  at  Mach  numbers  greater  than  3,  the  working  section  fluc- 
tuations were  almost  independent  of  the  settling  chamber  fluctuation  level. 

At  Mach  4,  the  overall  sound  pressure  level  in  the  working  section  20  log10  (A p/0. 0002 
dynes/cm2),  was  146*  dB93  at  the  stagnation  pressure  level  used  in  the  present  experiment.  At  Mach  2, 
on  the  other  hand,  where  the  pressure  drop  across  the  control  valve  was  greater  than  four  (a  feature  of 
our  present  tests),  the  overall  sound  pressure  level  increased  to  158  dB*.  This  is  a substantial  fluctuation 
level,  and  Green94  argues  that  such  noise  might  artificially  increase  the  effective  Reynolds  number  of 
the  flow  by  significant  amounts. 

If  we  assume  acoustic  waves  propagating  parallel  with  a uniform  flow,  then  the  ratio  of  pres- 
sure and  velocity  fluctuations  in  the  flow  is 

Ap  = p a Au 

The  local  speed  of  sound  is  denoted  by  a,  p is  the  mean  density  and  the  Ap,  Au  are  the  fluctuating  pres- 
sure and  velocity  terms.  Now  the  ratio  of  Ap  to  p(Au)2  is  typically  very  large  in  a sound  field.  In 
contrast,  the  same  ratio  becomes  small  for  fluctuations  of  spatially  homogeneous  and  isotropic  turbu- 
lence95 , with  the  result  that  Ap  ~ p(Au)2  • If  we  accept  the  former  definition  as  more  representative  of 
the  blowdown  wind  tunnel  flow,  especially  at  Mach  2,  the  overall  sound  pressure  level  of  158  dB  in  the 
working  section  corresponds  with  a fluctuation  velocity  of  about  1.4  percent  of  the  local  mainstream 
velocity,  which  although  being  an  upper  bound  estimate,  indicates  a substantial  fluctuation  level. 

As  Mabey  discusses  in  Reference  96,  the  spectrum  of  surface  pressure  fluctuations  for  a 
nominal  2D  turbulent  boundary  layer  approaching  separation  may  be  divided  into  high-frequency  and 
low-frequency  components.  The  high-frequency  fluctuations  are  much  the  same  as  those  measured 
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The  two  dB  levels  quoted  correspond  with  rms  fluctuating  pressure  levels  at  Mach  4 and  Mach  2 
of  0.03  and  0.7,  expressed  as  percentages  of  the  tunnel  flow  stagnation  pressure. 


under  a boundary  layer  in  zero  pressure  gradient,  and  are  generated  in  the  small-scale  inner  region  of 
the  boundary  layer  associated  with  the  law  of  the  wall.  On  the  other  hand,  the  low-frequency  pressure 
fluctuations  are  generated  in  the  large-scale  outer  region.  It  has  been  shown  by  McDonald1 09 , and  is 
further  discussed  by  Zwarts97 , that  increase  of  the  mainstream  turbulence  level  does  produce  higher 
shear  stresses  in  the  outer  region  and  does  change  the  mean  velocity  profile.  In  addition,  the  measure- 
ments of  Fiedler  and  Head98  demonstrate  that  pressure  gradient  has  a large  effect  on  the  distribution 
of  intermittency  in  the  outer  region  of  a turbulent  boundary  layer.  As  the  turbulent  boundary  layer 
approaches  separation  and  the  shape  parameter7 , 5,  */0n  , decreases,  the  intermittent  zone  decreases 
in  width  and  moves  farther  from  the  surface.  The  structure  of  a turbulent  boundary  layer  approaching 
separation,  then,  does  not  approximate  to  that  of  a free  wake  flow,  the  latter,  which  in  contrast,  is 
characterized  by  a very  wide  spread  of  intermittency. 

The  measurements  of  Fiedler  and  Head  would  seem  to  indicate  that  an  inner  region  does 
exist  as  the  boundary  layer  approaches  separation,  even  though  the  log  law  region  appears  to  shrink 
with  increase  of  adverse  pressure  gradient.  The  congruent  laws  of  the  wall  found  for  both  a turbulent 
boundary  layer  and  a fully  developed  pipe  flow  demonstrate  that  the  degree  of  intermittency  in  the 
outer  region  does  not  significantly  affect  the  mean  flow  of  the  small  scale  inner  region.  If  this  follows 
for  the  2D  turbulent  boundary  layer  approaching  separation,  we  might  expect  the  same  for  the  3D 
case,  where  the  flow  properties  of  the  inner  region  should  still  be  scaled  by  local  conditions  at  the  wall. 
Hence  we  would  expect  the  onset  of  separation  or  incipient  separation  to  be  dominated  by  the  same 
wall  properties. 

4.1.6  Cobra  Probe  Thermocouple  Measurements 

The  cobra  probe  assembly  was  placed  in  still  air  in  a controlled  oven  enclosure,  to  first  of  all 
check  the  integrity  of  the  iron/constantan  thermocouple  bead  circuit  and  its  junction  constant,  meas- 
ured in  degrees  F per  millivolt  output. 

The  probe  was  fitted  in  the  5 X 5-in.  blowdown  wind  tunnel  at  its  maximum  extension  of 
0.8-in.  (and  hence  in  the  mainstream).  The  oncoming  Mach  number  to  the  probe  was  then  changed 
by  adjusting  the  deflection  of  the  shock  generator  wedge  that  was  installed  upstream  of  the  probe 
station.  The  probe  was  allowed  to  align  itself  with  the  local  mainstream  velocity  vector. 

The  true  stagnation  temperature  was  read  on  the  glass  bulb  thermometer  in  the  settling 
chamber,  along  with  the  temperature  diode  in  the  same  location.  In  this  manner,  the  recovery  factor 
was  determined  at  Mach  numbers  between  1.6  and  2.0,  and  from  2.7  to  4.0,  at  Reynolds  numbers 
based  on  the  half-height  of  the  probe  and  the  local  stream  conditions,  in  the  range  6 to  12  X 104 . 

The  recovery  factor,  R,  of  the  probe  thermocouple,  is  defined  as  the  ratio  of  the  indicated 
temperature  rise,  TPM  - T,  to  the  adiabatic  temperature  rise,  T0  - T;  that  is 


T 


Rt  is  the  ratio  of  the  measured  total  temperature  at  the  thermocouple,  Tpm  , to  rue 

stagnation  temperature,  T0 , and  T is  the  local  static  temperature.  Clearly,  R and  RT  are  r<  h 

the  local  stream  Mach  number,  and  Figure  28  illustrates  this  relationship. 

Within  our  limits  of  measurement,  we  found  that  RT  = 0.981  ± O.OO4  for  Mach  l.i  .j  2.U, 
and  0.971  ± O.O2  from  Mach  2.7  to  4.0.  The  corresponding  recovery  factors  were  0.957  and  0.962 
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respectively,  being  slightly  less  than  the  factors  found  by  Bradfield  et  al81 , Figure  29,  and  due,  presum- 
ably, to  enhanced  conduction  losses  from  the  present  design.  There  was  no  conclusive  effect  of  Reynolds 
number  in  the  range  investigated. 

Once  the  probe  is  calibrated,  the  local  static  temperature  in  an  unknown  flow  may  be  written 
as 


T = 


provided  R does  not  vary  significantly  through  the  range  of  Mach  number  and  Reynolds  number  to  be 
encountered  in  the  new  flow  field. 

(In  a given  traverse  of  the  cobra  probe  at  Mach  2,  for  example,  the  excursion  in  total  temper- 
ature, Tpm  , measured  through  the  boundary  layer  was  usually  8-10  degrees  F,  while  superimposed  on 
the  mean  temperature  profile  we  observed  t a maximum  peak-to-peak  fluctuation  of  about  ±C.3  degree  F.) 

The  resultant  velocity  in  our  skewed  boundary-layer  flow  is 
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noting  7 = 1.4  and  the  boundary-layer  edge  is  at  the  temperature  of  the  mainstream. 

Alternatively,  in  terms  of  the  temperature  ratio,  where  TPC  is  the  corrected  probe  temperature, 
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4.1.7  Wall  Thermocouple  Measurements 


At  the  same  time  as  checking  out  the  cobra  probe  thermocouple,  the  wall  temperature  thermo- 
couple plugs  (Fig.  23)  were  also  placed  in  the  controlled  environment  of  a warm  oven  to  check  their 
circuit  continuity.  The  thermocouple  junction  constant,  measured  in  degrees  F per  millivolt  output, 
was  then  used  in  determining  subsequent  wall  temperature  measurements  with  the  flow  on.  As  the 
circuitry  of  the  temperature  measuring  system  for  the  cobra  probe  was  used  also  for  the  wall  thermo- 
couples, the  accuracy  of  measurement  as  mentioned  in  the  last  subsection  should  apply  here. 


Since  the  thermal  inertia  of  the  thin  flat  plate  surface  element  should  not  be  significant  on  the 
time  scale  of  the  present  measurements  (=  1 min.),  it  should  take  up  the  wall  temperature  relatively  soon 
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? r , T flOW-  We  might  dso  exP«*  that  temperature  gradients  across  the 
i ^ Sted  SUlfaCe  element  would  be  smaU  ^r  seconds,  or  minutes  ofTnnel  blowdown 
so  that  the  measurement  of  temperature  on  the  inner  surface  of  the  steel  element  should  yield  a value 

element  * **  N°  corTections  to  account  for  the  finite  thickness  of the  surface 

element  such  as  those  discussed  by  Naysmithi  >0  and  Harris82  have  been  included 


4.2  Data  Reduction 


4.2.1  Quantities  in  the  External  Flow  at  the  Tunnel  Reference  Station  1 

If  we  assume  that  the  expansion  of  the  flow  from  the  settling  chamber  to  the  working  section 

oTSon  pressure  raS  “S*  JOUndary‘layer  «>ge,  ME1 , is  given  by  the  local  static 
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where  P0  is  the  stagnation  and  pj  the  static  pressures. 
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and  aE  ( is  the  local  speed  of  sound. 
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so  that  the  local  unit  Reynolds  number  per  ft  may  be  computed,  where 
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Allowing  a recovery  factor  of  0.89  at  the  wall,  the  adiabatic  wall  temperature  is 


0.89  ug , 


(4.2.7) 


where  cp  is  the  specific  heat  at  constant  pressure  = 6006  ft2  /sec.2  °R. 

42.2  Measurements  Through  the  Boundary  Layer 

The  local  wall  static  pressure,  p,  measured  in  the  absence  of  the  probe  measuring  device,  was 
assumed  to  be  constant  along  a normal  to  the  wall  at  all  measurement  stations  through  the  swept  inter- 
action region.  Static  pressure  through  a boundary  layer,  where  shock  waves  are  present,  is  a difficult 
quantity  to  measure  with  a conventional  static  pressure  probe5 , for  the  interaction  between  the  shocks 
and  the  local  boundary  layer  on  the  probe  itself  will  lead  to  spurious  pressures  at  the  static  measuring 
orifices.  For  example,  in  the  very  complex  mixing  between  a wall  jet  and  a separated  shock/turbulent 
boundary- layer  interaction  region  at  Mach  1.8,  a static  pressure  variation  at  the  most  of  about  ten  per- 
cent of  the  wall  static  pressure  was  indicated  by  a cone/cylinder  type  static  pressure  probe5 . In  our 
swept  interaction,  except  when  large-scale  3D  separation  is  present,  we  might  expect  variation  from 
the  wall  to  the  mainstream,  to  be  substantially  less  than  four  percent  of  the  local  dynamic  pressure  at 
the  outer  edge  of  the  boundary  layer. 

The  wall  static  pressure  was  then  combined  with  the  measured  values  of  circular  or  cobra 
probe  pitot  pressure,  Pp , to  form  the  ratio  of  Pp  /p  from  which  the  local  Mach  number  in  the  flow  was 
computed.  The  subsonic  Mach  numbers  (Pp  /p  < 1.893)  were  calculated  from  the  isentropic  flow  equa- 
tion 


7 

PP  / 7 _ 1 2 \ ^ ” * / 7 \ 

— = ( 1 + 2~y~  M2  j = (l  + 0.2  M2  ) 


(4.2.8) 


whence  Mach  number  is  obtained  explicitly. 

The  supersonic  Mach  numbers,  where  Pp  /p  > 1.893,  were  calculated  via  an  iteration  routine 
from  the  Rayleigh  pitot  formula 
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The  resultant  velocity  with  respect  to  that  at  the  boundary-layer  edge  is 
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where  Ttot  is  the  local  total  temperature  in  the  boundary  layer,  where 
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When  using  the  circular  pitot  tube  in  the  empty  tunnel  boundary  layers,  the  distribution  of 
total  temperature  across  the  boundary  layer  was  assumed  constant  and  equal  to  the  settling  chamber 
value,  then 
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On  the  other  hand,  the  cobra  probe  thermocouple  provided  the  facility  to  measure  the  total 
temperature  distribution  through  the  boundary  layer  in  the  swept  interaction  region.  The  corrected 
probe  temperature,  TPC,  (see  Section  4.1)  is  now  equivalent  to  Ttot  in  the  equation  for  local  velocity, 
so  that 
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the  total  temperature  at  the  boundary-layer  edge  being  typically  no  greater  than  the  stagnation  tem- 
perature in  the  tunnel  settling  chamber,  T0*,  by  more  than  0.2  percent. 

The  mean  velocity  components  within  the  boundary  layer  in  streamwise  co-ordinates,  that 
is,  with  respect  to  the  local  flow  direction  at  the  boundary  layer  edge  are 
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for  the  streamwise  flow;  and  for  the  crossflow, 
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where  pCL  is  the  measured  angle  of  yaw  between  the  resultant  velocity  vector  at  the  cobra  probe  pitot 
orifice  and  the  (datum)  tunnel  centre-line  direction. 

Alternatively,  the  corresponding  local  velocity  components  through  the  boundary  layer, 
parallel  and  perpendicular  to  the  tunnel  centre-line  (the  nominal  undisturbed  mainstream  direction) 


Ul  UB 

— = COS  PCJ 

UE  UE 


(4.2.17) 


* A single  point  measurement 
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vl  UB 
— = — sin  pCL 

uE  uE 


(4.2.18) 


With  constancy  of  static  pressure  assumed  through  the  viscous  flow,  we  have  for  the  ratio  of 
local  density  in  the  boundary  layer  to  that  at  the  boundary -layer  edge 


that  is. 


p _ T E _ / Te\  /Tpce\  /Tpc\ 
pE  T \Tpce/  \ Tpc  / \ T/ 

_p_  _ /l  + 0.2  M2  \ Tpce 
pE  = ^1  + 0.21^  TpC 


(4.2.19) 


(4.2.20) 


4.2.3  Integral  Parameters 

The  experimental  values  of  stream  wise  and  crossflow  displacement  and  momentum  thickness 
parameters  were  calculated  according  to  the  expressions  in  the  list  of  symbols,  by  performing  numerical 
integrations  of  the  traverse  data  using  the  trapezoidal  rule.  No  special  curve  fitting  was  employed  be- 
tween the  wall  and  the  first  data  point.  As  we  see  in  the  list  of  symbols,  the  integrands  of  the  various 
integral  parameters  are  composed  of  combinations  of  the  terms  p/pE , u/uE  and  v/uE , so  that  by  recourse 
to  Equations  (4.2.14),  (4.2.15),  (4.2.16)  and  (4.2.20),  we  may  calculate  all  of  the  chosen  integral 
parameters. 

We  note  that  in  a 3D  boundary-layer  flow  the  thickness  of  the  true  displacement  surface 
cannot  be  related  directly  to  an  integral  of  the  mean  velocity  profile  at  one  position  on  the  test 
wall1 1 7-1 1 8 . Instead,  it  can  be  determined  from  the  normal  component  of  velocity  of  the  boundary- 
layer  edge  by  means  of  a partial  differential  equation:  see  the  calculation  method  of  Nash  in  Appendix  1 
(Eqn.  A.1.16). 

4.2.4  Determination  of  Local  Skin  Friction  Coefficients  Using  Preston-Tube  Correlations 

The  flow  in  the  inner-layer  (or  wall  law  region)  of  an  incompressible  turbulent  boundary  layer 
is  determined  entirely  by  the  local  wall  shear  stress  and  fluid  properties;  and  is  not  influenced  by  the 
conditions  further  from  the  surface  except  insofar  as  these  determine  the  skin  friction.  Dimensional 
analysis  then  shows  that  the  only  velocity  scale  is 


- ■ (j) 


(4.2.21) 


and  that  the  only  length  scale  is  u/uTW . Hence,  as  argued  first  by  Preston,  the  difference  between  sur- 
face pitot  pressure  and  the  local  static  pressure  (in  the  absence  of  the  pitot)  may  be  expressed  in  such 
forms  as 


or  more  conveniently 


= f /urWdc\ 

\ V / 

rwdc2  ^ / Apd2  \ 

pv 2 F \ pv2  ) 


(4.2.22(a)) 


(4.2.22(b)) 
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where  dc  is  the  outer  diameter  of  the  (circular)  surface  pitot  tube,  and  Ap  is  the  difference  in  pressure 
between  the  surface  pitot  and  a local  static  pressure,  measured  in  the  absence  of  the  pitot. 


In  compressible  flow  over  adiabatic  walls,  the  local  speed  of  sound  at  the  wall  is  a suitable 
parameter  as  a measure  of  the  compressibility  effects.  The  corresponding  dimensionless  parameter  is 
called  the  ‘friction  Mach  number’,  MT  w t where 

urW 

MtW  = (4.2.23) 

aw 

and  if  we  define  the  fluid  properties  at  the  wall, 


Ap 

Tw 


/UrWdc 
\ °W 


(4.2.24) 


is 


The  best  form  of  calibration  law  for  Preston  tubes  in  compressible  flow,  according  to  Bradshaw 1 0 1 , 


Ap 

Tw 


+ f. 


/UrWdc 
l °W 


(4.2.25) 


where  fj  (urWdc/uw  ) is  the  well-proven  formulation  of,  say,  Patel107  in  low  Mach  number  flow,  and  fc 
is  a compressibility  correction  function.  To  determine  the  form  of  fc , Bradshaw  re-plotted  Allen’s  data 
of  Reference  119  as  Ap/rw  versus  utV  dc/uw  for  various  values  of  MT  w , and  demonstrated  that  an  ade- 
quate correlation  was 


fc 


104  mtI 


2.0 


(4.2.26) 


in  the  range  50  < urWdc/uw  < 1000,  0 < MtW  < 0.1.  Bradshaw  thought  it  reasonable  to  claim  an 
accuracy  of  ±10  percent  for  fc  in  the  range  stated.  The  error  in  Patel’s  data  for  fj,  on  the  other  hand, 
was  about  ±2  percent  or  less,  so  that  using  Equation  (4.2.25)  with  (4.2.26),  we  should  expect  maximum 
errors  in  skin-friction  of  ±10  percent  for  the  above  ranges  of  uTWdc/uw  and  MtW  . 


Combining  Patel’s  and  Allen’s  results,  Bradshaw  wrote  an  implicit  equation  to  determine  the 
wall  shear  stress  as 


9 


in  the  ranges  50  < urWdc/uw  < 1000,  0 < MrW  < 0.1.  The  equation  was  left  deliberately  in  this  com- 
plicated form  to  facilitate  modifications  to  the  compressibility  factor,  or  extension  to  Patel’s  low-speed 
calibration  at  higher  Reynolds  numbers,  (see  Ref.  88).  Moreover,  the  calibration  law  expressed  in  Equa- 
tion (4.2.27),  will  clearly  only  apply  to  air  flow  over  adiabatic  walls,  for  heat  transfer  will  introduce  the 
additional  parameter  Hw  /p  w cp  Tw  uT  w , where  Hw  is  the  heat  flux  per  unit  area  of  wall.  At  this  point 
in  time,  we  have  no  reliable  theory  for  the  effect  of  heat  transfer  on  Preston  tube  calibrations. 


w 


I 


Some  other  less  recent  Preston  tube  correlations  have  also  been  used,  to  compare  with  the 
results  obtained  from  Equation  (4.2.27).  These  are  the  correlations  of  Hopkins  and  Keener1 00,  where 
in  one,  an  ‘intermediate  temperature’,  T',  and  in  the  second,  the  local  wall  temperature,  Tw , is  used. 

It  will  be  remembered  that  the  intermediate  temperature  hypothesis  proposes  that  one  relation  for  sur- 
face shear  stress  may  be  found  for  both  incompressible  and  compressible  flows  if  the  density  and  kine- 
matic viscosity  are  calculated  at  T',  a value  between  die  wall  and  external  stream  static  temperature. 
The  form  of  T'  given  in  Reference  120  has  been  used  herein 


T'  = T, 


1 + 0 035  M|  + 0.45  1 


(4.2.28) 


where  the  wall  temperature,  Tw  " Taw.  the  adiabatic  wall  temperature,  with  recovery  factor  of  0.89. 
As  an  alternative  representation  of  Equation  (4.2.7),  we  may  write 


Taw  = TE(l  + 0.178  m|) 


(4.2.29) 


Hopkins  and  Keener  arrived  at  their  correlations  by  collapsing  the  compressible  flow  Preston 
tube  data  onto  a relevant  incompressible  correlation 


Rd,  cpp.  - f^Cp.] 


(4.2.30) 


by  writing 


RDcCPPc  - 'KcFJ 


(4.2.31) 


the  subscripts  i and  c referring  to  incompressible  and  compressible  flow,  respectively.  RD  is  the 
Reynolds  number  based  on  the  Preston  tube  dimension  dc  (the  outer  diameter  of  the  forward  facing 
circular  wall  pitot,  for  example),  Cpp  is  the  Preston  tube  pressure  coefficient,  and  CF  is  the  local  skin 
friction  coefficient. 


The  resulting  functional  equations  for  compressible  flow  in  the  form  of  Equation  (4.2.31) 


Method 


Hopkins/Keener, 

T',  (Ref.  100,  Eq.  4) 

Hopkins/Keener, 

Tw,  (Ref.  100,  Eq.  9) 


Calibration  Factors 

loglot^Dj.  CppJ 

m log10[RDcCFc  ]+Constant 

log10[f2(T')R^(Mp/ME)2] 

1.132  log10[f2(T')RpCF  1+1.517 

loglol^fTw  )^D  Cpp  ] 

1.162  logj  0[f2(Tw  )RpCF  ]+1.552 

Equation 
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where,  with  T in  °R, 


f2(T')  = 


~PE‘ 

2 , 
P_ 

Y + 198.6  " 

2 

"TE" 

t 

_P 

Pe 

Te  + 198.6 

T' 

(4.2.34) 


Rd  = 


xEdc 


°E 


(4.2.35) 


and  Mp  is  the  Mach  number  calculated  from  the  local  static  pressure  and  the  Preston  tube  pitot  pressure, 


via  the  Rayleigh  pitot  equation. 


To  provide  a comparison  with  the  methods  using  the  wall  pitot  pressure,  the  Sommer/Short 
correlation1 20  was  also  tried,  in  which  we  used  the  momentum  thickness  Reynolds  number  of  the 
boundary-layer  profiles,  based  on  the  streamwise  momentum  thickness,  du.  The  Sommer/Short 
correlation  also  utilizes  the  intermediate  temperature  and  is  written 


Cp  - 


17.076 


T 2 


+ 25.112  log! 


»(“•"  7) 


(4.2.36) 


+ 6.012 


In  the  three-dimensional  boundary  layer,  we  have  assumed  that  the  above  2D  correlations 
(using  wall  pitot  pressure)  hold  for  the  cobra  probe  geometry,  provided  the  wall  pitot  is  aligned  with 
the  local  shear  stress  trajectory  at  the  wall.  Further  comments  regarding  the  measurements  from  the 
circular  and  cobra  surface  pitots  are  left  until  Section  4.3,  when  the  empty  tunnel  flow  is  under  dis- 
cussion. Moreover,  all  of  the  above  correlations  were  developed  in  nominally  zero  pressure  gradient, 
2D  flat  plate  turbulent  boundary-layer  flows.  The  effect  of  pressure  gradients  on  Preston  tube  calibra- 
tions in  incompressible  flow  was  discussed  by  Patel107,  where  the  general  conclusion  was  that  provided 
a logarithmic  law  of  the  wall  region  existed,  the  calibrations  could  be  used  as  long  as  the  Preston  tube 
was  suitably  restricted  in  size.  In  addition,  the  more  recent  calculations  of  surface  shear  stress  by 
Bradshaw,  to  compare  with  experimental  flows  with  significant  adverse  pressure  gradients  at  Mach  4, 
and  in  which  the  effects  of  bulk  compression  and  dilatation  are  included,  show  very  favourable  agree- 
ment with  Preston  tube  results35 . Since  also,  the  outside  vertical  dimensions  of  the  circular  Preston 
tube  and  rectangular  pitot  tube  of  the  cobra  probe  are  well  within  the  maximum  and  minimum  toler- 
ances on  size  specified  by  Hopkins  and  Keener1 00  for  good  operation,  we  believe  that  the  two  pitot 
geometries  used  to  measure  surface  skin  friction  in  our  empty  tunnel  and  swept  interaction  flows 
provided  satisfactory  results. 


4.3  Flow  in  the  Empty  Tunnel  at  Mach  2 and  4 
4.3.1  Surface  Pressures  and  Flow  Visualization  at  Mach  2 


1 


Some  initial  tests  in  the  empty  wind  tunnel  were  conducted  to  determine  the  two-dimensionality 
of  the  undisturbed  test  boundary  layer  on  the  tunnel  floor  at  nominal  mainstream  Mach  numbers  2 and  4. 
The  uniformity  of  the  flow  was  assessed  by  reviewing  wall  static  pressures,  surface  oil  flow  visualization, 
mean  velocity  profiles  at  several  lateral  positions  and  at  two  axial  stations  in  the  working  section,  and 
local  surface  skin  friction. 


Figure  30  shows  surface  pressure  and  oil  dot  flow  visualization  results  at  Mach  2,  with  a 
stagnation  pressure  of  the  mainstream  flow  equal  to  34  psia.  This  value  of  blowing  pressure  was  used 
for  all  of  the  Mach  2 traverse  measurements  through  the  swept  interaction  to  be  discussed  in  subsequent 
sections.  The  surface  pressures  were  made  non-dimensional  by  reference  to  the  static  pressure  at 
Station  1 in  the  upstream  part  of  the  working  section  as  shown  on  Figure  15.  On  the  top  half  of 
Figure  30,  the  surface  pressures  and  oil  flow  visualization  indicate  reasonable  uniformity  across  the 


J 


-32- 


width  of  the  test  section.  Reference  to  Figure  15  will  assist  in  placing  the  0.5-in.  apart  lines  of  static 
pressure  holes  (labelled  Y,  Z,  <£,  A,  B,  D and  E)  on  the  test  floor,  and  the  two  axial  stations  numbered  4 
and  18  where  pitot  traverses  were  taken.  We  should  note  that  the  origin  of  the  X-ab6cissa  is  arbitrarily 
set  at  the  position  of  the  turntable  centre-line  of  the  blowing  blc  experiment  that  we  see  on  Figure  24  — 
this  is  for  corvenience  in  referencing  the  later  data  with  blowing.  Such  origin  is  situated  8.625-in.  up- 
stream of  the  exit  of  the  tunnel  working  section. 

Small  ‘wiggles’  in  the  floor  pressure  distributions  across  the  tunnel  flow  at  X s 4-in.  would 
appear  to  be  caused  by  a compression-expansion  wave  disturbance  from  the  upstream  edge  of  a non- 
flush fitting  of  the  ceiling  window.  The  corresponding  maximum  peak-to-peak  Mach  number  deviations 
are  within  the  spread  AM  = ±0.022.  The  oil  flow  visualization  indicates  some  small  convergence  towards 
the  tunnel  centre-line  — a residue  from  the  lateral  pressure  gradient  in  the  nozzle,  but  the  overall  paral- 
lelism may  be  deemed  quite  satisfactory. 

The  effect  of  the  boundary  layers  on  the  tunnel  walls*  perpendicular  to  the  test  surface  is 
not  usually  significant  if  a weak  adverse  pressure  gradient  is  superimposed  on  the  floor  test  boundary 
layer.  With  increase  of  adverse  pressure  gradient,  such  as  will  occur  when  we  deflect  the  wedge  (Fig.  13) 
more  than,  say  6-degrees,  the  degree  of  influence  of  these  ‘boundary  constraints’  will  increase,  to  even- 
tually affect  the  test  flow  in  a gross  manner  when  strong  separations  are  present.  One  way  of  alleviating 
the  endwall  constraint  is  to  partially  remove  the  boundary  layer,  but  this  is  not  usually  practical  in  a 
supersonic  facility  where  local  static  pressures  within  the  flow  are  low.  An  alternative  means  of  removal 
is  by  a sweptback  boundary -layer  fence,  or  vortex  generator  (v.g.),  submerged  within  the  boundary  layer, 
whose  height  is  roughly  0.35.  We  experimented  with  a vortex  generator  whose  leading-edge  was  set 
close  to  the  intersection  of  the  floor  and  the  wall  of  the  tunnel  (facing  the  deflected  surface  of  the 
wedge)  at  XVGLE  = - 6.25-in.  and  0.0-in.,  where  the  angle  of  the  v.g.  to  the  floor  was  20-degrees  i.e., 
a sweepback  angle  of  70-degrees.  The  lower  illustrations  on  Figure  30  show  the  results  for  the  more 
upstream  v.g.  position. 

The  pressure  distributions  along  lines  Y to  B show  no  effect  of  the  installation  of  the  v.g. 

Only  along  lines  D and  E,  closes  to  the  comer,  is  there  seen  evidence  of  the  v.g.,  which,  however,  is  not 
evident  on  the  floor  oil  dot  flow  visualization  between  Stations  4 and  18.  The  mechanism  of  the  v.g. 
is  to  cause  a three-dimensional  separation  of  the  oncoming  boundary  layer  which  rolls  up  along  and 
from  the  sharp  free  edge  of  the  vertical  fence,  and  passes  downstream  in  the  comer  between  the  wail 
and  the  tunnel  ceiling.  In  so  doing,  a new  boundary  layer  commences  from  the  attachment  line  down- 
stream of  the  v.g.  to  produce  a substantially  reduced  local  boundary-layer  thickness  on  the  endwall  in 
the  region  of  the  test,  and  a consequent  reduction  in  induced  upstream /downstream  lengths  of  reflected 
shock/boundary-layer  interaction.  The  lower  right-hand  photograph  of  Figure  30  illustrates  a large 
divergence  in  the  limiting  streamline  pattern  on  the  tunnel  endwall  just  downstream  of  the  v.g.,  and 
we  interpret  this  as  the  attachment  region  where  the  growth  of  the  ‘new’  endwall  boundary  layer 
commences. 

At  Mach  4,  the  empty  tunnel  flow  on  the  test  floor  is  disappointingly  non-uniform  — a result 
that,  in  retrospect,  might  have  been  anticipated!  Figure  31,  on  the  top  left-hand  side,  shows  large  per- 
turbations along  all  lines  of  static  pressure  orifices  away  from  the  tunnel  centre-line.  This  is  clearly  the 
result  of  using  the  tunnel  nozzle  sidewall  as  the  test  floor,  where  the  transverse  pressure  gradients  in  the 
acceleration  from  the  tunnel  throat  to  the  Mach  4 working-section,  cause  the  test  boundary  layer  to 
converge  towards  the  tunnel  centre-line,  as  we  see  on  the  flow  visualization  on  the  top  right-hand  side 
of  Figure  31. 


* Because  of  the  orientation  of  the  wind  tunnel  nozzle  assembly,  see  Figure  13,  there  may  be  some 
confusion  with  respect  to  terminology  used.  The  surface  or  wind  tunnel  ‘floor'  on  which  we  have 
the  test  boundary  layer  is,  in  fact,  the  sidewall  of  the  nozzle.  The  profiled  nozzle  blocks  provide 
the  vertical  boundaries  to  the  test  section  flow.  In  the  ensuing  discussion,  these  vertical  walls  will 
< be  denoted  the  endwalls  or  sidewalls  to  the  test  floor. 
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In  an  attempt  to  alleviate  the  non-uniformity,  we  investigated  the  installation  of  a V-type 
vortex  generator,  of  height  typically  equal  to  0.36,  but  this  time  situated  symmetrically  on  the  test 
floor,  upstream  of  the  test  section.  A sketch  of  the  position  of  the  v.g.  is  shown  on  the  lower  leftrhand 
side  of  Figure  31,  where  the  semi-apex  angle  was  30-degrees  (i.e.  each  v.g.  was  swept  back  at  60-degrees). 
The  v.g.  apex  was  situated  on  the  tunnel  centre-line,  and  each  swept  ‘arm’  of  the  v.g.  was  sufficiently 
long  to  span  the  floor  between  the  centre-line  and  the  endwall. 

With  the  vortex  generator  in,  we  note  from  the  illustrations  on  the  lower  half  of  Figure  31, 
a gross  improvement  in  the  smoothness  of  the  surface  pressure  distributions,  and  a change  from  overall 
convergence  of  the  limiting  streamlines  towards  the  centre-line,  to  one  of  slight  divergence  in  the  centre- 
region  of  the  flow.  A hump  in  the  pressure  distribution  appears  on  the  centre-line  near  X = 0.0-in.,  the 
reason  for  which  is  unclear.  The  correspondinb  perturbation  in  Mach  number  is  AM  = 0.091.  Towards 
the  endwalls,  the  oil  flow  visualization  allows  the  induced  effects  of  the  rolled-up  vortical  fluid  from 
the  V-type  vortex  generator  to  be  detected  as  local  ripples  in  surface  shear  stress  direction. 

The  mode  of  operation  of  the  vortex  generator  is  to  again  cause  a three-dimensional  separa- 
tion, but  on  this  occasion,  of  the  test  boundary  layer  upon  the  floor.  The  thicker  centre-line  region  is 
caused  to  flow  outwards  towards  the  endwalls  by  means  of  the  induced  separated  flow,  to  provide  some 
improvement  in  uniformity  of  the  lateral  distribution  of  boundary-layer  thickness. 

On  the  lower  right-hand  side  of  Figure  31,  we  see  small  horse-shoe  type  three-dimensional 
separations  generated  by  the  Preston  tubes  at  Stations  4 and  18;  and  the  local  effects  of  the  v.g.  <1ow 
upon  the  limiting  streamlines  on  the  starboard  endwall. 

The  stagnation  pressure  level  was  about  170  psia  for  these  empty  tunnel  Mach  4 tests,  and 
this  value  was  used  in  the  cobra  probe  traverses  to  be  discussed  later. 

4.3.2  Velocity  Profiles  at  Mach  2 Measured  with  Circular  Pitot  Tube 

Figure  32  displays  velocity  profiles  at  Mach  2 as  deduced  from  the  traverse  measurements 
of  the  circular  pitot  tube  at  the  axial  Stations  of  4 and  18  (see  Figs  15  and  30  for  station  identification). 
The  velocity  profiles  were  computed  assuming,  through  the  boundary  layer,  constant  static  pressure 
equal  to  the  local  wall  condition  in  the  absence  of  the  pitot  tube;  and  constant  stagnation  temperature 
equal  to  the  tunnel  settling  chamber  temperature.  The  boundary-layer  thickness,  6995,  is  that  where 
the  local  velocity  is  99.5  percent  of  the  local  boundary -layer  ‘edge’  velocity,  and  a displacement  correc- 
tion for  the  circular  tube  corresponding  with  that  of  Section  4.1  has  been  incorporated. 

We  see  on  Figure  32,  that  the  uniformity  is  acceptable  at  both  Stations  4 and  18.  The  loca- 
tions C and  X are  1.25-in.  or  six  undisturbed  boundary-layer  thicknesses  from  the  tunnel  centre-line. 

At  Station  18,  the  nominal  two-dimensional  profiles  are  described  better  by  a l/7th  than  a 1 /9th  power 
law,  although  towards  the  edge  of  the  boundary  layer,  the  profiles  indicate  more  fullness  than  either 
power  law  would  predict.  The  high  mainstream  turbulence  level  might  be  a contributing  factor  here, 
producing,  in  conjunction  with  the  rapid  acceleration  of  the  boundary  layer  from  tunnel  throat  to 
Mach  2 conditions  and  associated  non-equilibrium  situation,  an  increase  in  width  of  the  zone  of  inter- 
mittency98. 

We  note  from  the  table  above  the  profiles  on  Figure  32,  that  the  nominal  undisturbed 
boundarv-layer  thickness,  6995 , is  0.2-in.  near  X = 0.0,  and  the  shape  factor  is  about  3.3.  The  latter 
value  is  ten  percent  higher  than  the  predicted  profile  shape  parameter  for  flat  plate  turbulent  boundary 
layers,  by  the  2D  momentum  integral  method  of  Stratford  and  Beavers' 03 . The  Reynolds  number 
based  on  5995  = 0.2-in.  and  the  local  mainstream  conditions  is  R6  s 1.5  X 10s . 

Some  confirmation  of  the  repeatability  of  measurements  of  the  mean  velocity  profiles  is 
shown  on  Figure  33  at  Stations  4 and  18,  before  and  after  a rebuild  of  the  tunnel  settling  chamber 
which  involved  the  replacement  of  baffles  in  the  wide  angle  diffuser  at  the  settling  chamber  entry. 

The  differences  between  profile  measurements  at  Station  4,  in  particular,  are  perhaps  associated  with 
the  recontouring  of  the  settling  chamber  wall  that  has  produced  small  changes  in  the  tunnel  nozzle 


boundary  layers;  rather  than  to  measurement  errors.  We  note  satisfactory  agreement  on  profile  shape 
factor,  H,  from  run  to  run. 

The  pitot  velocity  profiles  on  Figures  32  and  33  are  replotted  in  semi-logarithmic  (law  of  the 
wall)  form,  as  u/urW  * versus  hur  w /uw  on  Figure  34.  Linear  wall  law  regions  are  displayed  for  all  pro- 
files, and  their  slopes  are  compared  with  three  published  laws  of  the  wall  in  the  literature.  We  see  that 
the  slope  of  the  present  data  falls  between  Spalding’s  2D  incompressible  law1 04 , and  the  compressible 
formulations  of  Winter/Gaudet  (2D)105  and  Hall/Dickens  (3D)90.  The  latter  law  was  obtained  from 
profile  measurements  in  a specially  constructed  supersonic  nozzle  at  Mach  1.6  to  2.0,  where  there  were 
three-dimensional  test  boundary  layers;  Hall’s  formulation  appears  somewhat  closer  to  our  measurements 
than  the  other  two. 

43.3  Velocity  Profiles  at  Mach  4 Measured  with  Circular  Pitot  Tube 

Figure  35  shows  pitot  velocity  profiles  corresponding  with  the  Mach  4 empty  tunnel  flow 
conditions  displayed  on  Figure  31,  with  the  floor  vortex  generator  not  installed  (the  left-hand  side 
results);  and  then  with  the  v.g.  location  as  shown  on  the  bottom  of  Figure  31.  Again  constancy  of 
static  pressure  and  stagnation  temperature  through  the  boundary  layer  has  been  assumed  in  reducing 
the  data. 


At  Station  18,  the  convergence  of  the  floor  boundary  layer  towards  the  centre-line,  without 
the  v.g.,  is  well  illustrated  on  the  left-hand  plots  of  Figure  35.  The  transverse  pressure  gradient  is  fa- 
vourable towards  the  tunnel  centre-line  and  that  results  in  a .entre-line  boundary-layer  thickness  about 
twice  that  at  Stations  C and  X,  1.25-in.  from  the  centre-line. 

With  the  floor  v.g.  installed,  we  actually  ‘over-correct’  the  transverse  thickness  distribution 
of  the  test  boundary  layer,  for  now,  the  thickness  conditions  at  C and  X are  greater  than  on  the  centre- 
line, but  the  non-dimensional  velocity  profiles  show  somewhat  better  agreement  in  shape  across  the 
flow  than  without  the  v.g.  The  stagnation  pressure  level  of  about  171  psia  has  also  been  used  in  gather- 
ing cobra  probe  traverse  data  at  Mach  4,  and  this  pressure  level  corresponds  with  the  pitot  data  shown 
on  the  right-hand  graph  of  Figure  35.  The  centre-line  nominal  boundary -layer  thickness  of  5 = 0.2-in. 
(0.18-in.  with  cobra)  is  similar  to  that  at  Mach  2,  but  now  the  Reynolds  number,  R6  has  increased  to 
2.5  X 10s. 

There  would  appear  to  be  an  effect  of  Reynolds  number  on  the  pitot  profiles  measured  with 
the  floor  v.g.  installed.  The  profiles  at  Station  18  on  the  extreme  right  of  Figure  35  were  taken  with 
the  stagnation  pressure  reduced  to  151  psia  and  somewhat  closer  agreement  appears  between  the  nor- 
malized profiles  across  the  wind  tunnel  flow. 

The  Mach  4 pitot  profiles  with  the  floor  v.g.  installed  are  shown  in  semi-logarithmic  form  on 
Figure  36,  with  the  slopes  of  the  log  law  regions  compared  against  Spalding’s  incompressible  law104, 
and  the  compressible  formulation  of  Winter/Gaudet105 . Again,  the  slope  of  our  experimental  results 
falls  somewhere  between  these  two  wall  laws,  although  now,  the  measurements  provide  ‘bumpier’ 
distributions. 

4.3.4  Effect  of  Reynolds  Number  on  Pitot  Velocity  Profiles  at  Mach  2 and  4 

The  effect  on  the  pitot  measurements  at  a given  station,  of  changing  Reynolds  number,  is 
illustrated  on  Figure  37.  An  increase  of  Reynolds  number  per  foot  at  Mach  2,  from  roughly  9 to  14 
million  produces  a slight  increase  in  profile  fullness  and  small  reductions  in  boundary-layer  thickness, 
as  we  see  on  the  left-hand  plot.  At  Mach  4,  we  see  on  the  right-hand  graph  of  Figure  37,  that  the  effects 
of  a small  increase  in  Reynolds  number  from  12  to  14  million  appears  more  impressive  in  the  lower 
regions  of  the  normalized  profiles,  but  only  yields  a small  reduction  on  profile  shape  factor  from  8.0 
to  7.7. 


*ur  w Was  computed  using  Bradshaw’s  Preston  tube  correlation  (Ref.  101). 


4.3.5  Cobra  Probe  and  Pitot  Measurements  in  Empty  Tunnel  at  Mach  2 

Figure  38  provides  us  with  a comparison  between  empty  tunnel  Mach  2 profile  measurements 
taken  with  the  cobra  probe  and  with  the  circular  pitot  tube,  at  axial  stations  two  undisturbed  boundary- 
layer  thicknesses  apart,  near  X = 0.0-in.  On  the  left  of  Figure  38,  vertical  traverse  height  in  inches  of 
movement  is  plotted  versus  local  Mach  number  through  the  boundary  layer,  while  on  the  right-hand 
side  the  profiles  are  normalized  and  plotted  as  h/5,  *,  h/0j , and  h/5995  versus  u/uE . 

Although  the  cobra  and  pitot  profiles  were  not  measured  at  exactly  the  same  locations,  the 
discrepancy  in  Mach  number  profiles  is  perhaps  larger  than  one  would  have  thought  based  on  the  0.5-in. 
distance  apart  of  the  stations.  The  relevant  displacement  corrections  in  Section  4.1  have  been  included. 
The  upstream  (circular  pitot)  Mach  number  profile  demonstrates  less  retardation  in  the  inner  region  of 
the  boundary  layer  than  does  the  more  downstream  cobra  result,  and  this  is  further  emphasized  on  the 
normalized  profiles.  Changing  the  displacement  correction  of  the  rectangular-mouth  cobra  probe  with- 
in reason  near  the  wall,  does  not  collapse  the  cobra  data  on  to  the  circular  tube  data,  and  we  feel  that 
in  all  likelihood,  the  proximity  of  the  support  stem  of  the  cobra  head  to  the  measuring  pitot  orifice  is 
perhaps  the  cause  for  the  discrepancy.  Figure  19  shows  the  design  dimensions  of  the  cobra  head,  where 
the  leading-edge  of  the  support  stem  is  better  than  five  probe  head  heights  aft  of  the  pitot  orifice.  This 
should  be  satisfactory  for  most  measurements  away  from  the  wall,  and  the  comparison  between  cobra 
and  circular  tube  results  is  satisfactory  for  h > 0.050-iri.  It  is  conceivable  that  the  flow  passage  formed 
between  the  probe  head,  the  probe  stem  and  the  test  surface,  as  the  latter  is  approached,  provides  too 
much  internal  contraction,  and  the  shock  wave  about  the  wedge-shaped  support  stem  detaches  and  pops 
forward  to  disturb  and  attenuate  the  local  Mach  number  at  the  cobra  pitot  orifice.  We  note  all  meas- 
urements near  the  wall  yield  supersonic  flow.  Clearly,  the  local  three-dimensional  separation  about  the 
circular  tube  and  the  cobra  head  very  close  to  the  wall  will  show  differences,  which  will  also  affect  the 
results. 

If  we  accept  the  circular  tube  measurements  as  the  more  definitive,  because  of  the  previous 
vast  experience  with  that  frontal  geometry,  then  we  should  keep  in  mind  the  observed  discrepancies 
between  the  circular  and  rectangular  tubes  when  we  review  the  skewed  boundary -layer  profiles  in 
Sections  4.5,  4.6  and  4.7.  On  the  other  hand,  the  circular  pitot  tube  results  demonstrate  some  unusual 
‘kinks’  in  the  normalized  distributions  for  h/6995  < 0.2,  and  we  might  then  prefer  the  smoother  cobra 
probe  results.  However  we  argue  the  merits  of  either  probe  geometry,  the  profile  results  from  the  cobra 
probe  will  be  systematic  amongst  themselves,  and  better  than  satisfactory  results  appear  to  have  been 
obtained  with  this  cobra  probe  in  previous  experiments34-49,  when  compared  with  calculations. 

4.3.6  Integral  Parameters  for  Flow  in  the  Empty  Tunnel  at  Mach  2 and  4 

The  boundary-layer  integral  parameters  and  local  surface  skin  friction  results  as  deduced 
from  the  circular  pitot  and  cobra  probes,  are  plotted  on  Figures  39  to  42,  for  the  two  empty  tunnel 
cases  at  Mach  2 and  4.  The  Mach  2 flow  illustrates  a very  small  increase  in  boundary-layer  thickness, 
displacement  thickness  and  momentum  thickness,  along  more  than  3-in.  of  the  test  region  (equivalent 
to  150  momentum  thicknesses).  The  discrepancy  between  the  circular  tube  and  cobra  probe  profiles 
is  again  demonstrated,  with  slightly  smaller  thicknesses  given  by  the  circular  pitot. 

The  locally  measured  wall  temperature,  Tw , is  plotted  on  the  lower  graph  of  Figure  39,  with 
respect  to  the  stagnation  temperature  in  the  tunnel  settling  chamber.  The  wall  appears  to  be  warmer 
than  the  adiabatic  temperature,  TA  w , at  the  corresponding  test  Mach  number,  although  with  inc  'easing 
run  time,  the  wall  tended  to  cool.  Such  a change  in  temperature  with  run  times  from  20  seconds  up  to 
20  minutes,  is  indicated  by  the  change  in  wall  temperature  ratio  from  close  to  0.97  to  less  than  0.96. 
Most  run  times,  when  continuous  motion  cobra  probe  traverses  were  taken,  were  of  the  order  of  20  to 
30  seconds,  and  under  these  conditions,  we  would  expect  the  most  heat  flux  from  the  wall  to  the  flow. 

Pitot  pressures  at  the  wall,  expressed  as  pressure  coefficients,  Cpp , and  some  corresponding 
local  skin  friction  coefficients  from  Bradshaw’s  correlation1 01  * , are  shown  next  on  Figure  40.  The 


See  Equation  4.2.27. 


surface  pitot  measurements  from  the  cobra  probe  provide  us  with  a local  skin  friction  coefficient  con- 
sistently less  than  that  given  by  the  circular  tube  by  ACF  s 0.00035.  Although  the  measurements  of 
skin  friction  by  Quarmby  and  Das106  with  rectangular  pitot  tubes  in  incompressible  flow  demonstrated 
insignificant  difference  from  the  circular  tube  calibration  of  Patel  (where  the  height  of  the  probe  was 
used  in  correspondence  with  the  outer  diameter  of  the  circular  tube),  we  interpret  the  cobra  probe 
geometry  with  its  chamfered  flattened  side  tubes,  to  require  a calibration  in  a high  Reynolds  number 
pipe  flow.  It  is  worthy  of  note  that  a chamfered  three-tube  surface  yaw  probe,  but  with  three  circular 
tubes  affixed  together,  was  used  in  conjunction  with  Patel’s  calibration  in  Reference  108,  when  the 
surface  direction  of  the  flow  was  initially  unknown. 

We  further  note  in  Figure  40,  that  the  circular  tube  result  from  Bradshaw’s  correlation101  is 
less  by  ACF  as  0.0004  than  the  adiabatic  wall,  Spalding/Chi102  result.  The  discrepancy  is  in  all  likeli- 
hood bound  up  with  the  recovery  of  the  test  boundary  layer  from  its  large  degree  of  three-dimensionality 
in  the  expanding  nozzle  flow  to  one  with  a very  small  residual  convergence  (toward  the  centre-line)  in 
the  working  section.  This  would  cause  a lesser  CF  value  in  our  case  than  the  Spalding  flat  plate  answer 
at  corresponding  values  of  Rg  ^ ^ . The  small  increase  in  wall  temperature  above  that  of  the  adiabatic 
wall  condition  should  not  be  a reason  for  the  discrepancy.  No  wall  temperature  measurements  were 
taken  at  Mach  4.  The  Mach  4 centre-line  values  of  integral  parameters  that  we  see  in  Figure  41  again 
show  only  small  variation  with  distance  along  the  working  section.  The  enhanced  values  of  thickness 
to  the  sides  of  the  centre-line,  at  Stations  C and  X,  illustrate  the  effect  of  the  upstream  vortex  generator 
(Fig.  31)  on  the  floor  boundary  layer,  in  overcorrecting  the  massive  convergence  of  the  empty  tunnel 
flow.  A reduced  height,  plus  a reduction  in  sweep  of  the  v.g.  would  tend  to  flatten  the  transverse  dis- 
tribution of  boundary-layer  thickness. 

The  local  skin  friction  results  in  the  empty  tunnel  Mach  4 flow  are  shown  in  Figure  42.  There 
is  a small  increase  in  CF  as  we  move  downstream,  no  doubt  from  the  mild  ‘attachment  type’  flow  along 
the  tunnel  centre-line  that  is  caused  by  the  floor  vortex  generator.  Again  the  cobra  probe  provides  a 
local  skin  friction  coefficient  less  than  the  circular  tube  result  (both  using  Bradshaw’s  correlation)  of 
ACf  =5  0.0002.  In  this  Mach  4 flow,  the  adiabatic  wall,  Spalding/Chi  result  is  greater  than  Bradshaw’s 
circular  pitot  result  by  ACF  = 0.0005.  Any  discrepancy  here  between  our  nozzle  boundary  layer 
experiment  and  flat  plate  results  is  hardly  surprising  because  of  the  substantial  disturbances  introduced 
to  change  the  transverse  distribution  of  boundary-layer  thickness. 

We  see,  then,  that  with  the  uniform  discrepancy  between  the  CF  results  from  the  pitot  and 
cobra  probes,  we  can  tentatively  propose  that  the  cobra  probe  pitot  measurements  at  the  surface  may 
be  used  in  a circular  Preston  tube  correlation,  provided  the  probe  height  is  input  instead  of  the  diam- 
eter, and  the  offset  in  CF  is  noted,  that  was  measured  in  near  zero  pressure  gradient  flow  conditions 
between  the  cobra  probe  and  circular  pitot. 

4.4  Surface  Pressure  Distributions  and  Flow  Visualization  at  Mach  2 and  4 

The  static  pressure  distributions  along  the  boundary-layer  test  surface  and  a selection  of 
photographs  of  oil  streak  patterns  obtained  with  the  wedge  in  the  wind  tunnel  test  section,  are  shown 
on  Figures  43  to  50  for  the  flow  at  Mach  2;  and  on  Figures  51  to  54  for  the  flow  at  Mach  4.  The  limits 
on  axial  position  of  the  wedge  are  shown  on  Figure  15. 

The  test  surface  pressure  distributions  were  taken  at  closely  spaced  intervals  in  wedge  angle 
particularly  near  flow  deflections  where  incipient  separation  was  predicted  to  occur  from  various  ex- 
periments54-64. On  each  pressure  distribution,  we  shall  observe  that  the  theoretical  shock  pressure  rise 
in  inviscid  flow  corresponding  with  the  mainstream  Mach  number  has  been  plotted,  while  the  origin  of 
the  X-abscissa  of  each  pressure  plot  will  be  the  calculated  shock  wave  position  where  it  intersects  the 
given  line  of  pressures.  We  should  note  that  since  the  nominal  undisturbed  boundary-layer  thickness 
is  0.2-in.,  each  inch  along  the  test  surface  corresponds  with  five  undisturbed  boundary-layer  thicknesses. 
Because  of  the  ease  of  this  transformation,  all  distances  in  the  X-direction  remain  as  absolute  units  on 
the  figures.  The  static  pressure  was  measured  relative  to  the  static  pressure  at  position  one,  upstream 
of  the  interaction  region;  see  Figure  15.  Each  pressure  distribution  was  in  the  direction  of  the  undisturbed 
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stream  parallel  to  the  tunnel  centre-line  at  different  y values  or  lateral  distances  from  the  wedge  apex. 

For  clarity,  only  lines  Z,  (£  and  A pressures  were  usually  plotted.  The  photographs  of  the  tunnel  floor 
were  taken  with  the  optical  axis  of  the  camera  either  parallel  to  the  wedge  deflection  surface  or  just  to 
one  side  of  it.  The  flow  along  the  floor  is  from  left-to-right,  thus  we  frequently  view  a reflection  of  the 
flow  visualization  on  the  test  floor  in  the  smooth-machined  wedge  surface.  The  triangular  section  of 
the  wedge  furthest  from  the  floor  was  usually  out  of  focus,  and  has  consequently  been  outlined  in  ink 
to  assist  the  reader.  The  wedge  deflections,  5W , are  quoted  with  respect  to  the  tunnel  centre-line  direc- 
tion, or  in  other  words,  to  a nominal  stream  direction.  In  Sections  4.5  to  4.7  we  shall  see  small  angles 
of  skew  of  typically  1-degree  detected  by  the  cobra  probe  in  the  nominal  two-dimensional  flow  upstream 
of  the  swept  interaction  region. 

4.4.1  At  Mach  2 

Figure  43  displays  surface  pressure  results  when  the  incident  shock  wave  is  weak.  Here  we 
see  results  for  wedge  angles  of  0,  2 and  4-degrees.  The  surface  pressures  with  the  undeflected  wedge 
compare  reasonably  well  with  those  already  discussed  on  Figure  30  in  the  empty  tunnel  flow.  The 
one  difference  is  that  just  downstream  of  the  wedge  trailing-edge  at  X = 3.0-in.  we  observe  local  pertur- 
bations in  the  pressure  distributions,  which  are  largest  near  the  wedge,  and  which  perhaps  are  due  to 
effects  of  the  deflection  of  the  flow  about  the  wedge  surfaces  opposite  to  the  test  deflection  surface, 
feeding  around  the  trailing-edge. 

Where  the  shock  wave  meets  the  boundary  layer,  and  we  have  attached  flow,  we  see  on 
Figures  43  and  44  (5W  up  to  7.5-degrees)  a steep  rise  in  pressure,  not  dissimilar  to  the  measurements 
of  McCabe54.  The  pressure  rises  smoothly  over  typically  five  boundary-layer  thicknesses  upstream  of 
the  inviscid  shock  wave  position  and  this  length  increases  with  increasing  lateral  distance  from  the  wedge 
apex.  In  comparison  with  the  nominal  2D  interaction  measured  by  Seddon52,  we  see  from  Figure  6 
that  the  extent  of  the  upstream  influence  is  not  unlike  the  3D  case,  but  that  the  region  of  interaction 
downstream  of  the  inviscid  shock  position  appears  to  be  many  times  larger  in  2D  than  in  the  3D  flow. 

Above  5W  = 4-degrees,  the  maximum  pressure  downstream  of  the  shock  decreases  with  in- 
creasing distance  from  the  wedge,  and  this  may  be  due  to  the  upstream  influence  of  the  expansion  fan 
from  the  wedge  trailing-edge,  as  it  approaches  the  shock  with  increasing  lateral  distance  from  the  wedge. 
It  is  interesting  to  note  that  McCabe54  used  both  short  and  long  wedges,  but  still  found  the  attenuation 
of  the  maximum  pressure  with  increasing  distance  from  the  wedge.  He  suggested  that  the  reduction  in 
maximum  pressure  level  was  just  a result  of  the  increasing  extent  of  the  downstream  influence  of  the 
shock,  with  the  asymptote  in  pressure  level  not  reached  before  the  trailing-edge  expansion  causes  a 
massive  fall  in  the  pressure  field. 

As  the  wedge  angle  becomes  larger  (see  Fig.  44  for  6W  up  to  7.5-degrees)  the  continuing  rise 
in  shock  pressure  ratio  turns  the  limiting  streamlines  more  and  more  away  from  the  deflection  angle  of 
the  mainstream  flow.  We  note  the  significant  upstream  influence,  with  turning  of  the  surface  oil  flow 
well  upstream  of  the  calculated  shock  wave  location.  The  upstream  effect  increases  with  distance  from 
the  wedge  over  the  limited  region  investigated  (about  106  from  the  wedge*),  and  Stanbrook53  conjec- 
tures that  this  might  be  explained  by  a form  of  spanwise  influence.  He  proposed  that  the  effects  of 
the  shock  presence  would  be  propagated  along  rearward-facing  local  Mach  cones,  which  in  a large  part 
of  the  boundary-layer  thickness  possess  angles  greater  than  the  shock  wave  angle.  Figure  45  continues 
the  story  up  to  a wedge  deflection  angle  of  9-degrees.  The  spacing  of  the  wall  static  orifices  is  such 
that  the  beginning  of  inflections  in  the  pressure  distributions  is  difficult  to  observe  to  provide  evidence 
of  incipient  separation.  From  the  previous  experiments53-54  with  swept  shock  interactions,  we  expect 
incipient  separation  at  6W  = 7.5  to  8-degrees  (Stanbrook53 ) when  p2 /p j = 1.5;  or  6W  = 8.5-degrees 
(McCabe54). 

From  the  oil  flow  photographs  on  Figure  45,  in  particular,  we  see  no  evidence  of  major  coales- 
cence of  the  limiting  streamlines,  nor  of  the  streamlines  becoming  parallel  with  the  calculated  shock 
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wave  position.  Indeed,  only  as  we  approach  a wedge  angle  of  10-degrees  on  Figure  49  do  the  limiting 
streamlines,  from  the  wedge  apex  outwards,  align  themselves  with  the  calculated  shock  wave  direction 
in  the  inviscid  flow;  and  by  11-degrees,  the  occurrence  of  a prominent  3D  separation  line  is  obvious. 

The  limiting  streamlines  coalesce  into  an  envelope  (originating  at  the  wedge  apex)  that  lies  just  upstream 
of  the  calculated  shock  wave  position.  The  boundary  layer  separates  and  rolls  up  as  a free  shear  layer, 
while  the  indication  is  of  an  attachment  on  the  wedge  face  or  on  the  floor  somewhere  near  the  inter- 
section of  the  wedge  with  the  floor.  We  can  suggest  this,  for  virtually  all  of  the  limiting  streamlines 
downstream  of  the  3D  separation  line  emanate  from  a location  near  the  intersection  of  the  orthogonal 
surfaces. 


From  the  pressure  distributions  at  5W  = 11.5-degrees  on  Figure  49,  we  see  that  the  pressure 
rises  sharply  up  to  the  position  of  the  3D  separation  line.  Thereafter  the  pressure  gradient  falls,  since 
the  rapidly  thickening  3D  separated  flow  does  not  withstand  a large  adverse  pressure  gradient.  The 
gradient  increases  once  more  downstream  of  reattachment. 

Our  incipient  separation  result  at  Mach  2 of  about  10-degrees  wedge  angle  is  greater  than  the 
values  found  previously  and  correlated  by  Korkegi64.  It  is  noteworthy  that  although  McCabe  observed 
separation  to  commence  at  5W  = 8.5-degrees,  his  vorticity  theory54  predicted  separation  at  5W  = 11-degrees 
at  a nominal  Mach  number  of  2.  Moreover  a comparison  of  our  maximum  limiting  streamline  deflections 
downstream  of  the  shock  at  Mach  2,  for  wedge  angles  up  to  almost  12-degrees  demonstrate  reasonable 
agreement  with  McCabe’s  vorticity  model,  as  we  see  on  Figure  10. 

In  these  comparisons  with  previous  results,  one  distinctive  feature  is  evident.  Previous  inves- 
tigations used  wedges  with  deflection  angles  such  that  shocks  were  produced  on  either  side  of  the  wedge 
body.  Our  thin  triangular  section  shock  generator  allows  expansions  to  be  produced  from  the  leading- 
edge  in  the  channel  opposite  to  the  test  flow  for  5W  > 5-degrees.  Hence  the  influence  of  the  flow  about 
the  back -side  of  the  wedge  should  be  negligible  when  we  are  searching  for  incipient  separation  details. 

The  bluffness  of  the  wedge  leading-edge  coupled  with  relatively  large  included  angles  at  the  apex,  in 
previously  reported  work  has  also  often  provided  a protuberance  type  separation  about  the  wedge,  in 
the  junction  with  the  floor  (see  Ref.  54).  The  0.006-in.  leading-edge  radius  on  our  triangular-section 
generator  plus  the  small  included  angle,  caused  only  a very  small  horse-shoe  type  separation  about  the 
wedge  foot,  as  we  can  observe  on  many  of  the  flow  visualization  pictures. 

For  wedge  angles  much  greate  than  12-degrees,  the  internal  contraction  to  the  mainstream  is 
such  that  it  is  difficult  to  isolate  the  test  flow  on  the  floor  from  that  on  the  endwall,  and  no  further 
useful  results  can  be  obtained. 

4. 4. 1.1  Effect  of  Sidewall  Vortex  Generator  at  Large  Wedge  Angles 

No  sidewall  vortex  generator  was  present  in  the  flow  fields  discussed  on  Figures  43,  44,  45 
and  49.  In  contrast,  Figure  48  displays  the  surface  oil  flow  visualization  results  at  5W  = 9-degrees 
obtained  without  and  with  the  sidewall  v.g.  situated  at  XVGLE  = 0.0-in.  We  should  first  of  all  explain, 
however,  that  much  of  the  disturbance  to  the  endwall  boundary  layer  that  we  can  view  on  the  afore- 
mentioned figures  was  due  to  upstream  propagation  and  re-circulatory  flow  in  a small  slot  existing  at 
the  intersection  of  the  floor  with  the  endwall.  When  this  slot  was  filled,  the  improvement  in  the  corner 
flow  was  dramatic.  For  example,  if  we  compare  the  upper  photograph  in  Figure  48,  with  the  correspond- 
ing 6 w = 9-degree  case  on  Figure  45,  we  note  that  the  separation  bubble  in  the  comer  flow  on  the  latter 
figure  has  been  eliminated.  The  addition  of  the  sidewall  vortex  generator  does  not  appear  to  alter  the 
test  surface  flow  to  any  noticeable  extent.  We  remark  again,  however,  that  the  limiting  streamlines 
pass  through  the  calculated  shock  wave  position  at  5W  = 9-degrees  and  no  separation  is  perceived. 

The  effect  of  the  sidewall  vortex  generator  on  the  surface  pressures  closes  to  the  wall,  along 
line  D and  E,  is  to  restrict  the  upstream  influence  of  the  shock  wave  because  of  the  thinner  sidewall 
boundary-layer  flow  — see  Figure  46  for  <5W  = 8-degrees.  At  Sw  = 11.5-degrees,  the  sidewall  v.g.  now 
permits  useful  pressure  distributions  to  be  obtained  as  far  as  line  E as  we  note  on  Figure  50. 
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4.4.1.2  Effect  of  Reynolds  Number  at  5W  = 8-Degrees 


Figure  47  shows  surface  oil  flow  visualization  at  Sw  = 8-degrees  in  a Mach  2 flow,  for  a small 
change  in  Reynolds  number,  R6  , from  1.5  to  2.3  X 105 . In  an  overlay  of  the  surface  shear  stress  direc- 
tions, we  are  not  able  to  detect  any  significant  effect  of  this  limited  Reynolds  number  increase. 


4.4.2  At  Mach  4 


There  is  significant  correlation  between  the  qualitative  flow  features  at  Mach  4 and  those  at 
Mach  2.  The  experimental  arrangement  was  altered  in  the  Mach  4 tests,  however,  to  include  a crude 
Schlieren  system  plus  a transparent  floor,  to  provide  the  location  of  the  mainstream  shock  wave  super- 
imposed upon  the  oil  dot  flow  visualization.  All  of  the  pressure  and  flow  visualization  results  were 
taken  with  the  floor  vortex  generator  situated  as  shown  on  Figure  31. 


Figure  51  illustrates  surface  and  Schlieren  flow  visualization  at  small  wedge  angles.  The  cal- 
culated shock  wave  and  leading  expansion  wave  from  the  wedge  leading  and  trailing  edges,  are  now 
inclined  very  much  closer  to  the  wedge  deflection  surface  than  for  the  corresponding  Mach  2 wedge 
deflections.  The  flow  is  again  from  left-to-right.  The  shock  and  expansion  waves  on  the  Schlieren 
photographs  bear  close  resemblance  with  their  calculated  counterparts.  Moreover,  the  boundary  layer 
on  the  wedge  surface  and  its  wake  downstream  are  clearly  visible,  and  we  observe  the  significant  turning 
of  the  limiting  streamlines  on  the  floor,  far  ahead  of  the  actual  leading  wave  of  the  expansion  fan.  Such 
effect  causes  a drastic  fall  in  the  pressure  rise  through  the  shock  as  we  see  on  Figure  52.  The  upstream 
influence  at  5W  = 6-degrees  stretches  as  far  as  five  to  ten  undisturbed  boundary-layer  thicknesses  ahead 
of  the  calculated  shock  wave  position.  In  proceeding  from  static  pressure  orifice  lines  Z to  A,  there  is 
rapid  attenuation  of  the  peak  pressure  rise  as  distance  increases  from  the  wedge  apex,  while  the  calculated 
shock  pressure  rise  is  never  reached. 


From  the  flow  visualization  records  at  5W  =6  and  8-degrees  on  Figure  52,  we  note  parallelism 


of  the  limiting  streamlines  with  the  calculated  shock  wave  directions  at  both  angles,  but  more  particularly 
at  the  larger  wedge  deflection.  In  addition,  there  is  evidence  of  inflections  on  the  Z and  <£,  pressure 
distributions  at  both  6 and  8-degrees,  so  that  incipient  separation  is  deemed  to  occur  close  to  6-degrees. 


With  an  increase  of  flow  deflection  angle  to  10  and  14-degrees,  that  we  see  on  Figure  53, 
strong  three-dimensional  separations  are  caused  by  the  shock  wave.  The  upper  photograph  on  Figure  53 
demonstrates  limiting  streamline  motion  beneath  the  calculated  shock  wave  location,  towards  a separa- 
tion envelope  far  upstream  of  the  shock.  With  the  wedge  angle  increased  to  14-degrees,  the  leading 
expansion  from  the  wedge  trailing-edge  becomes  close  to  the  calculated  shock  location  and,  on  the 
tunnel  centre-line,  attenuates  the  pressure  rise  sharply  at  six  undisturbed  boundary-layer  thicknesses 
aft  of  the  calculated  shock.  There  is  a divergent  region  of  limiting  streamlines  just  ahead  of  the  leading 
expansion  wave,  where  we  anticipate  flow  reattachment.  The  view  of  the  wedge  surface  at  the  bottom 
of  Figure  53  shows  a downflow  of  limiting  streamlines  on  the  wedge  face  close  to  the  comer,  preceded 
by  an  attachment  region.  The  corresponding  pattern  on  the  tunnel  sidewall  of  the  so-called  nominal 
2D  interaction  is  very  three-dimensional  in  character! 


Finally,  an  extreme  case  of  three-dimensional  flow  separation  is  illustrated  on  Figure  54, 
for  5W  = 16-degrees.  The  attachment  on  the  wedge  face  is  particularly  noticeable  now,  near  its  inter- 
section with  the  tunnel  ceiling  (the  identical  flow  occurs  near  the  test  floor).  The  lower  attachment 
zone  would  appear  to  extend  off  the  wedge  face  on  to  the  test  floor  along  the  line  of  the  trailing-edge 
expansion.  The  sidewall  flow  demonstrates  a cellular  structure  not  unlike  that  found  on  2D  aerofoils 
in  buffet  in  transonic  flow19. 


4.5  Cobra  Probe  Traverses  at  Mach  2 and  5W  = 8-Degrees 


4.5.1  Moving  the  Swept  Interaction  Region  Relative  to  the  Cobra  Probe 

A secondary  objective  to  that  of  controlling  a shock-induced  3D  separation  was  to  obtain  a 


better  physical  knowedge  than  currently  exists,  of  the  mean  flow  details  through  a swept  interaction, 
near  or  at  incipient  separation,  and  in  the  presence  of  a shock  wave  sufficiently  strong  to  cause  a large- 
scale  3D  separation. 
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We  saw  in  Section  4.4,  that  at  Mach  2,  there  was  no  clear  evidence  of  separation  from  the  oil 
flow  and  surface  pressure  distributions  up  to  wedge  deflection  angles  of  about  10-degrees.  The  5^  = 8-degree 
case  was  chosen,  therefore,  as  indicative  of  flow  approaching  the  condition  of  incipient  separation.  (The 
corresponding  wedge  angle  at  Mach  4 was  also  chosen,  as  representative  of  a flow  just  past  incipient 
separation,  while  the  Sw  = 16-degree  case  at  Mach  4 was  selected  to  illustrate  the  details  of  the  rolled- 
up  free  shear  layer.) 

With  the  availability  of  only  four  traverse  stations  for  the  cobra  probe  (see  Fig.  22)  we  elected 
to  move  the  interaction  region  relative  to  the  probe,  by  utilizing  the  6-in.  range  of  available  axial  move- 
ment of  the  wedge  shock  generator.  In  so  doing,  we  would  have  the  option  of  traversing  at  any  axial 
station  through  the  pressure  rise  that  was  of  interest. 

We  have  already  seen  in  the  empty  tunnel  flow  discussed  in  Section  4.3,  that  especially  at 
Mach  2,  there  is  only  a very  small  change  of  mean  flow  boundary-layer  properties  with  distance  down- 
stream in  the  working  section  of  the  wind  tunnel.  But  if  we  cause  the  wedge  to  be  situated  anywhere 
in  the  test  section,  will  the  static  pressure  distributions  along  any  given  line  be  congruent,  when  plotted 
with  respect  to  the  calculated  shock  wave  location? 

We  resolved  this  question  by  first  of  all  looking  at  the  surface  oil  flow  visualization  with  the 
wedge  leading-edge  situated  at  various  stations  between  XWLE  = - 2.5-in.  and  +3.0-in.:  a range  of  280 
undisturbed  momentum  thicknesses  in  length  along  the  test  section. 

Figure  55  illustrates  these  oil  flow  visualization  results,  where  by  superposition,  we  can  detect 
no  significant  changes  in  the  surface  patterns  along  the  test  section.  With  this  encouragement,  we  meas- 
ured the  wall  static  pressure  distributions  for  5W  = 8-degrees,  with  the  wedge  leading-edge  at  XWLE 
equal  to  -2.0,  0.0,  1.0  and  2.0-in.  Figures  56,  57  and  58  show  these  wall  pressures  along  the  tunnel 
centre-line,  and  on  either  side  of  it  along  lines  Z and  A,  that  border  the  cobra  probe  traverse  locations. 

We  note  that  the  measurements  along  each  line  at  each  wedge  position  overplot  satisfactorily,  when 
referenced  to  the  calculated  shock  wave  location  and  we  detect  the  aforementioned  spanwise  attenua- 
tion of  maximum  pressure  level  with  increase  of  distance  from  the  wedge  apex. 

With  these  results  in  hand,  we  felt  it  a logical  procedure  to  investigate  the  swept  interaction 
by  moving  the  wedge  relative  to  the  probe  station. 

4.5.2  Cobra  Profile  Measurements  as  Pitot  Pressure,  Total  Temperature  and  Yaw  Angle 

Figure  59  provides  a map  of  the  cobra  probe  traverse  stations  with  respect  to  the  shock  wave 
position  in  a 6W  = 8-degrees  inviscid  flow  at  Mach  2.  We  remember  from  Figure  22  that  the  traverse 
stations  labelled  T1  and  T2  are  situated  on  the  tunnel  centre-line  while  T3  and  T4  are  at  lateral  dis- 
tances: y = 0.028  and  0.15-in.  respectively  away  from  the  centre-line.  Two  traverses,  downstream  of 
the  shock,  numbered  24806  and  24702,  were  positioned  to  be  on  the  projection  of  one  inviscid  flow 
streamline,  that  is,  along  a line  sloped  at  8-degrees  to  the  centre-line.  Constant  reference  to  Figure  59 
will  assist  the  reader  with  the  ensuing  chapter.  Additional  traverse  information  and  integral  parameters 
are  listed  on  Tables  7.1  to  7.3. 

The  measurements  from  the  cobra  probe  (before  their  manipulation  into  velocity  profiles) 
are  shown  on  Figure  60,  as  corresponding  pitot  pressures,  total  temperatures  and  yaw  angles  (with  re- 
spect to  the  tunnel  centre-line  direction)  plotted  against  vertical  distance  from  the  test  surface,  h.  The 
profiles  exhibited  are  those  taken  along  the  tunnel  centre-line  only,  where,  as  on  the  static  pressure 
plots,  the  calculated  shock  wave  location  is  the  origin  for  the  X-axis.  The  flow  is  from  left-to-right. 
Although  we  follow  the  tunnel  centre-line  in  presenting  this  profile  information,  we  can  infer  from  the 
integral  parameters  that  the  qualitative  features  to  be  discussed  may  equally  well  represent  the  flow 
along  the  path  of  an  inviscid  streamline,  because  of  the  relatively  small  wedge  deflection  angle. 

The  pitot  pressure  profiles  as  related  to  the  settling  chamber  stagnation  pressure  are  at  the 
top  of  Figure  60.  The  local  wall  static  pressure  at  eacn  profile  station  is  written  beneath  the  X-axis 
adjacent  to  the  start  of  the  profile.  Any  local  value  of  pitot  pressure  may  hence  be  determined  using 
the  given  horizontal  scale  and  the  pressure  at  the  wall. 
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The  upstream  profiles  exhibit  the  characteristics  of  the  nominal  two-dimensional  flow  existing 
there.  At  two  undisturbed  boundary-layer  thicknesses  upstream  of  the  calculated  shock  wave  location, 
we  observe  a movement  away  from  the  wall  of  the  sonic  line  and  a rapid  decrease  in  pitot  pressure  signal 
as  the  mainstream  is  approached.  This  pitot  pressure  deficit  would  be  consistent  with  passage  through 
a densely  packed  fan  of  compression  waves,  that  commences  ahead  of  the  location  of  the  calculated 
shock  front.  This  outer  gradient  of  pitot  pressure  disappears  suddenly  at  the  calculated  shock  position, 
denoting  a sudden  increase  in  local  Mach  number  as  h -*■  0.5-in.  With  further  distance  downstream,  the 
sonic  line  slopes  back  towards  the  wall  at  two  undisturbed  boundary-layer  thicknesses  downstream  of 
the  calculated  shock  position.  We  note  that  especially  through  the  interaction,  the  pitot  pressure  adjusts 
very  slowly  to  a constant  level  in  the  region  adjacent  to  the  boundary -layer  edge. 

With  these  non-uniformities  in  view  in  the  non-equilibrium  interaction  region  as  the  main- 
stream is  approached,  how  do  we  determine  the  ‘edge’  of  the  boundary-layer  flow?  A local  edge  condition 
is  desirable  to  provide  normalizing  parameters  for  the  traverses,  in  order  to  remove  any  small  inconsistencies 
that  may  appear  from  moving  the  interaction  region  with  respect  to  the  probe. 

We  selected  the  ‘edge’  as  that  location  away  from  the  wall  where  the  pitot  pressure  first  attains 
its  maximum  value,  and  this  is  shown  by  a small  horizontal  line  on  each  pitot  profile  on  Figure  60.  The 
corresponding  reference  boundary-layer  thickness  where  the  resultant  velocity  is  99.5  percent  of  the 
edge  condition  is  also  shown  as  a solid  dot. 


Some  previous  experiments  in  2D  flow  (see  Z warts9 7 , for  example)  have  estimated  the  boundary  - 
layer  thickness  by  writing  a power  law  velocity  profile  relation 


in  the  form 
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The  outer  60  percent  of  the  boundary-layer  traverse  points  were  fitted  via  a least  squares  routine  to  this 
logarithmic  equation,  to  yield  values  of  the  exponent  n and  boundary-layer  thickness  6.  We  felt  that 
such  a technique  for  the  flow  herein  was  questionable  because  of  the  non-equilibrium  appearance  of 
the  profiles.  Our  visual  approach  of  inspecting  the  pitot  profiles  appears  to  provide  valid  and  consistent 
results,  as  we  see  on  Figures  66,  67  and  68,  where  integral  parameters  and  local  streamwise  skin  friction 
coefficients  are  plotted  for  centre-line  as  well  as  off-centre  traverses. 

The  total  temperature,  referenced  to  the  stagnation  temperature  in  the  tunnel  settling  chamber, 
is  viewed  next  on  the  centre  display  of  Figure  60.  The  total  temperature  rises  rapidly  in  the  laminar  sub- 
layer immediately  adjacent  to  the  wall,  before  assuming  some  constancy  in  the  wall  law  region.  There- 
after it  continues  to  increase  and  to  overshoot  the  mainstream  flow  level  by  typically  one  half  of  one 
percent  of  T0 . The  distance  from  the  wall  of  the  peak  stagnation  temperature  value  increases  with 
progression  downstream  through  the  interaction.  The  asymptotic  approach  to  the  tunnel  mainstream 
flow  value  for  total  temperature  is  again  slow.  We  note  that  the  boundary-layer  edge  as  obtained  from 
the  pitot  profiles  is  usually  close  to  the  position  on  the  total  temperature  profile  where  unit  temperature 
ratio  is  reached  subsequent  to  the  overshoot  (see  Fig.  61). 

The  corresponding  yaw  angles  are  illustrated  on  the  lowest  of  the  three  graphs.  The  skew  at 
the  wall,  in  a direction  that  turns  away  from  the  wedge  deflection  surface,  begins  to  increase  from  its 
nominal  2D  value,  at  least  three  boundary-layer  thicknesses  upstream  of  the  calculated  shock  wave  loca- 
tion. As  we  move  downstream  (but  still  upstream  of  the  shock)  the  extent  of  the  yawed  flow  region 
normal  to  the  wall,  and  the  maximum  yaw  angle,  grow  progressively;  while  at  our  prescribed  boundary-layer 


edge,  the  yaw  angle  has  not  returned  to  the  nominal  mainstream  flow  direction  parallel  to  the  tunnel 
centre-line.  The  yaw  angle  continues  to  reduce  up  to  two  undisturbed  boundary-layer  thicknesses 
away  from  the  surface  (profile  25004,  for  instance).  Immediately  upstream  of  the  shock,  at  profile 
25101,  we  note  the  opposite  behaviour,  see  Figure  61:  above  the  boundary-layer  edge,  the  yaw  angle 
increases  slowly,  this  trend  carrying  over  into  the  recovering  flow  downstream  of  the  calculated  shock 
position. 


We  recollect  that  our  diagnosis  of  the  oil  flow  visualization  records  on  Figure  45  contended 
that  no  separation  existed  at  5W  = 8-degrees.  The  raw  profile  information  that  we  have  just  examined 
would  seem  to  support  this  view,  for  we  shall  see  in  Section  4.7  that  the  stagnation  temperature  profiles 
appear  to  be  sensitive  indicators  of  three-dimensional  separated  flow  regions.  In  contrast,  the  Mach  2, 
6W  = 8-degree  temperature  profiles  on  Figure  60  display  no  sudden  variation  in  characteristics. 

4.5.3  Crossflow  and  Stream  wise  Flow 

Let  us  now  examine  the  cobra  profile  information  from  the  four  traverse  stations.  Where 
profiles  at  nominally  the  same  distances  upstream  or  downstream  of  the  calculated  shock  position 
along  the  lateral  traverse  lines  y = constant  have  been  gathered,  they  will  be  overplotted. 

To  begin.  Figure  62  illustrates  the  yaw  angle  distributions  with  respect  to  the  boundary-layer 
edge  yaw  direction,  and  the  corresponding  total  temperature  plots  at  selectee,  intervals  through  the 
swept  interaction  at  5W  = 8-degrees  at  Mach  2.  The  pressure  ratio*  increases  in  the  direction  of  pro- 
ceeding ‘downstairs’,  to  the  left  for  the  yaw  cases  and  to  the  right  for  the  temperature  measurements. 
The  symbol  code  at  the  top  left-hand  corner  distinguishes  the  traverse  locations,  and  Figure  59  should 
be  sought  as  an  aid  in  identifying  these  traverse  positions  relative  to  the  caclulated  shock  wave.  The 
solid  lines  are  usually  drawn  on  the  graphs  through  the  mean  value  positions  of  the  data  from  one  or 
more  traverses. 

In  the  centre  of  the  page,  we  observe  some  small  residual  skew  of  about  1 -degree  in  the  up- 
stream boundary  layer  at  four  undisturbed  boundary-layer  thicknesses  ahead  of  the  calculated  shock 
wave  location,  that  swings  from  a direction  towards  the  wedge  at  the  test  surface,  to  one  of  away  from 
the  wedge  at  h/6 995  = 0.1.  In  contrast,  a significant  skewing  of  the  flow  very  close  to  the  surface  is 
established  at  one  to  two  boundary-layer  thicknesses  ahead  of  the  calculated  shock  wave  location 
(XT  - XSW1NV  = -0.3).  With  continuing  distance  dowstream,  the  yaw  angle  continues  to  increase 
very  close  to  the  wall,  indicating  some  collinear  flow  at  stations  close  to  the  swept  shock  wave.  Y aw 
data  from  traverses  at  different  (lateral)  y locations  at  corresponding  distances  on  either  side  of  the 
shock  overplot  to  a satisfactory  degree  of  comparison. 

On  the  right-hand  side  of  Figure  62,  the  overshoot  character  of  the  total  temperature  distri- 
butions in  the  outer  region  of  the  boundary  layer  is  very  evident  when  referenced  against  the  vertical 
bar  at  TPC/T0  = 1.  As  we  saw  in  the  empty  tunnel  flow  that  was  discussed  in  Section  4.3,  and  as  we 
see  on  Figure  67  for  the  swept  interaction  flow,  the  wall  temperatures  are  slightly  higher  than  those 
expected  on  an  adiabatic  surface,  implying  some  heat  transfer  from  the  wall  to  the  boundary  layer. 

The  presence  of  the  significant  overshoot  in  the  temperature  distributions  means  that  the  total  enthalpy 
defect  in  the  layer 
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should  be  a small  value,  tending  to  zero  for  a boundary  layer  on  a thermally  insulated  surface. 


Wall  static  pressure. 


I 


-43- 


We  might  also  offer  the  comment  that  if  the  total  temperature  distributions  are  expressed  in 
terms  of  the  quotient 
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they  can  be  compared  against  the  Crocco  relationship1 1 1 put  into  a 3D  formulation, 
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where  uB  is  the  resultant  velocity  in  the  viscous  flow.  We  noted  that  when  plotting  the  ratio  of  tem- 
perature differences  as  ordinate  versus  uB  /uE  as  abscissa,  the  results  were  situated  well  above  the  line 
of  100  percent  correlation  on  the  graph.  The  parabolic  relation  of  Walz1 1 2 , 


provided  even  worse  comparison.  Neither  of  these  frequently  assumed  2D  type  temperature  distribu- 
tions appear  to  exist,  therefore,  in  the  swept  interaction  flowfield  presently  under  discussion. 

Taking  the  measured  total  temperature  through  the  viscous  flow,  the  distribution  of  mean 
velocity  can  also  be  found.  (We  might  interject  at  this  stage,  that  since  at  any  given  Mach  number  in 
the  flow,  the  resultant  velocity,  uB , depends  only  upon  the  half  power  of  the  local  absolute  stagnation 
temperature,  whether  we  use  a constant  stagnation  temperature  through  the  boundary  layer  equal  to 
T0 , or  whether  we  use  that  measured  locally  in  the  flow,  causes  little  change  in  the  calculated  result 
for  uB.)  Knowing  the  flow  conditions  at  our  prescribed  boundary -layer  edge,  the  local  mean  velocity 
can  be  resolved  into  streamwise  and  crossflow  components,  u and  v,  with  u = uE  and  v = 0 at  the 
boundary-layer  edge. 

Representative  boundary-layer  profiles  showing  the  streamwise  and  crossflow  components  of 
mean  velocity  are  exhibited  on  Figure  63,  for  some  of  the  traverses  at  Sw  = 8-degrees  and  M,  = 2.  The 
traverse  height  perpendicular  to  the  wall  and  the  velocities  have  been  normalized  by  5995  and  uE  respec- 
tively. 


The  difficulties  of  finding  adequate  representation  for  the  crossflow  profiles  is  illustrated  on 
the  left-hand  side  of  Figure  63,  where  the  experimental  measurements  at  the  calculated  shock  location 
are  compared  with  the  parabolic  form  suggested  by  Mager90 , which  is 
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tan  0 


The  angle,  0,  on  the  test  surface  is  between  the  local  limiting  streamline  direction  and  the  projection  on 
to  the  surface  of  the  inviscid  streamline  at  the  boundary-layer  edge.  This  profile  representation  clearly 
forces  agreement  immediately  adjacent  to  the  surface*,  but  is  substantially  more  full  throughout  the 
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* There  is  a discrepancy,  however,  between  the  flow  near  the  surface  from  the  cobra  probe  traverse 
and  that  on  the  parabolic  distribution.  This  is  due  to  the  fact  that  the  local  wall  shear  stress  direc- 
tion from  the  oil  flow  visualization  has  been  inserted  for  0,  rather  than  the  cobra  probe  determined 
yaw  angle  closest  to  the  wall. 
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boundary  layer  than  its  experimented  counterpart.  The  maximum  value  of  crossflow  velocity  ratio  and 
flow  angle  tends  almost  to  a constant  value  downstream  of  the  calculated  shock  position. 

In  contrast  with  the  crossflow  profiles  on  Figure  63,  the  corresponding  streamwise  flow  de- 
monstrates an  almost  self-preserving  character  through  the  region  of  strong  interaction,  reminiscent  of 
a typical  2D-type  power  law  profile. 

As  we  have  discussed  in  the  introduction,  the  thesis  of  Lowrie56  dealt  almost  exclusively  with 
representing  chosen  velocity  profiles  in  a swept  interaction  region  by  Johnston’s61 -62,63  triangular  model 
for  crossflow,  that  is  sketched  on  Figure  11(a).  Figure  64  shows  our  Mach  2 swept  interaction  results 
at  5W  = 8-degrees  re-plotted  in  polar  or  hodograph  form.  Although  some  evidence  of  triangularity  exists 
in  the  flow  downstream  of  the  calculated  shock  position,  the  outer  regions  in  the  flow  at  stations  up- 
stream of  the  shock  may  clearly  not  be  approximated  by  a straight  line.  Moreover,  the  frequently 
assumed  collinear  inner  region  would  also  appear  to  be  a doubtful  proposition  here,  for  the  limiting 
streamline  angle  at  the  surface,  as  measured  from  the  oil  flow,  always  appears  larger  than  the  yaw  angle 
at  the  smeared-over  junction  point.  The  equation  established  by  Green,  based  on  Lowrie’s  results  for 
the  junction  point  velocity,  in  terms  of  the  undisturbed  upstream  velocity  in  the  mainstream  and  the 
local  shear  stress  velocity  at  the  wall  (see  Section  1.6),  does  not  provide  a suitable  relationship  to  explain 
the  junction  point  data  of  the  profiles  illustrated  on  Figure  64.  A test  of  our  experimental  results  against 
the  Lowrie/Green  junction  point  criterion  on  Figure  11(b)  indicates  poor  comparison  with  their  line  of 
100  percent  correlation. 

It  remains  to  be  seen  whether  the  family  of  hodograph  crossflow  models  suggested  by 
Shanebrook  and  Hatch1 13  and  discussed  in  terms  of  the  compressible,  3D  turbulent  boundary  layer  in 
Reference  114,  would  yield  any  satisfactory  degree  of  comparison  in  a non-equilibrium  boundary  layer, 
such  as  in  the  swept  interaction.  Comparisons  between  selected  members  of  the  hodograph  family  of 
Shanebrook  and  Hatch1 13  and  the  nozzle  wall  crossflow  profile  data  of  Hall  and  Dickens90  were  good. 
But  to  find  the  analytical  profiles,  the  method  of  least  squares  must  be  used  to  fit  the  experimental 
hodographs  with  the  crossflow  models  from  the  subfamilies1 14 . Such  technique  is  not  too  helpful 
when  attempting  to  calculate  a flowfield  before  experimental  results  are  taken.  For  reasons  such  as 
these,  the  use  of  a three-dimensional  boundary-layer  calculation  method  utilizing  the  turbulent  kinetic 
energy  equation1 15  would  appear  a more  worthwhile  pursuit,  in  which  no  assumptions  about  crossflow 
profile  shapes  have  to  be  made. 

4.5.4  Comparison  of  Experimental  Velocity  Profiles  with  the  3D  Compressible  Turbulent 
Boundary-Layer  Calculation  Method  of  Nash1 15 

Since  we  have  so  far  failed  to  find  a model  of  the  crossflow  that  compares  with  the  experi- 
mental results,  we  discussed  with  John  Nash1 16 , whether  the  existing  formulation  of  his  3D  compressible 
turbulent  boundary-layer  calculation  method  (see  Appendix  1)  could  cope  with  the  swept  interaction 
problem.  The  method  had  been  set  up  in  a form1 15  suitable  for  calculating  the  skewed  boundary  layer 
on  a swept  wing  in  transonic  flow,  where  the  mainstream  at  infinity  is  parallel  to  the  root  chord  of  the 
wing  and  no  flow  is  allowed  across  the  streamwise  wing  tip.  Nash  considers  his  formulation  as  adequate 
for  mainstream  Mach  numbers  up  to  at  least  2. 

Using  the  swept  wing  programme  as  supplied,  but  with  zero  sweep,  we  entered  the  Mach 
numbers,  calculated  from  the  experimental  static  pressure  distributions  as  measured  on  the  tunnel 
floor  at  6W  = 8-degrees  and  M,  = 2,  as  well  as  the  upstream  2D  experimental  momentum  deficit  thick- 
nesses across  the  width  of  the  test  flow.  In  the  method,  conditions  in  the  external  flow  are  considered 
to  be  isentropic  and  irrotational.  The  programme  was  then  run  as  though  the  line  through  the  wedge 
leading-edge  parallel  to  the  tunnel  centre-line  was  the  root  chord  line,  and  the  presence  of  the  deflected 
wedge  boundary  was  thus  ignored.  This  clearly  will  lead  to  under-turning  of  the  calculated  flow,  in  the 
direction  away  from  the  wedge  surface. 

As  an  example,  Figure  65  shows  the  measured  experimental  profiles  in  directions  parallel 
and  normal  to  the  tunnel  centre-line  at  traverse  25104  downstream  of  the  calculated  swept  shock  wave 
location,  compared  with  calculated  profiles  at  lateral  stations  0.5-in.  on  either  side  of  the  traverse  point. 
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Calculated  shear  stress  profiles  are  also  presented.  We  observe  qualitative  agreement  between  calculation 
and  experiment,  but  the  calculated  ‘crossflow’  is  substantially  less  than  that  measured,  while  the  calcu- 
lated ‘streamwise’  flow  evidences  a fuller  profile.  We  would  postulate,  however,  that  a satisfactory 
modelling  of  the  experimental  flow  might  well  result  if  the  bounding  constraint  of  the  deflected  surface 
of  the  wedge  could  properly  be  input  to  the  programme. 


4.5.5  Integral  Parameters,  Local  Skin  Friction  and  Wall  Conditions 

Figure  66  illustrates  the  growth  of  the  boundary  layer  through  the  6W  = 8-degree,  Mach  2 
swept  interaction  region.  All  of  the  data  from  the  profiles  are  plotted,  taken  along  both  the  tunnel 
centre  and  off/centre-line  stations,  at  y = 0.028  and  0.150-in.;  and  where  the  three  lines  of  results  are 
identified  with  different  symbols.  The  shaded  symbols  are  at  those  two  traverse  stations  along  the  pro- 
jection on  the  floor  of  one  inviscid  streamline,  downstream  of  the  calculated  shock  wave  position.  We 
note  that  the  qualitative  interpretation  of  the  results  is  the  same  whether  we  stay  on  one  line  of  traverse 
measurements,  or  whether  we  follow  the  inviscid  streamline  direction. 


All  of  the  plots  on  Figure  66,  whether  they  be  boundary-layer  thickness,  6995 ; streamwise 
and  crossflow  displacement  thickness;  or  streamwise  and  crossflow  momentum  thickness,  indicate  a 
thickening  of  the  viscous  flow  with  increase  of  lateral  distance  from  the  wedge  apex  along  the  calculated 
shock  wave  location  — even  with  the  small  lateral  change  of  five  to  six  undisturbed  boundary-layer 
thicknesses*  that  exists  between  the  borders  of  the  traverse  points  (see  Fig.  59).  As  such,  the  flow 
appears  more  quasi-conical  than  quasi-cylindrical,  as  we  discussed  in  Reference  75. 


In  the  axial  direction,  there  is  a lag  in  the  growth  of  the  streamwise  and  crossflow  thickness 
parameters,  until  about  one  undisturbed  boundary-layer  thickness  downstream  from  the  commencement 
of  the  pressure  rise  (Fig.  67),  where  the  local  static  pressure  gradient  is  at  its  maximum.  During  the 
next  two  undisturbed  boundary-layer  thicknesses  downstream,  the  streamwise  flow  displacement  and 
momentum  thickness  parameters  show  rapid  growth  until  the  calculated  shock  position  is  reached,  and 
the  pressure  gradient  begins  to  reduce  significantly. 

The  shape  factor,  H,  the  ratio  of  the  displacement  to  the  momentum  thickness  parameters 
of  the  streamwise  flow,  (see  Fig.  67),  is  non-responsive  to  the  imposition  of  the  pressure  rise,  but  does 
begin  to  reduce  where  the  pressure  rise  tends  asymptotically  towards  the  calculated  value,  and  where 
the  expansion  fan  from  the  wedge  trailing-edge  is  approached.  This  behaviour  of  H is  consistent  with 
the  almost  self-preserving  character  of  the  streamwise  flow  profiles  that  was  evident  on  Figure  63,  but 
would,  on  the  other  hand,  appear  an  insensitive  indicator  of  what  surely  is  a non-equilibrium  flow.  The 
distribution  of  streamwise  shape  factor  H,  however,  does  indicate  a small  dip  close  to  where  the  maxi- 
mum pressure  gradient  occurs.  This  dip  amplifies  with  increasing  shock  strengths  as  we  see  on  Figure  75 
at  Mach  4. 


Supplementing  the  performance  of  the  streamwise  flow  parameters,  we  see  that  the  corres- 
ponding crossflow  thicknesses  increase  throughout  the  interaction  region,  reducing  only  as  the  wedge 
trailing-edge  expansion  field  is  approached. 

On  Figure  66,  we  show  the  thickness  parameters  along  the  tunnel  centre-line,  calculated  by 
the  Nash  3D  turbulent  boundary  layer  method1 15 , for  this  Mach  2,  6W  = 8-degree  flow  case.  The  cal- 
culated results  offer  a near  invariancy  of  boundary-layer  thickness,  and  only  a small  increase  in  the 
growth  of  the  true  displacement  surface  (see  Appendix  1)  with  distance  downstream  through  the  swept 
interaction,  in  contra-distinction  with  the  rapid  growth  of  the  measured  individual  streamwise  and  cross- 
flow  components.  The  plot  of  displacement  thickness,  illustrates  clearly,  that  the  viewing  of  the  com- 
ponents in  only  the  streamwise  and  crossflow  directions,  does  not  necessarily  provide  the  correct  physical 
impression  of  the  true  displacement  effect.  In  a qualitative  sense,  the  true  displacement  thickness  differs 
from  the  streamwise  flow  value  because  the  boundary  layer  causes  some  flow  at  right  angles  to  the 
direction  of  the  external  streamlines1 18. 


Measured  from  wedge  apex. 
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On  Figure  67,  we  also  show  the  measured  wall  temperatures  plotted  with  regard  to  the  corre- 
sponding adiabatic  wall  condition.  At  the  slightly  elevated  temperature  of  the  test  wall,  the  measured 
temperature  rise  through  the  maximum  pressure  gradient  region,  is  not  as  noticeable  as  on  the  corre- 
sponding calculated  adiabatic  wall  temperature  distribution. 

Most  2D  compressible  calculation  methods  of  turbulent  boundary-layer  flows,  when  compared 
with  experiment,  will  determine  the  momentum  thickness  adequately,  the  displacement  thickness  satis- 
factorily, and  the  skin  friction  (very  often)  to  only  fair  or  poor  accuracy34 . Clearly,  the  crucial  test  for 
the  success  or  failure  of  a computation  method  in  compressible  flow  will  still  be  the  accurate  prediction 
of  the  most  sensitive  and  important  parameter:  the  local  skin  friction  coefficient,  CF . The  measured 
distribution  of  surface  pressure  shown  on  Figure  67  is  not  unlike  that  measured  in  an  unseparated,  two- 
dimensional  flow  — see  Figure  8.  The  distribution  of  local  skin  friction  coefficient,  CF,  that  is  shown 
on  Figure  68,  also  appears  to  possess  the  same  qualitative  form  as  in  2D  flow  — a rapid  fall  in  CF  as  the 
shock  induced  pressure  ratio  increases  to  its  maximum  pressure  gradient  condition,  followed  by  a more 
gentle  reduction  where  the  pressure  gradient  is  virtually  constant,  to  a near  exponential  return  to  equi- 
librium further  downstream. 

Let  us  examine  the  plot  of  local  skin  friction  coefficient  more  closely.  We  first  note,  as 
reference  datums,  the  empty  tunnel  CF  values  for  the  circular  pitot  tube  and  the  cobra  probe  alculated 
according  to  Bradshaw’s  correlation101 , and  the  offset  between  them  of  ACF  — 0.00035  that  we  dis- 
cussed in  Section  4.3.  The  various  wall  pitot  correlations  for  the  cobra  probe  that  were  introauced  in 
Section  4.2  are  plotted  as  the  large  bulk  of  data  points  that  are  referenced  against  the  ‘ZPG,  CFBRAD, 
COBRA  PROBE’  line.  Towards  the  upper  part  of  the  CF  plots,  we  s >e  the  CF  results  calculated  accord- 
ing to  the  correlation  of  Sommer  and  Short120  that  is  based  on  the  streamwise  momentum  thickness 
Reynolds  number  and  the  intermediate  temperature  hypothesis.  This  latter  correlation  produces  data 
points  from  the  cobra  probe  profiles  in  the  upstream  undisturbed  boundary  layer  that  are  in  good  agree- 
ment with  the  undisturbed  values  calculated  according  to  Bradshaw’s  correlation  with  the  circular  pitot 
tube. 

In  the  region  where  the  gradient  of  wall  static  pressure  remains  more-or-less  constant,  the 
crossflow  component  of  wall  shear  stress  increases  steadily  as  we  may  infer  from  the  lower  graph  on 
Figure  68,  where  the  divergence  in  limiting  streamline  direction  is  plotted  with  reference  to  the  direc- 
tion of  the  undisturbed  boundary-layer  flow  upstream  of  the  interaction.  On  this  basis,  we  would 
expect  that  the  magnitude  of  the  resultant  surface  shear  stress  would  fall  less  rapidly  (as  determined  by 
the  cobra  surface  pitot  measurements)  than  the  streamwise  component  alone  (computed  from  the 
cobra  probe  velocity  profiles)  and  displayed  on  the  top  of  the  CF  plot.  That  the  opposite  observation 
is  made  from  the  CF  measurements  on  Figure  68,  is  difficult  to  rationalize.  We  note  moreover,  that 
the  Nash  calculation  method  provides  quite  satisfactory  agreement  with  the  correlation  of  the  velocity 
profile  data  according  to  the  Sommer  and  Short1 20  CF  method. 

i 

Some  further  points  regarding  the  Bradshaw  and  Hopkins/Keener  Preston  tube  correlations 
are  in  order,  because  of  differences  obtained  depending  on  the  manner  in  which  the  Mach  number  is 
calculated  at  the  boundary-layer  outer  edge. 

The  graph  at  the  top  of  Figure  68  shows  three  distributions  of  Mach  number,  calculated 
according  to  (i)  the  Rayleigh  pitot  equation,  using  pitot  measurements  at  the  outer  edge  of  the  boundary 
layer;  (ii)  the  same  equation  with  pitot  measurements  at  roughly  two  undisturbed  boundary-layer  thick- 
nesses away  from  the  wall;  and  (iii)  according  to  the  ratio  of  local  static  pressure  to  tunnel  stagnation 
pressure  in  the  settling  chamber.  In  all  three  computations,  the  local  static  pressure  measured  at  the 
wall  was  assumed  to  be  constant  and  to  exist  along  a normal  to  the  wall  into  the  local  mainstream.  The 
three  computed  Mach  number  distributions  are  plotted  against  the  nominal  Mach  number  in  inviscid 
flow. 


The  Mach  numbers  calculated  according  to  the  local  static  to  stagnation  pressure  ratio,  are 
shown  by  square  symbols  and  fall  progressively  through  the  interaction  region,  consistent  with  the 
static  pressure  rise  illustrated  on  Figure  67.  But  the  pitot  derived  Mach  numbers  only  follow  the  same 
trend  up  to  2.5  undisturbed  boundary-layer  thicknesses  ahead  of  the  calculated  shock  position.  Thereafter 
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those  pitot  measurements  at  the  boundary-layer  edge  (the  circle  symbols)  result  in  a higher  level  of  Mach 
number  throughout  the  interaction  region;  while  the  pitot  measurements  taken  further  out  in  the  flow 
(circle  symbols  with  tails)  yield  a rapid  fall  in  local  Mach  number  up  to  one  undisturbed  boundary-layer 
thickness  ahead  of  the  calculated  shock  location,  followed  by  an  almost  instantaneous  jump  back  to  the 
Mach  numbers  calculated  at  the  outer  edge  of  the  boundary  layer.  In  shock-separated  2D  flows,  the 
work  of  Bogdonoff28  and  Green29 , for  instance,  has  indicated  that  in  the  vicinity  of  the  separation,  the 
outgoing  wave  field  may  be  seen  as  a ‘compression  field  generated  by  a strong  interaction,  but  interrupted 
by  an  expansion  fan  which  is  the  constant  pressure  reflection  of  the  incident  shock  wave’  (see  Ref.  29). 
But  even  before  2D  separation  commenced,  Bogdonoff  and  Kepler28  observed  an  expansion  generated 
at  the  foot  of  the  incident  shock  wave.  That  we  observe  the  same  expansion  phenomenon  in  the  3D 
swept  interaction,  where  there  is  no  3D  separation,  is  of  interest;  and  we  must  presumably  ascribe  the 
differences  in  Mach  number  distribution  from  the  two  sets  of  pitot  measurements  as  those  arising  from 
a non-uniform  static  pressure  normal  to  the  wall,  as  a result  of  the  axial  stretching  of  the  shock  foot 
into  a compression  fan. 

If  we  now  return  to  the  distributions  of  local  skin  friction  coefficient  in  the  centre  of  Figure  68, 
we  note  the  Hopkins/Keener  correlation  based  on  wall  temperature  (adiabatic),  yields  lower  CF  results 
with  the  Mach  number  from  the  pitot  measurement  at  the  outer  edge  of  the  boundary  layer  than  from 
the  Mach  number  calculated  according  to  the  ratio  of  static  pressure  to  undisturbed  mainstream  stagna- 
tion conditions.  However,  the  differences  in  the  CF  results  are  probably  not  significant.  Based  on 
M(P/PO),  the  Bradshaw  correlation  provides  a slightly  increased  value  of  CF  near  the  calculated  shock 
position  than  either  of  the  Hopkins/Keener  results. 

To  conclude,  the  lower  graph  on  Figure  68  shows  a rapid  initial  growth  in  the  divergence  of 
the  local  limiting  streamline  direction  with  respect  to  the  upstream  flow  direction,  followed  by  a more- 
or-less  constant  value  of  PCL  and  a subsequent  region  of  slow  decay.  We  note  that  the  limiting  stream- 
line direction  through  the  interaction  region  as  measured  from  the  oil  flow  visualization,  always  leads 
the  nearest  value  of  yaw  angle  to  the  wall  that  is  obtainable  from  the  cobra  probe  (the  vertical  bars  on 
the  measurements  from  the  oil  flow  indicate  the  maximum  variations  in  @CL  determined  from  seven 
separate  tunnel  runs);  but  that  in  the  relaxing  region  downstream  of  the  calculated  shock  position,  the 
measurements  recover  to  one  of  good  comparison.  We  infer  from  this  discrepancy  in  @CL  through  the 
interaction  region,  that  yaw  conditions  near  the  wall,  between  the  wall  and  the  first  data  point  meas- 
urable with  the  cobra  probe  are  still  changing  rapidly  and  that  regions  of  collinear  flow  are  very  small. 
Such  an  observation  has  already  been  implied  by  the  polar  velocity  plots  on  Figure  64  and  suggests  that 
a significant  law  of  the  wall  region  does  not  exist  in  this  interaction  region. 

As  we  would  expect,  the  omission  of  t..e  deflected  wedge  surface  as  a boundary  constraint 
in  the  3D  boundary-layer  calculation,  causes  the  computed  skew  angle  at  the  wall  to  be  less  than  that 
measured,  by  a value  roughly  equal  to  the  deflection  angle  of  the  wedge. 

In  general,  the  yaw  angle  at  the  outer  edge  of  the  boundary  layer  downstream  of  the  calculated 
shock  position,  is  somewhat  less  than  that  measured  well  away  from  the  surface  (at  two  undisturbed 
boundary-layer  thicknesses),  and  which  in  turn  is  less  than  the  deflection  angle  to  the  flow  provided  by 
the  8-degree  wedge  setting. 

4.6  Cobra  Probe  Traverses  at  Mach  4 and  5W  = 8-Degrees 

4.6.1  Cobra  Profile  Measurements  as  Pitot  Pressure,  Total  Temperature  and  Yaw  Angle 

In  discussing  the  oil  flow  visualization  and  wall  static  p^ssure  distributions  at  Mach  4,  we 
discovered  that  incipient  separation  occurred  at  a wedge  deflection  angle  close  to  6-degrees  (see  Fig.  52). 
Hence,  we  might  expect  to  observe  some  significant  qualitative  differences  between  the  unseparated 
5W  = 8-degree  flow  profiles  at  Mach  2 (Fig.  60)  and  those  at  the  same  wedge  angle  at  Mach  4 (Fig.  70) 
where  large  scale  3D  separation  is  imminent.  The  centre-line  pitot,  stagnation  temperature  and  yaw 
profiles  on  Figure  70  follow  the  presentation  arrangement  established  on  Figure  60  in  Section  4.5,  and 
the  qualitative  aspects  of  the  Mach  4 flow  are  analogous  to  those  interpreted  at  Mach  2.  The  Mach  4 
traverse  locations  are  shown  on  Figure  69,  and  it  will  be  noted  that  no  off /centre-line  measurements 
were  taken. 
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Two  distinctive  features  between  the  Mach  4 and  Mach  2 flows  at  5W  = 8-degrees  are  clear, 
however.  Close  to  the  calculated  shock  wave  location  at  Mach  4 in  Figure  70,  the  rate  of  change  of 
pitot  pressure  with  distance  normal  to  the  surface  is  small,  in  sharp  contrast  with  the  measurements  at 
Mach  2 (on  Fig.  60).  Such  behaviour  is  reminiscent  of  the  pitot  profiles  through  a 2D  shock  separated 
boundary  layer29 . Secondly,  the  rate  of  change  of  total  temperature  with  distance  from  the  wall  in- 
creases as  the  calculated  shock  position  is  approached,  implying  more  vigorous  exchange  of  stagnation 
enthalpy  in  the  flow. 


The  yaw  angle  at  the  prescribed  boundary-layer  edge  is  close  to  zero  degrees  up  to  the  calcu- 
lated shock  position,  and  a tendency  toward  the  wedge  deflection  angle  is  shown  only  farther  out  in  the 
flow. 

4.6.2  Crossflow  and  Streamwise  Flow 


In  comparison  with  the  6W  = 8-degree  flow  at  Mach  2 illustrated  on  Figure  62,  Figure  71  shows 
the  Mach  4 profile  distributions  of  yaw  (with  respect  to  the  resultant  velocity  vector  at  the  outer  edge 
of  the  boundary  layer)  and  total  temperature,  at  the  same  wedge  deflection  angle.  The  most  obvious 
qualitative  feature  of  difference  between  the  two  flows  is  that  the  characteristic  smoothness  of  the  dis- 
tributions at  Mach  2 appears  to  be  changed  to  one  of  ‘hills  and  valleys’  at  Mach  4,  particularly  on  the 
total  temperature  distributions.  Whether  the  ‘wiggles’  on  the  temperature  plots  are  introduced  from 
disturbances  propagating  from  the  floor  vortex  generator  upstream  of  the  test  section,  is  debatable. 

The  yaw  distributions  at  Mach  4 have  also  changed  in  character  in  comparison  with  their  Mach  2 counter- 
parts, the  middle  region  of  the  Mach  4 boundary  layer  holding  a more  constant  yaw  angle  up  to  the  cal- 
culated shock  wave  position. 
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The  streamwise  flow  and  crossflow  velocity  ratio  profiles  are  drawn  on  Figure  72.  The  fea- 
tures of  the  yaw  distributions  seen  on  Figure  71  are  carried  over  to  the  crossflow  plots,  while  the  absolute 
magnitude  of  the  non-dimensional  crossflow  velocity  is  about  the  same  as  at  Mach  2.  The  normalized 
streamwise  flow  velocity  profiles  at  Mach  4 have  also  changed  character  in  the  region  close  to  the  wall, 
tending  toward  a 2D  type  separated  flow  profile.  The  Mach  4 polar  plots  of  velocity  on  Figure  73,  as  at 
Mach  2,  do  not  resemble  a triangular  form  except  at  the  traverse  station  just  upstream  of  the  calculated 
shock  wave  location. 

46.3  Integral  Parameters,  Local  Skin  Friction  and  Wall  Conditions 

The  boundary-layer  thickness  parameters  are  shown  on  Figure  74  for  the  Mach  4,  §w  = 8-degree 
flow  case.  Again  the  locus  of  the  point  of  99.5  percent  of  local  mainstream  velocity,  5g95 , indicates 
little  change  with  distance  downstream.  But  at  roughly  seven  undisturbed  boundary-layer  thicknesses 
ahead  of  the  calculated  shock  wave  position,  there  is  a noticeable  decrease  in  the  growth  of  the  stream- 
wise  displacement  thickness,  5 , much  like  that  measured  in  the  2D  Mach  4 decelerating  flow  of  Peake 
et  al  in  Reference  34.  For  this  decrease  in  5]  * to  occur  in  a decelerating  flow,  the  normalized  density, 
pip E , must  be  increasing  at  a more  than  compensating  rate  than  the  normalized  velocity,  u/uE , is  de- 
creasing, in  the  inner  region  of  the  boundary  layer. 

While  5 , * shows  a small  dip,  the  effect  is  amplified  when  we  transfer  across  to  the  streamwise 
shape  factor,  H,  displayed  on  Figure  75.  The  reduction  in  H occurs  in  the  region  of  maximum  adverse 
pressure  gradient,  but  H recovers  again  towards  the  calculated  shock  wave  location.  No  wall  temperature 
measurements  were  made  at  Mach  4,  but  the  adiabatic  wall  temperature  is  shown  for  reference  as  the 
lower  graph  on  Figure  75. 


The  local  skin  friction  coefficients  derived  from  the  cobra  wall  pitot  measurements  and  the 
cobra  velocity  profiles  are  illustrated  on  the  centre  >aph  of  Figure  76.  The  arrangement  of  the  plotting 
of  the  various  CF  correlations  is  analogous  to  those  measurements  for  5W  = 8-degrees  at  Mach  2 that  we 
saw  on  Figure  68.  There  is  a continuing  fall  in  Cr  towards  the  calculated  shock  wave  location,  but  as  no 
wall  pitot  measurements  were  made  further  downstream,  we  see  no  subsequent  plateau  or  recovery  towards 
the  flat  plate  value.  The  Cr.  correlation  of  Sommer  and  Short1 20 , utilizing  the  Reynolds  number  based 
on  the  streamwise  momentum  thickness  at  each  profile  station  follows  the  result  obtained  in  the  empty 
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tunnel  in  nominally  zero  pressure  gradient  flow.  As  at  Mach  2,  the  relative  differences  between  the 
Sommer/Short  and  surface  pitot  correlations  are  in  the  opposite  direction  to  instinct:  that  the  resultant 
wall  shear  stress  reduces  at  a faster  rate  through  the  interaction  than  the  corresponding  component  of 
streamwise  wall  shear  stress. 


The  sudden  increase  in  local  Mach  number  that  was  determined  at  h = 0.5-in.  above  the  wall, 
as  the  calculated  shock  location  was  approached  in  the  Mach  2,  5W  = 8-degree  flow  case,  is  now  more 
evident  at  Mach  4 (Fig.  76).  The  trends  in  limiting  streamline  angle  at  the  wall  on  Figure  76  follow 
closely  those  explained  on  Figure  68. 

47  Cobra  Probe  Traverses  at  Mach  4 and  5W  = 16-Degrees 

47.1  Cobra  Profile  Measurements  at  Pitot  Pressure,  Total  Temperature  and  Yaw  Angle 

Figure  78  illustrates  the  cobra  probe  pitot,  total  temperature  and  yaw  profiles  along  the 
tunnel  centre-line  through  the  Mach  4,  5W  = 16-degree  flow,  where  as  we  have  already  seen  on  the  oil 
flow  visualization  of  Figure  54,  there  is  extensive  upstream  influence  and  a significant  three-dimensional 
separated  flow  region. 

A map  of  the  seven  stations  relative  to  the  calculated  shock  wave  location  at  which  the  cobra 
traverses  were  taken,  is  given  on  Figure  77.  The  stations  were  chosen  to  provide  profile  measurements 
from  upstream  of  the  three-dimensional  separation  line  to  a position  downstream  of  the  calculated  shock 
wave  location.  The  stations  are  situated  along  the  centre-line  only. 

The  pitot  pressure  profiles  on  the  top  graph  of  Figure  78  are  sufficiently  close  in  spacing 
through  the  interaction  to  allow  a construction  of  the  forward  oblique  leg  of  the  shock  bifurcation  at 
the  foot  of  the  wedge  shock.  The  position  of  the  oblique  shock  leg  is  given  by  the  noticeable  and  sud- 
den deficit  in  pitot  pressure  observable  on  the  plotted  distributions,  and  as  we  would  expect,  the  three- 
dimensional  separation  line  is  positioned  very  close  to  the  intersection  of  this  oblique  shock  wave  and 
the  test  surface.  Downstream  of  the  3D  separation  line,  the  presence  of  shear  layers  with  heat  transfer 
is  indicated  by  developing  troughs  in  the  pitot  distributions  that  are  closely  matched  in  position  by 
similar  troughs  on  the  total  temperature  plots.  Figure  79  shows  the  three  measured  profiles  at  the  cal- 
culated shock  wave  location  re-drawn  in  close  proximity.  The  arbitrariness  in  choosing  6995  is  parti- 
cular important  here  and  should  be  noted  on  Figures  78  a»id  79.  Along  with  the  correspondence  of 
peaks  and  troughs  on  the  pitot  and  total  temperature  distributions,  we  observe  a sudden  slowdown  in 
tne  change  of  yaw  angle  between  these  maxima  and  minima. 

From  these  distributions,  we  may  postulate  that  the  deficit  troughs  correspond  with  the  core 
flow  of  a flattened  vortex  structure  of  the  rolled-up  free  shear  layer  from  the  3D  separation  line,  much 
like  the  flow  measured  by  Cooper  and  Hankey  and  reported  in  Reference  74.  The  axis  of  the  vortex 
would  appear  close  to  the  line  of  the  calculated  shock  wave,  the  origin  of  which  is  the  wedge  apex.  A 
stream  surface  in  the  upstream  boundary  layer  will  pass  over  the  rolled-up  flow  and  attach  to  the  floor 
and  wedge  surface  at  the  divergence  locations  of  the  limiting  streamlines  already  pointed  out  in  Section  4.4. 
Clearly,  the  cobra  probe  can  only  resolve  the  local  angle  of  the  mean  flow  in  planes  parallel  to  the  tunnel 
wall.  If  there  is  significant  pitch  of  the  oncoming  streamlines  to  the  cobra  pitot,  where  the  resultant 
velocity  vector  resolves  into  a larger  component  in  the  pitch  plane  than  in  the  yaw  plane,  the  streamwise 
and  crossflow  profiles,  as  resolved  herein,  will  be  in  error.  The  largest  pitch  angles  relative  to  the  cobra 
probe  will  occur  in  a plane  parallel  to  the  floor  that  cuts  the  vortex  core  position.  In  other  words,  the 
cobra  probe  might  be  providing  low  values  of  pitot  pressure  at  the  troughs  of  the  profiles  illustrated  on 
Figures  78  and  79. 


The  rapid  development  in  skew  at  the  wall  as  the  3D  separation  line  is  approached  is  illustrated 
on  Figure  80,  where  the  yaw  at  the  wall  is  more  than  40-degrees  with  respect  to  the  direction  of  the  re- 
sultant velocity  at  the  assumed  outer  edge  of  the  boundary  layer.  Of  particular  relevance  is  the  overshoot 
in  yaw  along  the  outer  part  of  each  profile  approaching  the  calculated  shock  wave  location.  That  this 
region  of  the  flow  yields  the  maximum  overshoot,  at  the  calculated  shock  wave  location,  provides  ten- 
tative support  to  claim  that  concentrated  roll-up  of  the  free  shear  layer  is  somewhere  close-by.  On 
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Figure  81,  we  have  again  plotted  some  of  the  total  temperature  distributions  measured  in  the  Mach  4, 

5W  = 16-degree  flow,  to  show  the  relatively  small  movement  away  from  the  surface  of  the  troughs  with 
distance  downstream.  Further  out  in  the  viscous  flow,  other  sudden  reductions  or  troughs  in  total 
temperature  are  seen  (see  Fig.  78)  that  become  more  pronounced  with  increasing  distance  down- 
stream. 

4.7.2  Crossflow  and  Streamwise  Flow 

The  features  of  the  crossflow  and  streamwise  flow  velocity  profiles,  Figures  82  and  83,  follow 
closely  on  the  remarks  made  about  the  yaw  and  the  pitot  pressure  distributions.  Analogous  troughs  are 
viewed  on  the  streamwise  flow  profiles  consistent  with  the  troughs  in  pitot  pressure  on  Figure  78,  and 
S-shaped  crossflow  velocity  profiles  emerge  because  of  the  overshoot  in  yaw  angle  discussed  previously. 

At  the  3D  separation  line,  the  streamwise  flow  profile  resembles  in  shape  the  corresponding 
profile  in  2D  flow. 

The  polar  velocity  plots  on  Figure  84  for  the  Mach  4,  6W  = 16-degree  flow  case,  exhibit  sub- 
stantial differences  in  the  separated  flow  to  those  measured  at  the  lower  8-degree  wedge  angle,  near  to 
incipient  separation,  as  shown  on  Figure  73.  As  we  progress  downstream  in  the  former  flow,  from  the 
3D  separation  line  towards  the  calculated  shock  wave  location,  the  shape  of  the  hodograph  in  the  outer 
region  resembles  the  cresting  of  a wave  on  the  sea-shore! 

4.7.3  Integral  Parameters,  Local  Skin  Friction  and  Wall  Conditions 

Since  following  the  3D  separation,  a new  boundary  layer  must  commence  to  grow  from  an 
attachment  line  on  the  wedge  surface  close  to  the  intersection  of  the  wedge  foot  with  the  floor  (see  the 
oil  flow  visualization  on  Fig.  54)  this  new  boundary  layer  can  withstand  a relatively  large  pressure  gradient 
imposed  upon  it  (but  modulated  presumably  by  the  roll-up  of  the  free  shear  layer  above  it).  We  thus 
note  a rapid  increase  in  all  of  the  boundary-layer  thickness  parameters  downstream  of  the  3D  separation 
line  as  shown  on  Figure  85.  The  streamwise  shape  parameter,  on  Figure  86,  shows  a large  decrease  at  the 
3D  separation  line  position  compared  with  that  in  the  upstream  undisturbed  mainstream,  and  this  follows 
the  trend  established  in  the  weaker  shock  pressure  rises  discussed  on  Figures  67  and  75.  However,  the 
distribution  of  H shows  large  up-and-down  variations  as  we  approach  the  calculated  shock  wave  position. 

The  arrangement  of  the  local  skin  friction  information  on  Figure  87  is  consistent  with  that 
already  discussed  on  Figures  68  and  76.  Of  significant  interest,  however,  is  the  fact  that  the  wall  pitot 
(cobra  probe)  does  not  measure  the  lowest  pitot  pressure  at  the  3D  separation  line,  but  at  one  undisturbed 
boundary-layer  thickness  further  downstream.  The  local  resultant  skin  friction  as  calculated  from  the 
Preston  tube  correlations,  exhibits  a rapid  increase  as  we  move  away  from  the  minimum  towards  the 
calculated  shock-position,  in  a region  where  the  limiting  streamline  angle  at  the  wall  has  virtually  con- 
siant  deflection,  see  Figure  87.  In  contrast,  the  streamwise  flow  component  (Sommer/Short)  is  decreasing 
as  the  resultant  CF  is  increasing. 

Outside  of  the  boundary  layer  and  downstream  of  the  3D  separation  line,  we  again  see  the 
trough  in  the  Mach  number  distribution;  but  this  feature  now  seems  to  be  carried  over  to  the  three 
Mach  number  distributions  that  we  see  at  the  top  of  Figure  87. 

Downstream  of  the  3D  separation  line,  the  crossflow  angle  at  the  wall,  0CL,  tends  to  a plateau 
value  of  about  45-degrees. 

4.8  Air  Injection  into  a Swept  Shock/Turbulent  Boundary-Layer 
Interaction  in  a Mach  2 Mainstream  Flow 

4.8.1  Wall  Jet  Flow  in  the  Empty  Tunnel 

In  our  initial  statement  of  objectives,  we  stipulated  that  the  goal  of  the  present  investigation 
was  the  control  by  tangential  air  injection  of  a three-dimensional  separated  flow  region  induced  by  a 
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swept  shock  wave.  As  a follow-on  to  the  successful  two-dimensional  blc  experiment  reported  in  Refer- 
ence 5,  we  planned  to  use  a Mach  3 wall  jet  beneath  a Mach  2 mainstream,  to  provide  a jet  velocity 
greater  than  that  of  the  mainstream,  coupled  with  an  excess  of  jet  momentum  flux  slightly  more  than 
the  momentum  deficit  of  the  tunnel  boundary  layer  upstream  of  the  jet  slot. 

Section  2.6  described  the  jet  blowing  apparatus  consisting  of  a 2D  convergent /divergent  super- 
sonic nozzle  installed  in  the  wind  tunnel  floor  as  part  of  a turntable.  The  wedge  leading-edge  could  be 
situated  at  a distance  up  to  14  original  boundary-layer  thicknesses  downstream  of  the  jet  exit  (with  0j  = 
0-degree)*. 

Before  commenting  on  the  degree  of  control  exercised  by  blowing  into  a swept  interaction 
region,  we  first  of  all  review  the  downstream  penetration  qualities  of  the  Mach  3 jet  in  the  empty  tunnel 
flow.  Figure  88  illustrates  the  surface  pressure  distributions  and  oil  flow  visualization  from  one  such 
configuration,  where  the  jet  flow  (displayed  by  superposing  small  black  arrows)  exits  at  -12-degrees 
relative  to  the  normal  to  the  nominal  mainstream  flow  direction  (large  clear  arrow).  The  negative  sign 
given  to  the  jet  exit  setting  implies  anti-clockwise  rotation  of  the  turntable  from  the  nominal  zero  line 
with  the  jet  exit  plane  perpendicular  to  the  mainstream  flow,  when  we  look  at  the  floor  from  above  and 
the  mainstream  flow  is  from  left-to-right.  The  turntable  periphery  is  made  visible  by  the  dark  material 
inserts  used  in  its  construction.  The  rearward  facing  steps  situated  between  the  extremities  of  the  jet 
nozzle  and  the  tunnel  endwalls  introduce  regions  of  disturbed  flow  bordering  the  jet**. 

The  surface  oil  flow  visualization  on  Figure  88  indicates  negligible  change  in  direction  of  the 
jet  flow  about  its  redirected  centre-line  along  the  test  floor,  whose  limiting  streamline  directions  (not 
magnitudes  of  shear  stress)  are  plotted  on  Figure  89  to  compare  with  a flow  case  at  = 0-degrees.  The 
jet  stagnation  pressure  of  167  psia  yields  a calculated  static  pressure  in  the  Mach  3 jet  flow  about  equal 
to  the  Mach  2 mainstream  static  pressure. 

The  abrupt  change  in  slope  of  4-degrees  between  the  jet  nozzle  wall  and  the  tunnel  floor  in 
the  plane  of  the  jet  exit  (see  Fig.  26)  generates  a shock  wave  and  consequent  increase  in  local  static 
pressure,  which  is  almost  immediately  cancelled  by  the  reflection  of  the  shock  wave  within  the  jet,  as 
an  expansion  from  the  free  jet  boundary.  (Figure  103,  showing  some  Schlieren  results  in  a 2D  flow5 , 
may  assist  the  reader  to  visualize  the  disturbances  near  the  jet  exit.  We  should  note,  however,  that  the 
nomenclature  of  j3  and  Pi/P1N  refer  to  the  deflected  wedge  angle  and  the  ratio  of  stagnation  pressure 
in  the  jet  to  that  in  the  mainstream.)  At  one  inch  or  ten  jet  slot  heights  downstream  of  the  jet  exit, 
the  wall  static  pressure  has  returned  to  the  nominal  mainstream  value  within  the  limits  of  ±10  percent 
of  the  reference  undisturbed  static  pressure  of  the  mainstream  at  station  alZ***  (see  Fig.  24). 

What  of  the  uniformity  of  the  jet  flow  across  the  width  of  the  jet  nozzle  exit?  Because  of  the 
large  number  of  static  pressure  orifices  in  the  wind  tunnel  floor,  and  the  physical  obstruction  of  the 
supply  pipe  to  the  blowing  nozzle  beneath  the  jet  nozzle  insert  (Fig.  13),  it  was  not  possible  to  install 
either  the  cobra  probe  or  the  circular  pitot  probe  traversing  gear.  Instead,  as  well  as  relying  on  local 
static  pressure  distributions  and  accompanying  oil  flow  visualization,  we  chose  to  use  the  obstacle  block 
indirect  skin  friction  measuring  devices  that  Nituch1 25  had  investigated  (following  the  suggestion  of 
Rainbird)  in  a low  speed  tunnel  flow. 

The  optimum  dimensions  that  Nituch  found  for  the  obstacle  block  are  sketched  on  Figure  90, 
where  the  height,  width  and  length  of  the  block  are  typically  3,  4.5  and  9 times  the  diameter  of  the 
static  hole  against  which  the  block  is  to  be  abutted.  A semi-circular  cut-out,  positioned  at  the  centre 


* The  so-called  ‘mixing  length’  along  which  distance  the  blc  jet  and  original  boundary-layer  flows 

merge  together,  may  conveniently  be  defined  as  the  distance  between  the  jet  exit  and  the  calculated 
shock  wave  position  from  the  wedge  (in  the  Mach  2 flow)  along  the  line  of  the  jet.  Hence,  for  any 
given  rotation  of  the  jet  efflux,  and  shock  wave  angle  from  the  wedge,  the  local  mixing  length  will 
vary  along  the  length  of  the  shock  front. 

**  Table  8 lists  the  run  parameters  for  the  blowing  blc  tests. 

***  alZ  is  sometimes  written  as  AlZ. 
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of  the  forward  facing  side  of  the  block,  along  with  a drill  bit  in  the  static  hole,  facilitates  the  accurate 
location  of  the  block  with  respect  to  the  inside  of  the  static  hole.  As  with  the  geometrically  similar 
Preston  tube  devices,  the  disturbance  to  the  local  flow  produced  by  a congruent  family  of  obstacle 
blocks,  provided  they  are  immersed  within  the  law  of  the  wall  region  of  a turbulent  boundary  layer, 
will  be  governed  only  by  wall  variables  and  a characteristic  size  of  the  obstacle  (the  height).  The  blocks 
can  be  calibrated  in  a pipe  flow  and  Vinh126  has  demonstrated  that  under  high  Reynolds  number  con- 
ditions, the  calibration  for  the  geometry  of  block  described  is  similar  to  that  of  Patel’s  calibration  for 
the  Preston  tube. 

For  a qualitative  determination  of  the  variation  of  local  wall  shear  stress  in  a flow,  we  can 
clearly  observe  the  disturbance  pressure  due  to  the  obstacle  block  and  compare  it  with  the  undisturbed 
static  pressure,  pt , at  our  reference  Station  1. 

Figure  90  shows  the  static  and  obstacle  block  pressures  (with  reference  to  p, ) across  the  wind 
tunnel  at  Line  A2*  upstream  of  the  jet  slot  — see  Figure  24,  with  the  Mach  3 wall  jet  on  and  off.  Neither 
the  changing  or  yawing  of  the  blocks,  by  10-degrees,  nor  the  status  of  blowing  affected  the  nominal 
uniformity  of  both  the  block  and  static  pressure  distributions  across  the  tunnel  upstream  of  the  jet  slot. 

The  static  and  obstacle  block  pressures  are  next  shown  across  Line  10  (see  Fig.  24),  at  six  slot 
heights  or  about  three  undisturbed  boundary-layer  thicknesses  downstream  of  the  jet  exit.  Figure  91 
shows  that  the  distribution  of  static  pressure  across  the  Mach  3 jet  flow  is  nominally  uniform,  but  the 
block  pressures  display  a triangular  distribution  reminiscent  of  a preferential  cellular  type  flow,  and 
thus  a very  uneven  shear  stress.  (Note  the  relevant  scales  for  pressure  ratio!)  That  the  tunnel  flow  was 
not  itself  responsible  for  the  unevenness  is  shown  by  the  jet-off  pressure  distribution  measured  with  the 
same  blocks,  which  again  indicates  relatively  good  spanwise  uniformity. 

We  tried  the  same  measurements  across  Line  17  with  the  Mach  3 jet  on,  where  the  location 
was  17.5  slot  heights  or  almost  nine  undisturbed  boundary-layer  thicknesses  downstream  of  the  jet 
exit  — see  Figure  92.  The  same  qualitative  results  were  obtained  but  the  relative  dispositions  of  the 
peaks  and  troughs  in  the  block  pressure  distributions  had  moved  by  0.5-in.  (or  five  slot  heights)  in  the 
lateral  direction!  Again  the  flow  with  the  jet  off  showed  only  a small  variation  in  the  lateral  distribution 
of  block  pressures. 

We  postulated  that  perhaps  the  thermometer  or  the  screen  in  the  blowing  pipe  upstream  of 
the  jet  nozzle  contraction  was  responsible  for  the  anomalies  in  the  spanwise  distribution  of  shear  stress, 
and  these  were  consequently  removed.  Figure  93  shows  the  block  pressures  again  across  Line  10  (Fig.  24) 
which  demonstrate  once  more  the  aforementioned  triangular  shape  to  the  block  pressure  distributions, 
irrespective  of  whether  or  not  the  screen  and  thermometer  were  in  the  settling  chamber**. 

We  thus  concluded  that  the  jet  flow  provided  a high  mean  level  of  local  shear  stress  at  the 
wall,  but  was  modulated  substantially  by  a triangular  waveform  across  the  span  of  the  jet.  As  the  local 
skin  friction  coefficient  is  very  roughly  proportional  to  (pB/p  - 1),  all  other  parameters  being  equal,  we 
can  see  that  the  jet  flow  has  introduced  a wall  shear  stress  about  five  times  that  existing  in  the  empty 
wind  tunnel  (compare  Figs  90  and  91).  However,  the  amplitude  of  modulation  generated  by  the  trian- 
gular waveform  across  the  jet  is  typically  0.7  times  the  level  of  the  wall  shear  stress  in  the  undisturbed 
flow. 

4.8.2  Effect  of  0j  with  5W  = 8-Degrees,  with  XWLE  = 1.5-in. 

Section  4.4  and  4.5  have  included  a comprehensive  discussion  of  the  flow  in  the  swept  inter- 
action region  at  Mach  2,  generated  by  an  8-degree  wedge  deflection,  where  it  was  considered  that  in  spite 
of  substantial  skewing  of  the  boundary  layer,  three-dimensional  separation  had  not  yet  begun.  Hence, 
this  flow  regime  provides  an  initial  datum  against  which  to  check  the  effects  of  adding  the  Mach  3 jet 
flow. 


* A2  means  a2  on  Figure  24. 

**  There  were  no  ‘nicks’  or  scratches  across  the  width  of  the  jet  nozzle,  to  cause  the  flow  imperfections. 
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Figure  94  is  comprised  of  a series  of  illustrations  depicting  the  jet  efflux  issuing  beneath  and 
at  various  angles  to  the  direction  of  the  undisturbed  boundary-layer  flow,  of  between  -20  and  +30-degrees, 
where  the  mixing  distance  between  the  jet  exit  and  the  shock  front  varies  from  about  2 to  20  upstream 
boundary-layer  thicknesses,  depending  on  0 ■ . 

In  the  early  days  of  planning  the  blc  experiment,  it  was  thought  that  directing  the  blowing  at 
say  0}  = 30-degrees  would  perhaps  offer  most  benefit  for  controlling  the  3D  separated  flows,  for  in  so 
doing,  the  direction  of  the  jet  surface  shear  must  be  turned  through  more  than  60-degrees  for  it  to  follow 
the  line  of  a typical  3D  separation  at  Mach  2,  generated  by  a wedge  standing  normal  to  the  floor.  It  was 
suggested  that  the  penetration  ability  of  the  jet  flow  should  be  sufficient  to  erase  the  separation  line. 

In  this  simple-minded  approach,  the  effect  of  the  solid  wedge  obstruction  was  neglected  when, 
of  course,  it  acts  as  a bluff  protuberance  to  produce  a near  normal  shock  to  the  oncoming  jet  flow,  and 
a consequent  massive  3D  separation  that  we  view  on  the  lowest  photograph  on  Figure  94.  We  have  now 
obtained  a flow  with  substantial  separation  with  the  jet  on  as  compared  with  the  swept  interaction  at 
6W  = 8-degrees  with  no  jet,  where  no  separation  existed!  If  we  now  rotate  the  jet  efflux  such  that  it  is 
parallel  with  the  undisturbed  boundary  layer  at  the  jet  slot,  we  still  see  that  there  is  a 3D  separation. 
Incipient  separation  of  the  new  boundary  layer  growing  beneath  the  jet  peak  will  occur  at  a wedge  angle 
corresponding  with  the  local  jet  peak  Mach  number,  provided  the  latter  is  greater  than  the  Mach  number 
of  the  mainstream.  We  remember  from  Sections  4.5  and  4.6  that  the  wedge  angle  for  incipient  separation 
decreases  as  the  local  oncoming  Mach  number  to  the  wedge  increases.  Evidently,  with  the  jet  peak  Mach 
number  close  to  three,  the  wedge  angle  for  incipient  separation  is  less  than  8-degrees.  But  the  enhanced 
surface  shear  stress  provided  by  the  jet  is  inadequate  to  destroy  the  3D  separation  at  the  given  mixing 
distance  of  about  four  undisturbed  boundary-layer  thicknesses,  between  the  jet  exit  and  the  wedge 
leading-edge. 

If  we  follow  the  photographs  in  the  centre  of  Figure  94  to  the  top  of  the  page,  we  see  that  by 
aligning  the  jet  efflux  somewhere  between  the  wedge  deflection  angle  of  8-degrees  and  the  shock  wave 
deflection  angle  of  37-degrees  (at  Mach  2),  the  three-dimensional  separation  line  is  caused  to  disappear. 

The  flow  in  the  junction  of  the  wedge  foot  with  the  floor  when  0j  = - 12-degrees  is  well  displayed  on  the 
photograph  at  the  top  left-hand  comer,  where  no  coalescence  of  the  limiting  streamlines  is  observed. 

In  contrast,  similar  photographs  on  the  right-hand  side  of  the  page  show,  very  clearly,  the  3D  separation 
lines  at  0]  = 20  and  0-degrees.  The  mainstream  and  jet  flow  directions  are  again  indicated  with  the 
superposed  transparent  and  solid  arrows. 

On  the  left-hand  side  of  Figure  94,  we  observe  the  surface  static  pressure  distributions  along 
the  centre-line  and  Line  A for  the  range  of  wedge  angles  used  on  the  flow  visualization.  We  see  that  one 
aspect  of  the  controlling  mechanism  of  the  jet  at  0j  = - 12  and  - 20-degrees  (the  half-shaded  circular  sym- 
bols) is  to  cause  a local  expansion  fan  around  the  wedge  leading-edge.  This  reduces  the  maximum  value 
of  shock  pressure  rise  generated  by  the  wedge  deflection  angle  to  a value  less  than  the  calculated  shock 
pressure  rise.  In  addition,  the  jet  is  now  blowing  in  a direction  along  which,  from  the  wedge  apex,  there 
exists  a favourable  pressure  gradient  (in  the  absence  of  the  jet)  due  to  the  attenuation  of  maximum  pres- 
sure rise  with  distance  from  the  shock  generator.  For  these  reasons,  coupled  with  the  substantial  wall 
shear  stress  introduced  by  the  jet  flow,  the  3D  separation  would  appear  to  be  erased.  The  corresponding 
pressure  distributions  along  each  line  of  static  holes  at  the  respective  6}  angles  shown  on  Figure  94  are 
displayed  on  Figures  95(a)  to  95(e). 

48.3  Effect  of  9]  with  6W  = 11.5-Degrees,  with  XWLE  = 1.5-in. 

As  we  have  seen  on  Figure  94,  the  jet  flow  may  be  orientated  with  respect  to  the  wedge  deflec- 
tion angle  to  alleviate  the  shock  pressure  field.  In  a swept  interaction  with  5W  = 11.5-degrees  where  a 
large  scale  separation  was  generated  (see  Section  4.4),  we  see  on  the  oil  flow  visualization  on  Figure  96 
that  a jet  angle  of  0j  = -20  or  - 25-degrees  would  again  appear  to  remove  the  three-dimensional  separa- 
tion. The  mechanism  of  jet  control  is  once  more  indicated  by  the  reduction  in  shock  pressure  ratios 
displayed  in  the  pressure  distributions  along  Line  A and  (£,  on  the  left-hand  side  of  Figure  96.  It  is  note- 
worthy that  the  rolled-up  shear  layers  from  the  3D  separation  lines  at  0j  = 0 and  - 12-degrees  induce 
substantial  troughs  in  the  pressure  distributions  that  are  consistent  with  the  presence  of  vortex  flows. 

A full  set  of  pressure  information  for  each  turntable  angle  shown  on  Figure  96  is  given  on  Figures  97(a)  to  97(d). 
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4.8.4  Effect  of  Jet  Stagnation  Pressure  in  Empty  Tunnel;  and  with  6W  = 11.5-Degrees 

On  the  left-hand  side  of  Figure  98,  we  illustrate  the  effect  on  the  empty  tunnel  wall  limiting 
streamline  patterns  of  increasing  the  jet  blowing  pressure  progressively  from  0 to  214  psia.  As  the 
blowing  pressure  escalates,  the  local  static  pressure  in  the  jet  at  the  exit  plane  increases  from  a value 
below  (with  P0.  = 114  psia)  to  a value  above  that  in  the  tunnel  mainstream  (where  P = 214  psia.  The 

manifestation  of  this  increase  in  blowing  pressure  is  a progressive  outward  lateral  movement  of  the 
mixing  region  between  the  side  extremities  of  the  jet  and  the  base  flows  of  the  rearward  facing  steps 
bordering  the  jet  exit. 

At  a wedge  angle  of  11.5-degrees,  a doubling  of  the  blowing  pressure  (and  hence  the  jet  excess 
momentum)  from  an  overexpanded  to  an  underexpanded  jet  flow  condition  with  0j  = 0-degrees  and 
XWLE  = 3.0-in.  has  no  appreciable  effect  on  destroying  the  3D  separation,  as  we  see  on  the  right-hand 
side  of  Figure  98.  Pressures  along  the  test  wall  corresponding  with  the  flow  visualization  photograph 
at  the  design  blowing  pressure  conditions  of  P0.  = 171  psia,  are  shown  on  Figure  99,  where  evidence 

of  a pressure  trough  along  the  centre-line  downstream  of  the  calculated  shock  position  indicates  the 
presence  of  a rolled-up  vortical  flow  and  no  meaningful  control  of  the  3D  separation.  Clearly,  to  under- 
stand the  character  of  the  new  wall  boundary  layer  beneath  the  jet  peak,  and  the  mixing  region  above 
the  jet  peak  in  these  highly  skewed  flows,  we  require  boundary-layer  profile  measurements  to  supple- 
ment the  static  pressure  distributions  and  oil  flow  visualization  results. 

4.8.5  The  2D  Shock  Wave/Turbulent  Boundary-Layer  Interaction  and  Its  Control  by  Blowing 

As  a postscript  to  blowing  into  the  3D  swept  interaction,  we  thought  it  meaningful  to  repeat 
as  near  as  possible,  one  of  the  2D  flow  cases  reported  in  Reference  5.  To  the  glass  window  in  the  ceiling 
of  the  empty  wind  tunnel  (see  Fig.  13),  we  affixed  an  8-degree  wedge  as  depicted  in  the  top  left-hand 
drawing  on  Figure  100,  and  which  causes  a gross  2D  separation  in  the  absence  of  blowing5,29 . The  wedge 
was  positioned  in  such  a way  that  its  oblique  shock  wave  would  strike  the  test  flow  on  the  tunnel  floor 
somewhere  between  six  to  ten  original  boundary-layer  thicknesses  downstream  of  the  jet  exit,  when  set 
at  0 j = 0-degrees,  and  in  line  with  the  recommendations  for  jet/boundary-layer  mixing  length  reported 
in  Reference  5.  Two  injection  Mach  numbers  were  tried:  Mj  2 2.5,  a value  close  to  the  2.4  value  used 
in  Reference  5;  and  Mj  2 3. 

The  oil  flow  visualization  patterns  for  the  8-degree  wedge  angle  in  the  Mach  2 mainstream, 
at  each  jet  Mach  number,  are  shown  on  the  right  of  Figure  100.  The  mainstream  and  jet  flows  are 
parallel  to  the  tunnel  centre-line  in  a direction  from  left-to-right. 

At  each  jet  Mach  number,  we  observe  a nominal  2D  separation  line  at  the  ‘collision’  location 
of  the  incident  shock  with  the  turbulent  boundary  layer,  with  almost  immediate  reattachment  of  the 
flow  thereafter.  We  note  on  the  left-hand  pressure  distributions  of  Figure  100  the  static  pressure  ratio 
through  the  combined  incident  and  reflected  shocks  is  approximately  2.4  and  decreases  quickly  to  a 
level  nearer  to  the  calculated  pressure  rise  of  2.2,  just  downstream  of  the  separation.  The  jet  would 
appear  to  be  of  benefit  in  providing  a substantial  wall  shear  stress  downstream  of  the  interaction,  in 
accordance  with  the  findings  of  Reference  5. 

The  experiment  was  then  repeated  with  a wedge  angle  of  11.5-degrees,  a jet  Mach  number  of 
three,  and  a mixing  region  between  the  jet  and  original  boundary  layer  of  about  six  original  boundary- 
layer  thicknesses.  A 2D  separation  was  again  obtained,  as  we  see  on  Figure  101,  with  a more  extensive 
separation  bubble  length  than  at  the  8-degree  wedge  angle.  The  experimental  shock  pressure  ratio  is 
now  3.6. 

Downstream  of  the  separation  bubble,  we  again  view  a substantial  length  to  the  oil  streaks, 
implying  a satisfactory  level  of  wall  shear  stress.  Clearly,  the  expansion  fan  from  the  wedge  trailing- 
edge  will  in  part  be  responsible  for  the  re-establishment  of  the  flow  downstream  of  the  separation 
bubble,  as  will  the  jet  flow. 
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A series  of  profiles  through  the  jet/boundary -layer  combination  is  shown  on  Figure  102  and 
taken  from  Reference  5 for  a shock  pressure  rise  just  slightly  less  than  three,  but  in  most  other  respects 
analogous  to  the  flow  on  Figure  101.  The  velocity  profile  that  was  taken  nearest  to  the  peak  of  the 
pressure  rise  demonstrates  retardation  to  the  wall  boundary  layer  beneath  the  jet  peak  and  to  the  mixing 
region  between  the  jet  peak  and  the  tunnel  mainstream,  when  compared  with  the  velocity  profile  at  the 
commencement  of  the  pressure  rise.  Nevertheless,  the  corresponding  Schlieren  flow  visualization  photo- 
graph from  Reference  5,  that  is  shown  as  Test  N on  the  right-hand  side  of  Figure  103,  would  appear  to 
illustrate  that  control  of  the  test  flow  had  been  maintained  through  the  pressure  rise.  The  shock/expansion 
train  in  proximity  to  the  jet  exit  will  introduce  a considerable  increase  in  turbulence  into  the  flow  field, 
and  we  can  detect  no  large-scale  regions  of  separated  flow  beneath  the  Mach  reflection  type  of  shock 
formation.  The  jet  has  apparently  provided  sufficient  momentum  flux  to  ‘carry’  the  original  boundary 
layer  through  the  shock  pressure  rise  of  almost  three,  although  further  downstream,  the  expansion  from 
the  trailing-edge  of  the  wedge  shock  generator  in  the  ceiling  of  the  wind  tunnel  will  also  support  the 
flow  re-establishment. 


5.0  CLOSING  DISCUSSION 

5.1  About  the  Swept  Interaction 

What  then  have  we  found  from  this  investigation  into  swept  interaction  phenomena  in  the 
absence  and  in  the  presence  of  the  boundary-layer  control,  that  is  a new  contribution  to  the  under- 
standing of  such  flows? 

In  the  swept  interactions  investigated  at  Mach  2 and  Mach  4,  we  have  measured  for  the  first 
time,  to  the  author’s  knowledge,  the  mean  flow  and  total  temperature  distributions  in  a compressible 
three-dimensional  turbulent  boundary  layer  prior  to,  and  just  past  incipient  separation,  as  well  as  pro- 
files where  a large-scale  3D  separation  existed.  The  temperature  profiles  appeared  to  be  particularly 
sensitive  at  indicating  the  presence  of  rolled-up  free  shear  layers  from  3D  separation  lines. 

We  have  also  managed  to  document  the  near  incipient  separation  flow  (§w  = 8-degrees  at 
Mach  2)  with  closely  spaced  traverses  through  the  interaction  region,  using  a servocontrolled  self- nulling 
three-tube  yawmeter  that  provided  pitot  pressure,  yaw  angle  and  total  temperature  in  a continuous 
motion  traverse  from  the  test  wall  to  the  local  mainstream.  The  large  number  of  traverse  stations  was 
facilitated  by  the  design  of  the  wedge  shock  generator  employed,  standing  normal  to  the  test  wall,  that 
possessed  an  axial  movement  in  the  test  flow  of  30  undisturbed  boundary-layer  thicknesses  (6  = 0.2-in.). 

To  a lesser  degree  of  completeness,  we  also  sampled  the  post-incipient  separation  flow  and  a strong 
interaction  case,  both  at  Mach  4 (5W  =8  and  16-degrees)  where  in  the  latter,  the  existence  of  a flattened 
rolled-up  free  shear  layer  within  the  confines  of  the  boundary-layer  depth  was  postulated  from  the  pro- 
file measurements,  and  the  corresponding  surface  oil  flow  visualization.  We  did  not  observe  any  closed 
separation  bubbles  that  have  been  suggested  for  infinite  swept  wing  flows30.  The  rolled-up  vortical 
flow  from  the  3D  separation  line,  as  discussed  by  Cooper  and  Hankey74  would  appear  applicable  to  our 
separated  flow.  The  interaction  between  the  viscous  and  the  external  flow  field  is  noticeable  in  all  the 
traverse  measurements  through  the  swept  interaction  region  by  the  very  slow  return  towards  the  cal- 
culated external  flow  field  outside  of  our  prescribed  boundary-layer  edge. 

By  comparison  with  the  extensive  number  of  flow  profiles  obtained  through  the  interaction 
region,  neither  the  Johnston  triangular  model  for  crossflow61  *6 2 , that  Lowrie  adopted56,  nor  the  Mager 
parabolic  representation49*90,  appeared  to  fit  our  measurements.  On  the  other  hand,  our  limiting  stream- 
line deflections  at  the  wall  compared  very  closely  with  the  simple  vorticity  model  suggested  by  McCabe54. 
Rather  than  pursue  the  fruitless  task  of  trying  to  find  a suitable  model  for  the  velocity  profiles  in  the 
swept  interaction  region,  we  felt  that  a relevant  analytical  tool  resided  in  the  use  of  a turbulent  kinetic 
energy  boundary-layer  calculation  method,  where  the  history  of  the  boundary  layer  is  properly  taken 
account  of,  and  the  appropriate  boundary  conditions  can  be  inserted.  In  fact,  the  only  3D  compressible 
turbulent  boundary  layer  method  readily  available  (to  the  author)  that  is  based  on  turbulent  kinetic 
energy  principles  is  that  due  to  Nash1 15*1 2 1 , wherein  no  assumptions  regarding  families  of  velocity  pro- 
file shapes,  nor  the  magnitude  of  the  crossflows,  require  to  be  inserted.  The  results  from  Nash’s  method1 1 5 
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when  compared  with  the  Mach  2,  5W  = 8-degree  experimental  measurements,  appeared  promising.  The 
programme  was  formulated  to  treat  swept  wing  flows  and  consequently  requires  a suitable  modification 
to  the  side  domain  boundary  condition  to  allow  for  the  presence  of  the  deflected  surface  of  the  wedge. 

A two-stage  calculation,  where  the  initial  boundary  layer  is  first-of-all  computed  up  to  the  line  of  the 
rearward  facing  Mach  cone  from  the  wedge  leading-edge  utilizing  a rectangular  co-ordinate  system, 
followed  by  a second  calculation  based  on  a conical  co-ordinate  system  whose  origin  is  the  wedge  apex, 
should,  in  principle,  improve  the  prediction  of  the  flow.  The  solution  of  the  first-stage  calculation 
would  be  used  to  start  the  second  stage. 

The  empty  wind  tunnel  provided  a uniform  test  boundary  layer  at  Mach  2 on  the  nozzle 
sidewall.  A single  sweptback  fence-type  vortex  generator  was  used  to  control  the  thickness  of  the  end- 
wall  boundary  layer  opposite  to  the  deflected  surface  of  the  wedge,  and  to  reduce  the  upstream  influence 
of  the  interaction  of  the  test  shock  wave  with  that  boundary  layer.  At  Mach  4,  on  the  other  hand,  the 
wind  tunnel  nozzle  sidewall  boundary  layer  was  highly  non-uniform  and  a V-type  fence  vortex  generator 
was  situated  upstream  in  the  nozzle  on  the  test  floor  to  redistribute  the  converged  boundary  layer.  The 
redistribution  provided  only  a qualitatively  acceptable  flow  in  terms  of  uniformity  of  boundary-layer 
thickness  across  the  tunnel  width. 

As  McCabe  and  other  workers  have  shown54,  the  calculated  shock  wave  position  provides  a 
satisfactory  origin  against  which  to  compare  the  test  wall  static  pressure  distributions,  in  our  case, 
effectively  plotted  on  a base  of  x/5.  In  the  lateral  extent  of  flow  region  available  to  us,  of  typically 
12  undisturbed  boundary-layer  thicknesses  from  the  wedge  apex,  the  flow  appeared  to  develop  on  a 
quasi-conical  basis.  Incipient  separation  at  Mach  2,  in  particular,  was  not  diagnosed  until  Sw  — 10-degrees, 
a wedge  deflection  angle  much  higher  than  previous  experiments  had  found54 . 

In  the  Mach  4 tests,  we  installed  a transparent  floor  in  the  wind  tunnel  and  we  were  able  to 
obtain  Schlieren  records  of  the  wedge  generated  shock  wave  superimposed  on  the  oil  dot  flow  visuali- 
zation record.  The  experimental  shock  angle  was  in  good  agreement  with  that  calculated  from  oblique 
shock  tables. 

Vas1 27  has  queried  whether  the  effects  of  bluntness  of  the  wedge  leading-edge  could  propa- 
gate laterally  and  downstream,  to  modulate  the  flow  field  generated  from  a sharp  wedge  leading-edge. 

The  blunt  leading-edge,  of  course,  acts  as  a bluff  protuberance  to  the  oncoming  boundary 
layer,  producing  a horse-shoe  shaped  local  3D  separation  about  the  foot  of  the  wedge.  Westkaemper1 28 
has  concluded  that,  in  supersonic  flow  about  circular  section  protuberances  and  for  cylinder  length 
to  diameter  ratios  >1,  the  distance  of  the  singular  separation  point  ahead  of  the  cylinder  leading-edge 
(on  the  axis  of  symmetry)  is  approximately  2.7  times  the  diameter  of  the  cylinder.  That  is,  this  separa- 
tion distance  is  virtually  independent  of  cylinder  length  to  diameter  ratio,  the  ratio  of  the  cylinder 
diameter  to  the  undisturbed  momentum  deficit  thickness  of  the  oncoming  boundary  layer,  and  of 
mainstream  Mach  number  and  Reynolds  number  (based  on  cylinder  diameter).  This  independence 
would  appear  due  to  the  ‘swamping’  influence  of  the  very  strong  pressure  field.  In  our  experiment, 
however,  the  radius  of  the  wedge  leading-edge  is  0.006-in.,  and  the  (undisturbed)  boundary-layer 
thickness  is  0.2-in.  The  ratio  of  ‘effective’  cylinder  diameter,  D,  to  boundary-layer  thickness  is  only 
0.06  in  our  case,  and  we  would  expect  that  the  3D  separation  would  be  relatively  much  closer  to  the 
leading-edge  than  Westkaemper’s  result  due  to  the  very  small  D/5.  (Clearly,  in  the  limit  of  vanishingly 
small  D/5,  the  separation  should  disappear  altogether.)  As  far  as  we  can  determine  from  our  oil  dot 
flow  visualization,  the  effect  of  this  localized  3D  separation  is,  in  fact,  confined  to  the  immediate 
neighbourhood  of  the  wedge  meeting  the  wall,  and  is  masked  by  the  larger  scale  interaction  produced 
by  the  wedge  shock  wave.  Nevertheless,  one  does  sight  convergence  of  the  limiting  flow  very  close  to 
the  wedge  apex  adjacent  to  the  test  deflection  surface  before  the  angle  for  incipient  separation  of  the 
test  boundary  layer  is  reached;  and  this  could  be  caused  by  a lateral  movement  of  the  protuberance 
3D  separation  line  towards  the  wedge  shock  wave. 

Of  greater  significance  to  the  stand-off  distance  of  this  protuberance  type  separation  line, 
we  feel,  is  the  deflection  angle  of  the  wedge  shock  generator  on  the  side  opposite  to  the  test  deflection 
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surface.  Spurious  effects  from  strong  3D  interactions  here  may  well  feed  around  the  wedge  leading- 
edge,  to  cause  some  change  in  the  flow  development,  and  to  perhaps  incite  earlier  incipient  separation 
of  the  test  flow  than  would  be  obtained  otherwise. 

5.2  About  Blowing 

Before  attempting  the  blowing  experiment,  it  was  thought  that  the  most  beneficial  direction 
of  the  blc  air  was  at  a sweep  angle  set  such  that  the  jet  efflux  must  turn  through  more  than  90-degrees 
to  become  parallel  with  the  direction  of  the  3D  separation  (in  the  absence  of  blowing).  However,  as  we 
have  seen,  the  jet  is  directed  against  the  deflected  surface  of  the  wedge,  which  acts  as  a bluff  protuber- 
ance, and  causes  a substantial  3D  separation. 

The  present  results  illustrate  that  in  a junction-type  problem,  such  as  we  have  investigated, 
tangential  air  injection  introduced  with  a jet  excess  momentum  slightly  more  than  the  momentum  de- 
ficit of  the  undisturbed  boundary  layer,  and  with  a mixing  length  of  about  two  to  fifteen  undisturbed 
boundaiy-layer  thicknesses  between  the  jet  and  the  shock  wavet , should  be  orientated  to  blow  in  a 
direction  somewhere  between  the  deflected  wedge  surface  and  the  angle  of  the  shock  wave  (i.e.  blowing 
in  a direction  roughly  along  the  shock).  In  this  manner,  the  mode  of  operation  of  the  jet  is  to  substan- 
tially reduce  the  pressure  rise  of  the  wedge  shock  that  is  experienced  on  the  surface,  by  the  amount  of 
turning  of  the  jet  allowed  around  the  wedge  leading-edge,  coupled  with  the  significantly  enhanced  wall 
shear  stress. 

We  have  measured  no  flow  profiles  in  the  3D  jet/boundary-layer  combination,  and  conse- 
quently we  can  only  rely  on  our  2D  results5  as  a qualitative  guide.  In  the  few  limited  experiments 
with  blowing  that  we  have  managed,  it  is  not  clear,  for  example,  whether  the  advised  direction  of 
blowing  herein  would  carry  over  to  controlling  a 3D  shock-induced  separation  on  a swept  wing. 

More  work  clearly  needs  to  be  done  on  a flow  model  producing  a 3D  separation  without  the 
constraint  of  a solid  vertical  wall  in  proximity  to  the  jet. 


6.0  CONCLUSIONS 

Based  on  measurements  of  the  near-zero  heat  transfer  turbulent  boundary  layer  on  a flat 
wall  in  a swept  interaction  at  Mach  2 and  4,  generated  by  a wedge  standing  normal  to  the  test  wall, 
and  at  a Reynolds  number  of  ^2  X 105  based  on  the  undisturbed  boundary-layer  thickness  of  0.2-in., 
we  conclude  that: 

(1)  The  development  of  the  viscous  flow  in  the  swept  interaction  region  leading  eventually  to  3D 
separation  with  increasing  deflection  angle  of  the  wedge,  is  a gradual,  progressive,  relatively 
steady,  and  essentially  quasi-conical  process  near  the  wedge,  in  which  the  flow  leaves  the  3D 

1 separation  line  as  a free  shear  layer  to  roll-up  into  a flattened  vortex  roughly  within  the  depth 

of  the  original  undisturbed  boundary  layer.  There  appears  to  be  no  sudden  eruption  of 
vortical  fluid  from  the  test  surface.  In  fact,  the  precise  wedge  angle  at  which  incipient  sepa- 
ration occurs  is  particularly  elusive,  unless  the  observation  of  parallelism  of  the  limiting 
streamlines  with  the  calculated  (or  visualized)  shock  wave  is  adopted  as  the  incipient  separa- 
tion criterion. 

(2)  The  optimum  direction  of  a Mach  3 supersonic  wall  jet  to  control  or  destroy  the  shock- 
induced  separation  (produced  by  the  aforementioned  wedge  in  a flow  at  Mach  2)  is  along 
a line  between  the  wedge  deflection  angle  and  its  shock  wave  angle. 


t From  the  wedge  leading-edge,  in  a direction  along  the  shock  wave. 
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2.05 

2.44 

5.13 

11.8 


TABLE  1 

TYPICAL  RUN  CONDITIONS 


PQ  atm.  | psia 


30.15 

36.8 

75.3 

174.0 


Rgll 

4 

x 106 

4 

X 

O 

ON 

5 

X 106 

7 

X 106 

Nominal  (min.) 
Starting 
Pressure 
(atmos. ) 


1.16 

1.58 

3.9 

10.0 


Limited  by  strength  of  widge 


TABLE  2 

INJECTION  NOZZLE  PROFILE  DIMENSIONS  (SEE  FIG.  26) 


Jet  Mach  No. 


Throat  Ht. 
hfc  in. 

Dimension 
'A'  in. 

Dimension 
'B'  in. 

0.0850 

. 3454 

.0073 

0.0695 

.4567 

.0151 

0.0593 

.5300 

.0202 

0.0379 

.6829 

.0309 

0.0236 

.7855 

.0380 

0.0093 

.8879 

.0452 

0.0060 

.9105 

.0469 

1278  5160;  35.85  .00568  0.2301  ; .001308  98.067;  1684  .01838  0.0  0.0  0 

0585  11,280;  78.4  .01242  0.1317  ; .001637  109.64  ; 1882  .02565  198  5.08  1.52' 

02722  24,220;  168.5  .0267  .07623;  .002038  117.94  ; 2023  .0344  339  11.67  3.50: 


TABLE  3 (Cont’d) 


A 


- - — ■ 
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TABLE  4 


BOUNDARY  LAYER  MOMENTUM  DEFICIT  A 


M 


?pu2 
po  0 


atmos 

P, 


psf 


lbwt/f t2 
(spu2} 


lbwt/f t2 
{pu2} 


in. 


611 


ft. 


lbwt/f t 
span 

{ P u2 } ^11 


lb.wt. 
{pu2}  eus 


1.5 


0.4290 


2.05 


4335 


1862 


3724 


.01519 

.001267 


4.72 


1.416 


2.0 


0.3579 


2.87 


6065 


2.44 


5160 


2173 

1848 


4346 

3696 


.01690 

.001408 

.01730 

.001442 


6.12 

5.33 


1.836 

1.599 


3.0 


0.1715 


5.94 

12,580 

5.13 

10,830 


2157 

1860 


4314 

3720 


.01750 

.001458 

.01796 

.001497 


6.28 

5.57 


4.0 


0.07376 


9.72 

20,550 

11.8 

24,960 


1517 

1842 


3034 

3684 


.01727 

.001440 

.01675 

.001396 


4.37 

5.15 


1.884 

1.671 


1.311 

1.545 


20-degrees  to  floor;  HVG  = 0.1-in.;  on  endwall;  see  Figure  30,  for  example. 


TABLE  5 (Cont’d) 


TABLE  5 (Cont’d) 


30-degrees  (half-angle  of  'V');  HVG  = 0.1-in.;  on  floor;  see  Figure  31. 


Run  No.  6W  XWLE  XVGLE+  Po  T0  P1  ^1  qEl  ^Tl  X 10  COMMENTS 


egrees  (half-angle  of  V );  HVG  = 0.1-in;  on  floor;  see  Figure 


REFERENCE  MAINSTREAM  CONDITIONS  AT  STATION  1 FOR  CIRCULAR  PITOT  PROFILES  MEASURED  AT  MACH  2 
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(i)  No  sidewall  V.G. 

(ii)  No  flow  visualization. 
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LOCAL  SKIN  FRICTION  COEFFICIENTS  FROM  COBRA  PROBE  MEASUREMENTS  AT  MACH  2 


mir>fsjHO'0uno\O0NMOM0'00HONaMnff)NC00N 

«jNHHnn(nnHon>}HOOHOH>jMHOO 

HHHHHHHHHHHHHHHHHHrHHHrlH 

ooooooooooooooooooooooo 

ooooooooooooooooooooooo 

ooooooooooooooooooooooo 


HvOCOninminH^HNHinNONHMHvO^OH 

fnOO'^f'JMHHO'O'Mfn^O'ONONONOfOO^O'O' 

*H»HOO«H»HrH»HOOiHrHOOOOOrHrH»HOOO 

OOOOOOOOOOOOOOOOOOOOOOO 

OOOOOOOOOOOOOOOOOOOOOOO 

ooooooooooooooooooooooo 


or^oo^iHf-ir^\or^oo,HfH»HONOOiHcomcNioo*d' 

HOOHNHCNNNCONCOCN^rNf^HnCNNinON^in 

ONOOOOOOr^rsNOOONOCOOOOOOOOO^OONOOCOOO 


(NrHHHCNCNHHHHHNHH 


HrHHNHHHH 


ONO>moor^r^<ror^mor>»*d-omvomrHOvomoNO 

r^vovovor^r^r^r^vovooor^r^r^vovOv£>oooooor>*sor^ 

HHHrlHHHHHHHHHHHHrlHHHHHH 

OOOOOOOOOOOOOOOOOOOOOOO 

ooooooooooooooooooooooo 

ooooooooooooooooooooooo 


r^ON(y»'3’voomrHr^r^m'd,inr*^moovOiHC'")0<rrHvo 

vCHONO'J*^(n(NOO'^^HONHO'0'(N»}CnOHa 

»HrHOfHfHiH»H»H»HOiHrH»HOrHOOrHrHrHrHiHO 

OOOOOOOOOOOOOOOOOOOOOOO 

OOOOOOOOOOOOOOOOOOOOOOO 

OOOOOOOOOOOOOOOOOOOOOOO 


HcjNCNNO'srnHooocooorNHiriHONfOrNCNco^O' 
mC^OOCOCNCNrHOCOOOrHCNO>OOC^OOr^OCNi-HCX)C^r^ 
»H  O O O »H  rH  r— I i-H  O O »H  rH  O O O O O H H H O O O 
OOOOOOOOOOOOOOOOOOOOOOO 
OOOOOOOOOOOOOOOOOOOOOOO 

OOOOOOOOOOOOOOOOOOOOOOO 


co<rcMCMiHininooooococMcon-oor^vo<-v£>OrHvOiH 

c^fvjcovj-NncoOsjHoOsd-nnHoocNOcNHsrinm 

ONOOOOO'CONOOOOON(J'O^OM3NOO^a\OON^ON^ 


»HvDor^incooocM<rr^<r<rooinr^a>cocoooo><vocMCM 

coooinmrifNjHO'vD^Hcnr^^rs^^oocsaurir^'j 

MHHHMMf'JHHHNMHHHHHHNHHHH 

ooooooooooodoooddoooooo 


rH  rH  O O O »H  rH  t— H «— I O O O O O 

r^cNcooococM<^inOP-.o<rcoocMOin 


on  o>  o> 

in  CO  cm  rH 


OOOOrH»HOOOOrHOOOOOO 
III  I I I I I 


o o o o 

I I I I 


<r<r<j-<r<r<r^^-<r^-corororororororocMCMCN(MCsj 


o 

z 

rH  H CM 

3 

4 

5 
1 
2 

z 

I 

7 

€ 

3 

4 

5 

6 
1 
2 

.H  cm  co  <r  in 

O O 

O 

o 

O 

o 

c 

m vo 

r- 

oc 

ON 

o 

3 
a i 

<r 

<r 

•'tf 

m 

CM  CM 

CM 

CM 

CM 

CM 

0.012-in.  d = 0.012-in 


TABLE  7.3  (Cont’d) 


0.012-in.  d = 0.012-in 


REFERENCE  MAINSTREAM  CONDITIONS  AT  STATION  1 FOR  COBRA  PROFILES  MEASURED  AT  MACH  4 
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h IS  DISTANCE  PERPENDICULAR 
TO  WALL 


DEFINITION  OF  JET  PEAK  AND  WAKE  TROUGh 


FIG.  2(b):  SLOT  BLOWING  NOZZLE  ON  AEROFOIL  II 
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FIG.  3:  PRESSURES  AROUND  SUPERCRITICAL  AEROFOIL  SHOWING  EFFECT  OF 
SLOT  BLOWING  AND  JET  FLAP  ON  AEROFOIL  II  AT  MACH  0.85 
AND  CHORD  REYNOLDS  NUMBER  OF  30  MILLION 


WEDGE  AS  SHOCK  GENERATOR 


rrimii— in  ■iiunri  imrirwiT n 


V/ 


TWO  AND  THREE-DIMENSIONAL  SHOCK  WAVE/BOUNDARY-LAYER  INTERACTIONS  GENERATED  BY  A WEDGE 

IN  A RECTANGULAR  DUCT 


-99- 


= 1-47' 


✓Normal  shock 


/^Supersonic 
Unshaded  subsonic 
M = 0-7I  Eddy  flow 


Vortex  sheet 


Flow  field  through  the  interaction 


Separation 

poin^ 


Normal 
, shock 


I 


p = Static  pressure  on  plate 
pQ=  Free  stream  total  pressure 

x = Distance  downstream  measured 
from  start  of  interaction 

8 =■  Thickness  of  undisturbed 
boundary  layer  at  x = 0 


Surface  pressure  distribution 


FIG.  6:  POSSIBLE  INTERPRETATION  OF  SWEPT  NORMAL  SHOCK/TURBULENT 
BOUNDARY-LAYER  INTERACTION  WHEN  VIEWED  ALONG  SHOCK, 
AFTER  SEDDON,  REFERENCE  52 


- 100- 


(a)  THREE-DIMENSIONAL  SEPARATION 
LINE  AS  'ENVELOPE'  OF  UPSTREAM 
AND  DOWNSTREAM  SURFACE  SHEAR 
\ STRESS  TRAJECTORIES 


FIG.  7:  PLAUSIBLE  PATTERNS  OF  SURFACE  SHEAR  STRESS  TRAJECTORIES  FOR  A 
3D  SEPARATION  BENEATH  A SWEPT  SHOCK  WAVE 
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FIG.  9:  LIMITING  STREAMLINES  ON  UPSWEPT  FUSELAGE 


UNDISTURBED  TURBULENT  B/L  THICKNESS 


= 0-2in ; Rg  - l-5XI08ATMACH2 
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$ IS  MEASURED  WITH  RESPECT  TO  DEFLECTION  AT  B/L  EDGE 
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.]  PEAKE 


McCABE'S  THEORY,  APPROX. 
RANGE  FOR  THE  MACH 
NUMBERS  CONSIDERED 


AIR  SUPPLY  FROM  50,000  FT0  VARIABLE  DIFFUSER 

STORAGE  VESSELS  OF  5-FT.  TUNNEL  ON  WHEELS 


FIG.  12:  SCHEMATIC  ARRANGEMENT  OF  5 X 5-IN.  BLOWDOWN  WIND  TUNNEL 


FIG.  13:  TYPICAL  ASSEMBLY  OF  WORKING  SECTION  OF  5 X 5-IN.  WIND  TUNNEL  FOR  BLC 


1C  VIEW  OF  WEDGE  ASSEMBLY 


ED  BOUNDARY  QEFINBS  LIMITS  WITHIN  WHCH  TUNNEU  MAY  BE  OPERATED  WITH  SAFE  STRESS  LEVELS  IN  WEDGE  SUPPORTS 


FIG.  17:  TUNNEL  PERFORMANCE  CHART 
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FIG.  18(b):  PITOT  pROBE  UNIT  FOR  TRAVERSE  OF  2D  BOUNDARY  LAYER,  END  VIEW 


SECTION  IN  PLAN  VIEW 


0 002  DlA  UNSHIELDED  PLAN  VIEW 

THERMOCOUPLE  BEAD 


SIDE  ELEVATION 


FIG.  19:  HEAD  OF  COBRA  PROBE,  DESIGN  DIMENSIONS 


FRONTAL  VIEW  FROM  STARBOARD  SIDE  OF  PROBE  SHOWING 
THERMOCOUPLE  BEAD  ABOVE  CENTRE  PITOT  ORIFICE 


0278 


DIMENSIONED  FRONTAL  VIEW 


FIG.  20:  MAGNIFIED  VIEW  OF  COBRA  PROBE  HEAD  AS  MANUFACTURED 
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STAGNATION  PRESSURE 
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MOTOR  GENERATOR 
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TWO  STATIC  PRESSURES 
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FIG.  22:  FLOOR  INSERT  TO  SUPPORT  COBRA  PROBE 


LAST  CHARACTERISTIC 


FIG.  25:  ISOMETRIC  VIEW  OF  SUPERSONIC  BLOWING  NOZZLE 
IN  FLOOR  OF  5 X 5-IN.  WIND  TUNNEL 


(SEE  TABLE  2 ) 


FIG.  26:  INJECTION  NOZZLE  PROFILE 
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FIG.  27:  CALCULATED  RANGE  OF  AVAILABLE  JET  EXCESS  MOMENTUM  FOR  CORRECTLY  EXPANDED  JET  FLOW 
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FIG.  30:  EFFECT  OF  SIDEWAl  L VORTEX  GENERATOR  ON  EMPTY  TUNNEL  FLOW  AT  MACH  2 
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FIG.  31:  SURFACE  PRESSURE  DISTRIBUTIONS  AND  FLOW  VISUALIZATION  IN 

EMPTY  TUNNEL  AT  MACH  4 
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REPEATABILITY  OF  VELOCITY  PROFILES  IN  EMPTY  TUNNEL  AT  MACH  2 MEASURED  WITH  CIRCULAR  PITOT 


FIG.  35:  SPANWISE  SYMMETRY:  VELOCITY  PROFILES  IN  EMPTY  TUNNEL  AT  MACH  4 MEASURED  WITH  CIRCULAR  PITOT 

WITHOUT  AND  WITH  FLOOR  VORTEX  GENERATOR 


FIG.  36:  SEMI-LOGARITHMIC  VELOCITY  PROFILES  IN  EMPTY  TUNNEL  FROM 
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FIG.  38:  COBRA  AND  CIRCULAR  PITOT  TUBE  VELOCITY  PROFILES  ALONG  <t  AT  MACH  2 IN  EMPTY  TUNNEL 
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FIG.  40:  LOCAL  SKIN  FRICTION  FOR  FLOW  IN  EMPTY  TUNNEL  AT  MACH  2 
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FIG.  42:  LOCAL  SKIN  FRICTION  FOR  FLOW  IN  EMPTY  TUNNEL  AT  MACH  4 
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FIG.  49:  SURFACE  PRESSURE  DISTRIBUTIONS  AND  FLOW  VISUALIZATION  FOR  5 

AND  11.5-DEGREES  AT  MACH  2 
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FIG.  51:  SCHLIEREN  AND  SURFACE  FLOW  VISUALIZATION  AT  SMALL  WEDGE  ANGLES 

AT  MACH  4 (FLOOR  V.G.  IN) 
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FIG.  53:  SURFACE  PRESSURE  DISTRIBUTIONS  AND  FLOW  VISUALIZATION  FOR 
6W  = 10  AND  14-DEGREES  AT  MACH  4 (FLOOR  V.G.  IN) 
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FIG.  54:  SURFACE  PRESSURE  DISTRIBUTIONS  AND  FLOW  VISUALIZATION  FOR  STRONG 
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FIG.  55:  SURFACE  FLOW  VISUALIZATION  WITH  CHANGE  OF  WEDGE  AXIAL  LOCATION  FOR 

= 8-DEGREES  AT  MACH  2 
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FIG.  61:  PITOT,  TOTAL  TEMPERATURE  AND  YAW  ANGLE  DISTRIBUTIONS  FOR 
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CROSSFLOW  AND  STREAMWISE  FLOW  VELOCITY  PROFILES  FOR  Sw  = 8-DEGREES  AT  MACH  2 


FIG.  69:  MAP  OF  COBRA  PROBE  TRAVERSE  LOCATIONS  AT  MACH  4,  5W  = 8-DEGREES 
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FIG.  74.  INTEGRAL  PARAMETERS  FOR  6W  = 8-DEGREES  AT  MACH  4 
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FIG.  84:  POLAR  PLOTS  OF  VELOCITY  PROFILES  FOR  6W  = 16-DEGREES  AT  MACH  4 
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FIG.  85:  INTEGRAL  PARAMETERS  FOR  5W  = 16-DEGREES  AT  MACH  4 
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FIG.  87:  LOCAL  SKIN  FRICTION  AND  FLOW  DEFLECTIONS  FOR 
= 16-DEGREES  AT  MACH  4 
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FIG.  95(b):  PRESSURES  ALONG  TEST  WALL  AT  MACH  2 AND  WEDGE  ANGLE  OF  8-DEGREES  AT  A FIXED  WEDGE  POSITION, 

XWLE  = 1.5-IN.,  WITH  WALL  JET  ON:  AT  0.  = - 12-DEGREES 
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FIG.  97(b):  PRESSURES  ALONG  TEST  WALL  AT  MACH  2 AND  WEDGE  ANGLE  OF 
11.5-DEGREES  AT  A FIXED  WEDGE  POSITION,  XWLE  = 1.5-IN., 

WITH  WALL  JET  ON:  AT  0,  = -20-DEGREES 
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FIG.  97(c):  PRESSURES  ALONG  TEST  WALL  AT  MACH  2 AND  WEDGE  ANGLE  OF 
11.5-DEGREES  AT  A FIXED  WEDGE  POSITION,  XWLE  = 1.5-IN., 

WITH  WALL  JET  ON:  AT  0=-  12-DEGREES 
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FIG.  97(d):  PRESSURES  ALONG  TEST  WALL  AT  MACH  2 AND  WEDGE  ANGLE  OF 
11.5-DEGREES  AT  A FIXED  WEDGE  POSITION,  XWLE  = 1.5-IN., 

WITH  WALL  JET  ON:  AT  0:  = 0-DEGREE 
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FIG.  98:  AIR  INJECTION:  SURFACE  FLOW  VISUALIZATION  TO  DEMONSTRATE  INCREASING 
STAGNATION  PRESSURE  OF  JET  FLOW  WITH  Oi  = O-DEGREE  AND  Mj  s 3:  IN  EMPTY 
TUNNEL,  AND  AT  Sw  = 11.5-DEGREES  (XWLE  = 3.0-IN.) 
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FIG.  99:  PRESSURES  ALONG  TEST  WALL  AT  MACH  2 AND  WEDGE  ANGLE  OF  11.5-DEGREES 
AT  A FIXED  WEDGE  POSITION,  XWLE  = 3.0-IN.,  WITH  WALL  JET  ON:  0t  = O-DEGREE 
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FIG.  101:  AIR  INJECTION:  SURFACE  PRESSURES  AND  FLOW  VISUALIZATION  THROUGH 
2D  SHOCK  BOUNDARY-LAYER  INTERACTION  FOR  5W  = 11.5-DEGREES  AT  MACH  2, 

WITH  M,  as  3,  AND  6{  = 0-DEGREE 
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APPENDIX  X 

A RESUME  OF  THE  THREE-DIMENSIONAL,  COMPRESSIBLE 
TURBULENT  BOUNDARY-LAYER  CALCULATION  METHOD  OF  NASH 


As  the  3D  compressible  turbulent  boundary-layer  calculation  method  of  Nash  has  not  been 
published  in  full,  we  have  taken  the  liberty  of  excerpting  the  relevant  equations,  especially  those  treating 
the  Reynolds  stresses,  from  a NASA  Contract  Report1 15  that  has  seen  only  limited  distribution.  We 
have  lollowed  the  notation  of  Nash  in  this  appendix  to  facilitate  reference  to  his  3D  turbulent  boundary- 
layer  calculation  method  in  incompressible  flow121 . 

Briefly,  Nash’s  method  provides  the  capability  of  calculating  profiles  of  mean  velocity  and 
turbulent  shear  stress  in  a fully  three-dimensional  flow  field  over  a generally  curved  surface,  in  arbitrary 
pressure  gradients  and  with  no  restrictions  to  small  crossflows.  The  method  involves  the  numerical  in- 
tegration of  the  equations  of  motion  and  the  continuity  equation,  together  with  an  empirical  equation 
for  the  E.eynolds  stresses,  that  is  based  on  the  turbulent  kinetic  energy  equation.  The  extension  of  the 
incompressible  method  to  treat  compressible  flows  is  made  with  the  assistance  of  the  Crocco  relation 
for  temperature  (rather  than  including  the  energy  equation)  which  should  be  valid  for  Mach  numbers 
at  least  up  to  two  in  the  presence  of  an  adiabatic  wall. 

In  the  formulation  of  the  programme  supplied  to  us  by  Nash,  the  general  curvilinear  co-ordinate 
system  (x,  y,  z)  is  set  up  to  calculate  the  flow  over  a flat  plate  swept  wing  with  prescribed  pressure  dis- 
tribution, where  the  wing  root  and  tip  are  stream  surfaces  parallel  to  the  undisturbed  mainstream  at 
infinity.  The  plane  y = 0 represents  the  surface  of  the  wing,  the  x-  and  z-axes  are  respectively  parallel 
and  perpendicular  to  the  approaching  mainstream  flow,  and  the  swept  planform  is  situated  in  the  posi- 
tive x-z  quadrant*.  In  view  of  the  tunnel  geometry  in  the  swept  interaction  experiment,  we  adopted 
the  present  formulation  and  simply  ‘unswept’  the  wing.  The  plane  x = 0 is  far  ur"  cream  in  the  nominal 
undisturbed  boundary  layer,  while  x = 1 is  situated  somewhere  near  the  wedge  trailing  edge.  The  nor- 
malizing distance  for  x is  that  absolute  length  between  the  i tart  and  *>nd  points  of  the  calculation, 
parallel  to  the  tunnel  centre-line  (equivalent  to  the  stream  wise  chord  dimension  of  the  swept  wing 
flow).  The  test  case  chosen  was  with  the  deflected  wedge  angle,  5W  = 8-degrees  at  Mach  2,  and  we 
situated  the  deflected  wedge  leading-edge  in  the  plane  of  the  ‘wing  root’.  However,  since  the  boundary 
condition  at  the  wing  root  prescribed  a flow  parallel  to  the  oncoming  mainstream,  the  deflected  wedge 
surface  was  neglected,  and  the  measured  shock  pressure  rise  pi  us  the  upstream  undisturbed  momentum 
deficit  thicknesses  were  the  essential  inputs  to  the  programme.  The  mean  velocity  and  shear  stress 
profiles  at  the  commencement  of  the  calculation  were  generated  in  the  computer  programme. 

The  symbols  required  to  follow  the  general  formulation  are  as  follows:  these  are  used  only 
in  this  appendix,  not  in  the  main  report. 

at  empirical  constant,  identified  with  the  ratio  of  the  resultant  shear  stress  to  the  resultant 

turbulent  kinetic  energy,  see  Equation  (A-1.9) 

empirical  function  of  (y/S),  a diffusion  parameter,  see  Equation  (A-1.9) 
metric  coefficients 
curvature  parameters 
dissipation  length;  function  of  (y/6) 


* 


a2 

h,,h3 
K13>  K3l 

L 


N _>te  that  the  co-ordinates  (x,  y,  z)  in  the  calculation  correspond  with  (x,  h,  y)  in  the  experiment, 
when  the  leading-edge  is  unswept. 
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M 


Mjch  number 


P 

Q 

s 

U,  V,  w 
u,  v,  w 
x,  y,  z 


r 

8 

8* 

P 

(b  <b 


7Wx>  rWz 


static  pressure 

resultant  mean  velocity,  Q2  = U2  + W2,  to  the  boundary-layer  approximation 
length  of  line  in  space 

mean  velocity  components  in  the  x,  y,  z directions 
fluctuation  velocity  components,  corresponding  with  U,  V,  W 

orthogonal  curvilinear  co-ordinates,  with  y measured  normal  to  wing  surface.  U,  V,  W 
are  in  the  respective  directions:  x,  y,  z 

constant  in  Equations  (A-1.15) 

boundary-layer  thickness 

displacement  thickness,  Equation  (A-1.16) 

air  density 

functions  in  Equations  (A-1.9)  and  (A-1.10) 
magnitude  of  shear  stress  vector  within  boundary  layer 
skin  friction  components  in  the  x,  z directions 


Subscripts 

e outer  edge  of  boundary  layer 

“ freestream  at  infinity 

We  choose  an  orthogonal  curvilinear  co-ordinate  system  (x,  y,  z)  with  the  wing  surface  at 
y = 0.  h,  and  h3  are  metric  coefficients  for  distances  measured  in  the  x and  z directions,  such  that 
the  length  of  line  element,  ds,  is 


(ds)2  = h?(dx)2  + (dy)2  + hf(dz)2  (A-l.l) 


The  curvature  parameter  K,  3 , K3 , are  defined  by 


1 3h,  ! 3h3 

K'3  hj h3  0z  ’^31  “ h,h3  9x  (A-1.2) 


and  represent  the  geodesic  curvatures  of  the  lines  x = constant,  z = constant.  The  metric  coefficient, 
h2,  associated  with  distances  measured  in  the  y direction  is  chosen  to  be  unity,  with  the  result  that 
the  curvature  parameters  K2 , , K23  are  identically  zero' 2 ' . 
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u a v a w a r 3U 

h,  S <w)  * £ 57  (uv)  55  m * (K'*U-  K»W)  ™ * 2a>  p2 a™2>’'2  57 


a uv  

+ 4>x  + — (a2uv)  + — (uv2  +vw2)*/4  = 0 

o y L 


(A-1.9) 


u a % v a w a r 3W 

(VW)  + £ * (vw)  + £Jfr  (^)MK31W-K13U)uv  + 2a1^(uv’-+vw2)2/2  — 


a vw 

+ 4>z  + — (a2vw)  + — (uv2  + vw2)*/4  = 0 


(A-1.10) 


where  <l>x , 4>z  have  the  property  that 


uv  <l>x  + vw  4>z  = 0 


(A-l.ll) 


If  Equations  (A-1.9)  and  A-1.10)  are  combined  into  one  for  the  magnitude  of  the  turbulent 
shear  stress  vector, 


(uv2  + vw2)1/2  = r/p 


(A-1.12) 


the  resulting  equation  is  of  the  same  form  as  the  empirically  modified  turbulent  kinetic-energy  equation 
used  by  Bradshaw  et  al1 22 , which  is 


u 3t  v 3t  w a t r au  aw  a r3^2  ~\ 

h,  55  * i7  55 + £ 55 + 2a>|/uv  57  * p™  57  * 55  (ajT)  + THlJ  ■ 0 (A-113) 


On  this  basis,  a]  is  taken  to  be  the  same  constant,  and  a2  and  L the  same  functions  of  y/6  as  in  the  2D 
incompressible  case;  a,  is  identified  with  the  ratio  of  the  resultant  shear  stress  to  the  intensity  of  the 

turbulence,  pq2 , L is  identified  with  the  dissipation  length  and  a2  is  a diffusion  parameter. 

The  terms 

2a,  (uv2  + vw2)1/2  ^ , 2a,  (uv2  + vw2)>/2  — 

3y  3y 

are  assumed  to  represent  the  combined  effect  of  production  and  the  part  of  the  pressure-strain  corre- 
lation which  opposes  it,  and  the  terms 


uv  (uv2  + vw2 )’ /4  , vw  (uv2  + vw2 )* /4 

L L 


are  identified  with  the  remaining  part. 
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By  virtue  of  Equation  (A-l.ll),  the  functions  4>x , 4>z  do  not  contribute  to  the  magnitude  of 
the  shear  stress  vector  but  determine  its  directional  response.  The  method  of  Nash  can  be  regarded  as 
corresponding  with  a particular  choice  of  4>x , 4>z , which  leads  to  a rapid  direction  response  of  the  shear 
stress  vector  to  a rotation  of  the  rate-of-strain  vector,  such  that  the  two  vectors  are  everywhere  co- 
directional 


_/dU  /d  W 

/ ay  v 7 ay 


(A-1.14) 


The  essential  difference  between  the  approaches  of  Bradshaw123  and  Nash  is  that  Bradshaw  makes 
4>x  = <J>Z  = 0,  which  yields  a relatively  slow  direction  response. 

4>x , <l>z  are  given  by 

5L  ..|(K)'  .(=)■]'"] 


57  * ™{(^)  ’ (It)  I ] 


(A-1.15) 


in  which  T is  some  large  number  (from  experience,  any  number  greater  than  about  four). 

Equations  (A-1.9),  (A-1.10),  and  (A-1.15)  are  identical  with  their  incompressible  counter- 
parts. It  is  assumed  that  they  can  be  carried  over  to  compressible  flows  without  modification,  and 
that  the  empirical  constant  a,  (=0.15)  and  the  empirical  functions  a 2 and  L are  the  same  as  in  incom- 
pressible flow  (see  Ref.  121).  There  is  no  direct  experimental  support  for  these  assumptions  in  com- 
pressible 3D  flow.  However,  Nash  claims1 15  that  they  are  probably  valid  at  least  up  to  Mach  2.  In  2D 
flow,  they  appear  approximately  valid  up  to  higher  Mach  numbers,  for  comparisons  between  calculations 
performed  by  Nash  using  the  Nash  method  with  the  experimental  data  of  Peake  et  al34  at  Mach  4, 
showed  reasonably  good  agreement. 

The  method  of  solution  is  relatively  straightforward,  if  complex  in  programming  detail. 
Equations  (A-1.3),  (A-1.4),  (A-1.5),  (A-1.9)  and  (A-1.10)  form  a hyperbolic  set,  and  are  integrated  in 
a 3D  domain  (x,  y,  z)  by  an  explicit  finite  difference  method  based  on  a one-step  staggered  mesh 
scheme  (see  Ref.  124).  The  calculation  proceeds  in  the  x-direction,  and  values  of  the  five  dependent 
variables:  U,  V,  W,  uv,  vw  are  computed  on  successive  surfaces  x = constant.  Twenty  collocation 
points  are  typically  used  in  the  y-direction  (normal  to  the  wing  surface)  and  six  or  seven  in  the  z (or 
spanwise)  direction. 

The  numerical  solution  is  matched  at  y = 0.05  5 (approximately)  to  a separate  solution  for 
the  law  of  the  wall  region.  In  the  inner  layer,  it  is  assumed  that  the  resultant  velocity  Q,  and  the  skin 
friction,  rw , mutually  satisfy  the  law  of  the  wall.  This  assumption  then  provides  an  equation  for  tw 
in  terms  of  Q at  the  matching  station.  The  direction  of  the  skin  friction  vector  is  determined  from  an 
extrapolation  in  the  polar  plane  (U,  W)  making  use  of  the  correct  limiting  value  of  32W/9U2  as  U -*■  0 
(or  in  other  words,  from  the  pressure  gradient). 

The  boundary  conditions  at  the  wing  surface  (y  = 0)  and  at  the  outer  surface  of  the  integration 
domain  (at  y = 1.25  5 approximately)  are  handled  as  in  Reference  121.  The  boundary  conditions  on 
the  sides  of  the  integration  domain  call  for  special  treatment.  In  the  programme  supplied  by  Nash, 
these  were  stream  surfaces  at  z = constant,  parallel  to  the  undisturbed  freestream  velocity  at  infinity. 


— " ■ ,J 
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On  Displacement  Surface  Thickness,  6* 

In  a fully  3D  boundary  layer,  6*  cannot  be  related  directly  to  an  integral  of  the  velocity 
profile  at  one  position  on  the  body  surface1 17>118.  Instead  it  can  be  determined  from  the  normal 
component  of  velocity  at  the  boundary-layer  edge  by  means  of  the  partial  differential  equation 


Ue  05*  We  05* 

h,  0x  h3  0z 


{ V-  (y  - 5*) 


y > s 


(A-1.16) 


Once  the  velocity  field  is  computed,  this  equation  can  be  integrated  to  give  5*.  We  see  that  Equation 
(A-1.16)  is  compatible  with  the  2D  definition  of  6* 


5*  = /“  ^1 

either  where  the  velocity  profiles  are  collateral  (W  proportional  to  U)  or  where  deviations  with  respect 
to  x or  z vanish,  as  on  an  ‘infinite’  swept  wing. 

The  appropriate  boundary  conditions  required  to  integrate  Equation  (A-1.16)  are  provided 
by  the  conditions  at  the  ‘wing  root’  and  ‘wing  tip’  planes  of  symmetry,  where  05*/0z  = 0,  and  by  the 
initial  velocity  profiles. 


PQ  \ 
PeQc/ 


dy 


(A-1.17) 


